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FOREWORD 


This  technical  report  covers  an  investigation  on  unsteady  supersonic 
flutter  carried  out,  starting  on  January  1,  1974  and  ending  on  December 
31,  1977.  The  research  was  sponsored  by  the  Air  Force  Office  of  Scientific 
Research  under  Contract  F 44620-74-C-0040. 

For  most  of  the  contract  period,  the  research  was  carried  out  at  the 
Contractor's  site.  General  Electric  Corporate  Research  and  Development, 
Schenectady,  New  York;  due  to  the  unexpected  departure  of  the  principal 
investigator,  M.  Kurosaka,  from  the  General  Electric  Company  to  the  Univer- 
sity of  Tennessee  Space  Institute  on  September  1,  1977,  it  was  then  sub- 
contracted to  UTS I,  where  it  was  continued  and  completed. 

Aiding  the  principal  investigator,  I,  H.  Edelfelt  of  the  General 
Electric  CRD  performed  most  of  the  computational  task  and  C.  E.  Danforth, 
then  chief  consulting  engineer  of  aeromechanics  of  G.  E.  Aircraft  Engine 
Business  Group,  Evendale,  Ohio,  served  as  a consultant  to  the  program; 

J.  Q.  Chu  of  UTSI,  a graduate  research  assistant,  also  helped  the 
principal  investigator. 

The  contract  was  monitored  by  Lt.  Colonel  R.  C.  Smith,  Program 
Manager,  Directorate  of  Aerospace  Sciences,  Air  Force  Office  of  Scientific 
Research,  United  States  Air  Force,  Bolling  Air  Force  Base,  D.  C. 

This  final  report  supersedes  all  the  previous  interim  reports 
and  includes  all  the  technical  papers,  both  published  and  under  preparation, 
which  were  written  under  the  research  contract. 
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OBJECTIVES  OF  RESEARCH 

Overall  objective  was  to  develop  an  analytical  tool  which  enables 
one  to  define  instability  boundaries  for  supersonic  unstalled  flutter  of 
aircraft  engine  fans  and  compressors  and  also  to  provide  insights  into  its 

prevention.  t 

! 

The  region  of  the  supersonic  instability  often  falls  on  the  operating 

i 

line  of  a high  speed  fan;  consequently,  the  existence  of  such  flutter  pre-  I 

sents  itself  as  a barrier  problem  for  the  design  of  high  Mach  number  fans.  ? 

In  order  to  overcome  this  serious  obstacle,  it  is  necessary  to  improve  at 
an  accelerated  pace  our  capability  to  define  and,  at  the  same  time,  to 
minimize  the  region  of  blade  instability. 

The  present  program  was  comprised  of  three  phases.  In  Phase  1,  a 
simplified  model  of  flat-plate  airfoils  was  adopted  and  its  objective 
was  directed  to  extending  our  previous  low  frequency  analysis  to  a higher 
frequency  of  practical  interest  so  as  to  predict  oscillatory  forces 
acting  on  unsteady  supersonic  airfoils  in  a cascade  subject  to  subsonic 
axial  velocity.  Phase  2 was  conducted  so  as  to  remove  the  restriction  of 
flat-plate  airfoils;  this  extension  is  necessitated  because  there  are 
substantial  experimental  data  which  indicate  the  Importance  of  airfoil 
shape.  Phase  3 addressed  the  effect  of  flow  three-dimensionality, 
particularly  the  influence  of  radial  velocity  gradient  in  steady  base 
flow  and  the  means  of  flutter  suppression  by  the  provision  of  liners  on 
outer  casing  walls. 
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SVM4ARX  OF  SIGNIFICANT  RESULTS 


Since  all  seven  technical  papers,  which  are  either  already 
published  or  under  preparation  for  publication,  are  attached  in  this 
report,  we  relegate  the  necessary  details  to  them  and  here  summarize 
only  the  important  results  and  conclusions. 

In  Phase  1,  airfoils  were  assumed  to  be  flat  plates.  An  expres- 
sion of  unsteady  pressure  distribution  valid  for  the  range  of  frequency 
parameter  of  practical  interest  was  derived  in  closed  form  which  is 
applicable  to  any  cascade  geometry  and  arbitrary  motion  of  airfoils; 
this  was  an  extension  of  our  previous  low  frequency  analysis  (ASME 
Transactions,  Journal  of  Engineering  for  Power,  1974,  Vol.  95,  January, 
pp.  13-31).  When  the  unsteady  pressure  distribution  was  used  as  input 
to  flutter  prediction,  the  following  two  major  conclusions  were  establish- 
ed: the  zone  of  torsional  instability  tends  to  shrink  as  the  frequency 
parameter  increases,  and  at  the  value  of  frequency  parameter  of  convention- 
al design,  the  bending  instability  is  predicted,  the  latter  being  in 
agreement  with  the  experimental  data.  The  details  of  these  were 
presented  at  the  IUTAM  Symposium  on  Aeroelasticity  in  Turbomachinery , 
Paris,  October,  1976,  and  published  in  Revue  Francaise  de  Mecanique. 

Numero  Special,  1976,  pp.  57-64,  appended  here  as  Appendix  1.  In 
addition,  our  discussion  on  the  issue  of  'resonance'  in  a supersonic 
cascade  has  been  published  in  the  AIAA  Journal,  vol.  13,  No.  11, 


November,  1975,  pp.  1514-1516  (attached  as  Appendix  2). 
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In  Phase  2,  the  restriction  of  flat  plate  was  removed.  The  effect 
of  airfoil  shape  was  examined  first  for  an  isolated  airfoil  and  the 
cumulative  effects  of  nonlinearity  on  the  unsteady  pressure  in  the  far 
field  was  studied.  For  an  oscillating  airfoil  whose  contour  is  of  para- 
bolic arcs,  an  analytical  expression  of  unsteady  flow  field  was  obtained. 

The  results  clearly  indicate  that  the  effect  of  airfoil  shape  deeply 
affects  the  unsteady  flow  in  the  far  field.  This  was  published  in  the 
Journal  of  Fluid  Mechanics,  vol.  83,  part  4,  pp.  751-773,  1977, appearing  here 
as  Appendix  3.  The  obvious  implication  of  this  result  for  a single  airfoil 
was  indeed  borne  out  in  the  subsequent  extension  to  the  cascaded  airfoils 
where  the  importance  of  the  airfoils  shape  was  in  fact  confirmed,  the 
details  of  which  are  described  in  Appendix  4. 

In  Phase  3 the  effect  of  the  flow  three-dimensionality  on  flutter 
boundary  was  examined.  The  objective  here  was  not  so  much  on  the  precise 
prediction  of  the  three-dimensional  unsteady  flow  field,  which  would 
amount  to  an  extremely  complicated  task;  rather  it  was  deemed  more  bene- 
ficial to  focus  attention  on  the  effect  of  the  surrounding  casing  walls 
and  specifically  to  examine  whether  the  use  of  sound  absorbing  material  for 
the  wall  might  significantly  relieve  the  fluctuation  of  pressure  on  the 
airfoils.  Since  the  critical  section  of  the  fluttering  airfoils  in  super- 
sonic flutter  is  located  near  the  very  tip  of  the  airfoils,  the  provision 
of  the  acoustically  treated  wall  on  the  adjacent  outer  casing  would  be 
directly  effective  in  relieving  the  unsteady  fluctuation  and  eventually 
reduce  the  level  of  unsteady  pressure  acting  on  the  airfoil  surface. 

Based  on  this  idea,  an  analysis  was  carried  out  for  a model  problem 


i 
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where  an  isolated,  three-dimensional  airfoil  is  placed  in  a duct 
whose  upper  wall  is  acoustically  treated.  The  flow  in  the  duct  is  super- 
sonic and  initially  assumed  to  be  uniform;  the  airfoil  is  oscillating  in 
the  transverse  direction.  The  solution  indicates,  as  expected,  the  un- 
steady pressure  field  produced  by  the  motion  of  the  airfoil  is  significantly 
affected  by  the  presence  of  the  acoustically  treated  "soft"  wall,  as 
described  in  Appendix  5.  During  this  phase,  an  additional  analysis  was 
made  to  examine  the  effect  of  the  spanwise  variation  of  the  incoming 
steady  flow.  This  was  done  in  order  to  assess  the  refraction  effect  of 
the  acoustic  signal  due  to  its  interaction  with  the  non-uniform 
steady  velocity  field  where  the  fans/compressors  bladings  are  immersed. 

This  presented  a challenging  task  of  investigating  the  nature  of  the 
acoustic  wave  propagating  through  non-uniform  flow,  which  is  mathematically 
described  by  governing  equations  with  varying  coefficients.  From  our 
effort  to  surmount  this  difficulty  evolved  a powerful  but  simple  technique 
which  was  found  to  be  highly  effective  in  solving  a general  class  of  problems 
involving  unsteady  disturbances  propagating  through  non-uniform  media; 
this  is  reported  in  Appendix  6.  This  general  technique  was  then  applied 
to  examine  specifically  the  effect  of  velocity  gradient  upon  the  unsteady 
pressure  on  supersonic  airfoils,  where,  as  shown  in  Appendix  7.  the 
three-dimensional  effect  was  quantified. 
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APPENDIX  1 


"Some  Recent  Developments  in  Unsteady  Aerodynamics  of  a Supersonic  Cascade . " 
presented  at  the  Symposium  of  International  Union  of  Theoretical  and  Applied 
Mechanics  on  Aero  elasticity  in  Turbomachines;  published  in  Revue  Prancaise  de 
Mecanique 3 Numero  Special 3 19?6y  pp.  57-64. 


SOME  RECENT  DEVELOPMENTS  IN  UNSTEADY  AERODYNAMICS 
OF  A SUPERSONIC  CASCADE 

by 

M.  KUROSAKA  and  I.  H.  EDELFELT 


1.  INTRODUCTION 

Quite  unexpectedly,  the  advent  of  high  speed  turbo-fan  aircraft  engines  unfolded  a generlcally 
new  type  of  blading  Instability  — supersonic  installed  flutter.  The  Instability,  which  until  re- 
cently had  been  virtually  unheard  of,  is  found  to  occur  when  the  tip  speed  exceeds  sonic  velocity; 
since  the  unstable  regions  spread  over  the  operating  line  where  the  Incidence  Is  small,  the  flow 
over  the  bladings  Is  not  stalled.  In  order  to  develop  a consistently  reliable  means  of  averting 
this  trouble,  the  unsteady  aerodynamic  forces  acting  on  the  oscillating  airfoils  In  a fan  must  be 
accurately  predicted.  In  the  present  paper,  we  shall  attempt  to  present  recent  developments  of  the 
unsteady  aerodynamics  of  supersonic  cascades  carried  out  for  the  past  several  years  at  the  Research 
and  Development  Center,  General  Electric  Company;  this  Is  essentially  an  extension  of  our  previous 
work.  [lj.  The  other  related  works  In  this  area  which  have  appeared  until  now  are  [2]  to  [10], 

While  It  will  be  attempted  to  make  the  paper  as  self-contained  as  possible  starting  from  the  descrip- 
tion of  the  relevant  background,  prevailing  emphasis  will  be  focused  on  the  recent  and  hitherto  un- 
published analytical  results,  which  we  hope  to  be  of  both  sufficient  Interest  and  practical  importance. 

2.  CHARACTERISTICS  AND  PHYSICAL  MECHANISM  OF  SUPERSONIC  UNSTALLED  FLUTTER 

When  In  flutter,  the  bladings  vibrate  sinusoidally  at  their  natural  frequencies;  both  predomi- 
nantly bending  (or  flexural)  and  torsional  motion  of  the  airfoil  tip  section  have  been  observed  In 
unstalled  supersonic  flutter.  The  contour  of  the  airfoil  Is  known  to  Introduce  a significant,  first- 
order  effect  on  the  flutter  boundary.  As  an  example,  for  an  airfoil  which  had  experienced  unstalled 
supersonic  flutter,  apparently  minor  modifications  In  the  airfoil  shape  succeeded  In  removing  the 
Instability.  The  back  pressure  also  affects  the  unstable  region  but  lowering  the  back  pressure  It- 
self does  not  cure  the  flutter  problem. 

As  regards  the  physical  mechanism  of  the  unstalled  supersonic  flutter,  the  first  obvious  cause 
that  comes  to  mind  Is  the  self-sustained  shock-boundary  layer  oscillation.  The  shocks  Impinging  on 
the  boundary  layer  often  cause  Its  separation.  For  an  Isolated  airfoil  In  the  transonic  range  where 
a shock  appears  over  the  airfoil  surface,  such  shock-boundary  layer  Interaction  has  often  asserted 
to  be  capable  of  sustaining  self-induced  oscillation.  If  this  Is  Indeed  true,  similar  phenomenon 
can  certainly  be  expected  even  in  a cascade.  In  fact,  visualization  studies  of  the  fluttering  air- 
foils In  a supersonic  cascade  apparently  show  such  oscillatory  movement  of  the  shock-boundary  layer. 
Since  the  shock-boundary  oscillation  Is  the  phenomenon  most  conspicuously  visible  In  such  a study, 
one  Is  Irreslstably  tempted  to  conclude  that  this  Is  the  driving  mechanism  of  the  supersonic  unstalled 
flutter,  where  the  oscillatory  movement  of  large  pressure  rise  associated  with  the  Impingement  of  the 
shock  may  be  directly  responsible.  However,  before  one  hastily  jumps  to  the  conclusion,  we  have  to 
address  the  question  of  whether  the  shock-boundary  layer  oscillation  Is  the  real  cause  of  the  flutter 
or  the  aftereffect  of  the  flutter,  which  Itself  Is  Initiated  by  some  other  mechanism.  Indeed,  a 
careful  review  of  the  Issue  leads  to  what  appears  to  be  an  Inescapable  conclusion  — the  self- 
sustained  oscillation  between  the  shock  and  the  boundary  layer  does  not  play  a major  role  In  the 
present  Instability.  For  example.  In  (11],  Llepmann  and  Ashkenas  found  that,  although  the  oscilla- 
tion of  a shock  was  observed  for  a transonic  airfoil  placed  In  a wind  tunnel,  the  removal  of  the 
sources  of  disturbance,  which  otherwise  always  exist  In  any  wind  tunnel,  virtually  eliminated  the 
oscillation  of  the  shock.  Other  available  experimental  evidences  tend  to  support  the  similar  con- 
clusion that  the  self-sustained  oscillation  between  the  shock  and  the  boundary  layer  Is  a weak 
fluctuation. 

Then,  what  Is  the  predomlnent  cause  of  the  supersonic  unstalled  flutter?  It  Is  well-known  that 
for  an  Isolated  airfoil  In  supersonic  flow  where  the  shocks  are  attached  to  both  the  leading  and 
trailing  edges  and  do  not  appear  on  the  airfoil  surface,  the  motion  of  the  airfoil  could  become  un- 
stable under  certain  circumstances.  Such  an  Instability,  which  occurs  In  the  obvious  absence  of 
self-sustained  shock-boundary  layer  oscillation.  Is  caused  by  the  phased  lag  response  of  the  unsteady 
flow  outside  of  the  boundary  layer  to  the  motion  of  the  airfoil.  At  the  debut  of  supersonic  un- 
stalled flutter  In  the  late  sixties,  we  Identified  this  Invlscld  mechanism  as  the  primary  cause  of 
the  Instability  and  the  validity  of  the  diagnosis  has  been  Increasingly  buttressed  by  supporting 
evidences. 
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3.  DIFFICULTIES  ASSOCIATED  WITS  THE  AHALISIS 


For  the  reason  just  stated.  In  the  analytical  formulation,  we  assume  the  flow  to  be  Invlscld. 

Also  the  flutter  tends  to  become  more  severe  as  the  pressure  ratio  Is  lowered.  Hence  we  examine  the 
situation  of  the  low  pressure  ratio  where  the  flow  can  be  approximately  taken  to  be  purely  supersonic 
from  Inlet  to  discharge.  In  addition.  In  this  first  attempt,  the  flow  Is  assumed  to  be  two-dimensional 
Even  after  these  Idealization,  the  analysis  of  the  problem  on  hand  Is  a difficult  one.  One  source  of 
the  difficulties  results  from  the  fact  that  In  a current  design  of  the  fan,  the  axial  velocity  Is  sub- 
sonic even  when  the  relative  velocity  Is  supersonic.  This  axial  subsonic  velocity  Implies  that  the 
Mach  waves  emanating  from  the  leading  edges  extend  upstream  of  the  cascade  front  line.  Figure  (1), 
and  there.  In  what  Is  called  preinterference  zone,  one  has  to  take  the  mutual  Interference  Into  ac- 
count. Here  the  axial  velocity  supersonic,  the  Mach  waves  would  be  confined  within  the  blade  passage 
and  consequently  there  would  be  no  region  of  mutual  Interference  upstream  of  the  cascade  front  line. 
Furthermore,  If  one  looks  at  the  wake  region  off  the  airfoil.  Figure  (1),  we  note  that  the  wake  velo- 
city between  Q and  R does  Influence  the  trailing  edge  section  of  the  adjacent  blade  between  points  S 
and  T.  Since  one  does  not  & priori  know  the  wake  velocity.  It  must  be  sought  as  a part  of  the  solu- 
tion. (Again  for  supersonic  axial  flow,  the  wake  velocity  does  not  Influence  the  pressure  distri- 
bution on  the  airfoil  surface.) 

Another  source  of  difficulty  Is  masked  In  delicate  subtlety  and  more  fundamentally  hard  to 
overcome  — the  breakdown  of  the  acoustic  theory  In  the  far  field.  As  Is  well-known,  acoustic  theory 
In  a moving  medium  Is  based  on  two  major  underlying  assumptions:  a disturbance  propagates  at  a 
uniform  acoustic  velocity  and  Is  swept  downstream  at  a constant  freestream  speed.  Although  this 
approximation  Is  sufficiently  accurate  In  the  vicinity  of  the  body,  the  acoustic  theory  for  a super- 
sonic flow  Is  manifestly  unfit  for  the  description  of  the  far  field;  It  falls,  for  example,  to  re- 
produce the  fanning  out  or  coalescence  of  Mach  waves.  The  reasons  for  the  breakdown  have  long  been 
understood:  as  a wavelet  spreads  out,  two  nonlinear  effects  Ignored  In  the  acoustic  theory  — the 
nonuniform  acoustic  and  flow  velocities  which  vary  with  both  place  and  time  — emerge  to  exert  their 
Influence  over  a large  distance.  The  nonlinear  effects  are  locally  small  everywhere.  Including  the 
far  field.  However,  not  only  Is  the  disturbance  at  a given  point  Influenced  by  the  localized, 
slightly  perturbed  flow  properties  but  It  has  been  undergoing  a continual  distortion  while  propaga- 
ting through  a nonuniform  flow  field.  It  Is  this  cumulative  distortion  or  "memory"  content  of  the 
signal,  which  encroaches  upon  the  result  of  the  acoustic  theory  and  eventually  alters  It  In  the  far 
field.  The  breakdown  of  the  acoustic  theory  raises  a serious  concern  for  a supersonic  cascade  If  one 
takes  the  approach  along  the  same  line  as  adopted  In  the  subsonic  cascade  aerodynamics  where  the  up- 
wash  generated  from  all  the  airfoils  are  sunned  up  within  the  framework  of  the  acoustic  theory. 

4.  PASSAGE  APPROACH 

The  complication  arising  from  the  breakdown  of  the  acoustic  theory  In  the  far  field  can  be  cir- 
cumvented by  resorting  to  a simple  stratagem  of  what  we  call  the  passage  approach.  In  the  passage 
approach,  we  focus  our  attention  to  a strip  of  the  flow  field,  the  L-shaped  region  bounded  by  y,,  xi, 
xg  and  y2  axes  of  Figure  (2).  Then,  the  Infinite  cascade  arrangement  Is  replaced  by  a physically 
equivalent  requirement  of  flow  periodicity.  In  the  case  of  subsonic  flow,  this  approach  would  not 
offer  any  additional  advantage.  However,  In  the  present  supersonic  situation  where  the  range  of  In- 
fluence of  a given  point  In  the  flow  Is  limited  within  its  Mach  cone,  the  attractive  feature  of  the 
passage  approach  Is  that  It  enables  one  to  formulate  the  problem  Involving  only  those  points  close  to 
the  reference  airfoil;  for  example,  the  Mach  cone  emanating  from  a point  along  the  yj  axis  of  Figure 
(2)  and  located  far  from  the  leading  edge  completely  misses  the  flow  passage  between  the  first  and 
second  airfoil  and  consequently  In  this  formulation  the  pressure  distribution  upon  the  airfoil  Is 
unaffected  by  It.  Since  the  formulation  Involves  only  the  near  field,  the  flow  field  of  Interest  can 
be  closely  approximated  by  the  acoustic  theory,  provided  the  thickness  of  the  airfoil  Is  sufficiently 
thin  and  can  be  neglected.  (We  will  later  see  that  for  supersonic  airfoils  of  conventional  design, 
although  their  contour  Is  relatively  thin.  It  is  not  quite  slender  enough  to  be  described  adequately 
by  the  acoustic  theory.  Nevertheless,  even  under  such  a circumstance,  the  passage  approach  still 
extends  Its  convenience  In  circumventing  the  trouble  of  sunrolng  up  all  the  upwash.)  Restricting  our 
attention  for  the  time  being  to  flat-plate  airfoils,  let  us  examine  how  we  can  apply  the  passage  ap- 
proach to  determine,  within  the  framework- of  the  acoustic  theory,  the  flow  In  the  preinterference  zone. 

If  we  focus  our  attention  on  the  flow  entering  the  passage,  It  Is  of  course  not  equal  to  the  up- 
stream flow  and  must  be  determined  as  a part  of  the  solution.  If  the  Initial  conditions  along  the 
yi  axis  of  Figure  (2)  are  assumed  to  be  known  as 

?'(x  - 0,  y)  « f(y)  ; (x  * 0,  y)  « g(x) 
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where  ♦'  Is  the  spatial  amplitude  of  the  perturbed  velocity  potential,  then  along. the  boundary  condi- 
tions, It  Is  quite  straightforward  to  determine  the  flow  downstream  of  It.  The  solution  of  the  acous 
tic  equation  subject  to  the  Initial  conditions  and  the  boundary  condition  for  the  first  airfoils  Is 
given  by  equation  (3.4)  of  [1].  The  expression  of  4'  given  therein  Is  applicable  to  the  preinter- 
ference zone  where  the  effect  of  the  second  alrfbll  1$  not  felt  (the  shaded  region  of  Figure  (2)) 
and  It  obviously  Involves  the  Initial  conditions,  f and  g,  where  are  unknown.  However,  the  unknown 
Initial  conditions  must  be  equal  to  the  velocity  potential,  4',  at  the  corresponding  point  on  the  y^ 
axis,  with  the  exception  of  the  phase  lag.  Consequently  the  initial  conditions  may  be  replaced  by 
the  velocity  potential  4*  and  we  obtain  the  following  equation  for  the  velocity  potential: 

x - my 

4'(x,  y)  ■ jj  J V(t)  exp[-  ^ (x  - t)]J0(^  /(x  - x)s  - mV ]dx 
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where  V(x)  Is  the  normal  component  of  the  fluid  velocity  on  the  reference  airfoil,  m » (M*  - l)1* 

where  M Is  the  freestream  Mach  number,  c Is  the  chord,  k Is  toe  >£/U„(l£  - 1)  where  u Is  the  fre-  | 

quency  of  oscillation  and  Is  the  freestream  velocity,  jQ  and  jj  are  given  by 

\ 

J0(x,  y.  x)  - Jo[fiT  T)  ■ 2 Jl,iTc  ; * * /x*  - (x  - yV  m * j 

where  Jq  and  Ji  are  Bessel  functions  of  zeroth  and  first  order,  respectively,  p Is  the  Interblade  j 

phase  angle  and  S|  and  t are  parameters  defining  cascade  geometry  as  given  In  Figure  (2).  (Strictly  i 

speaking,  the  above  equation  (2)  1$  valid  only  for  points  belonging  to  the  doubly  hatched  region  of  1 

Figure  (3)  of  (1)  but  this  restriction  Is  not  essential.)  We  note  that  In  (2),  the  first  term  on  the  J 

right  hand  corresponds  to  the  solution  for  an  Isolated  airfoil  and  known.  The  bracketed  expression  j 

In  the  second  term  Involves  4'  In  difference  form,  while  the  third  term  contains  4'  under  the  Inte-  j 

gral  sign.  Thus  we  call  the  equation  (2)  as  the  Integro-dlfference  equation.  We  also  observe  that  ; 

the  Integro-dlfference  equation  Is  'two-dimensional'  In  the  sense  that  there  are  two  Independent 
variables,  x and  y.  We  have  to  solve  (2)  for  the  unknown  4'  to  determine  the  preinterference  zone. 

6.  LOU  FREQUENCY  SOWTIOH 

In  our  earlier  Investigation,  (1],  it  was  found  that  by  confining  our  attention  to  the  low  fre- 
quency solution.  It  was  possible  to  obtain  a solution  of  the  Integro-dlfference  equation  In  closed- 
form.  According  to  the  scheaa,  we  expand  all  the  variables  In  terms  of  the  frequency  parameter  and 
retain  only  tanas  linear  to  frequency.  Thus,  for  Instance, 

4’  - 4(0  ♦ IB  4<‘)  ♦ ...  ; V • V(0  ♦ 16  V<‘)  , 

where  B » wc/U„.  Then  the  Integro-dlfference  equation  (2)  becomes  considerably  sliqpllfled;  the  j 

aquations  corresponding  to  4O  and  4 v‘)  are  given  as  equation  (3.17)  and  (3.18)  of  (1),  respectively.  ! 

Thf  solution  for  the  simplified  Integro-dlfference  equation  Is  given  by  equation  (3.26)  of  the  same  j 

paper,  where  the  first  term  corresponds  to  an  Isolated  airfoil  solution  and  the  other  term  repre-  I 

sents  the  correction  due  to  cascade  arrangement.  For  the  limit  of  sonic  leading  edge,  the  second 

— 59  — 


♦ j exp(-1k  |)elM 

♦ j exp(-1k  jr)e*y 

- mT*'(si*  t 

00 

y 

u*i(sr 

- Wc  ♦,{si*  T 


ten*  is  reduced  to  zero,  leaving  only  the  term  corresponding  to  an  Isolated  airfoil  solution.  This 
is  what  should  be  expected  because  of  the  lack  of  any  preinterference  in  the  sonic  leading  edge  limit. 

Pursuant  to  this  determination  of  the  preinterference  zone,  a solution  for  the  flow  downstream 
of  it  can  be  obtained,  including  those  regions  affected  by  the  wake.  The  complete  expression  for  the 
unsteady  pressure  distribution  is  given  in  closed  form  by  the  equations  (4.6)  and  (4.7)  of  Part  2, 

11).  (There  Is  a typographical  error  in  the  equation  (4.6)  and  the  last  term  in  the  bracket  should 

to  00 

read  l H(x.  - 2n  tm)  instead  of  [ H(x.  - 2n  tm).)  The  unsteady  pressure  given  there  is  in 
n « 1 n = 0 

general  form  and  applicable  to  any  cascade  geometry.  For  cascade  geometries  of  practical  Interest, 
however,  it  can  be  reduced  to  considerably  simpler  form.  For  Instance,  if  one  chooses  the  cascade 
geometry  of  Figure  (3),  which  is  characterized  by  the  requirement  that  the  bow  shock  emanating  from 
the  leading  edge  of  the  airfoil  misses  the  preceding  airfoil,  the  unsteady  pressure  distribution  is 
given  in  the  following  form:  if  we  express  the  pressure  distribution  as 

ACp  * 16  4l)  0»n«i1ng)  ; Acp  = e1u)t  fc[°)  + ^jr  (torsion) 


where  hg  and  ag  are  the  amplitude  of  bending  and  torsional  motion,  respectively,  then  for  0 < x < 
c - Sj  + tm. 
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where  xn  is  the  position  of  the  torsional  axis  measured  from  the  leading  edge.  The  point  x ■ c - 
si  + ^corresponds  to  a point  on  the  pressure  surface  where  the  shock  emanating  from  the  trailing 
edge  of  the  preceding  airfoils  Impinges  (point  P of  Figure  (3)).  For  the  other  cascade  configura- 
tions, once  the  shock  pattern  is  specified,  the  pressure  distribution  of  [1]  can  be  reduced  to 
similarly  simple  form.  Applying  the  pressure  distributions  to  stability  calculation  based  on  a two- 
dimensional  model,  the  flutter  Is  shown  to  occur  at  the  torsional  motion  and  the  flutter  boundary  is 
strongly  affected  by  cascade  parameters  (see  Figures  (4),  (5)  and  (6)  of  Part  2,  11]);  bending  motion 
is  always  predicted  to  be  stable  according  to  low  frequency  analysis. 


For  the  purpose  of  practical  applications,  however,  an  extension  of  the  analysis  was  called  for, 
because  the  frequency  parameters  of  the  aircraft  engine  fan  is  of  the  order  of  unity  and  cannot  be 
categorized  as  belonging  to  the  low  frequency  case. 


6.  EXTENSION  TO  HIGHER  FREQUENCY  RANGE 


Hence,  subsequent  to  this  low  frequency  analysis,  we  embarked  on  the  extension  of  the  analysis 
to  the  frequency  range  of  practical  Interest.  It  was  found  that  the  best  way  to  solve  the  integro- 
dlfference  equation,  (2),  is  to  write  4'  as 

x - my 

4'  » i f V(t)  exp(-  (x  - T)IJntii*-  Ax  - t)2  - m*y*]dT  + A + 8x  + Cx*  + Dy*  + Ex*  + Fxy*  ♦ ... 
m J c 0 M«C  (4) 

where  the  first  term  corresponds  to  an  Isolated  airfoil  solution  and  the  rest  of  the  terms  represents 
the  cascade  correction.  We  note  that  the  odd  terms  of  y,  like  y,  Ay  and  x*y  are  missing  in  the  cas- 
cade correction  term  and  this  is  due  to  the  asymnetry  of  the  integro-dlfference  equation  with  respect 
to  y.  The  coefficients.  A,  B,  ....  can  be  determined  by  substituting  them  into  the  integro-dlfference 
equation.  It  turns  out  that  due  to  the  very  nature  of  the  difference  form  involved,  one  cannot 
determine  coefficients  in  a successive  way;  they  have  to  be  determined  simultaneously  by  solving  an 
Infinite  number  of  linear  algebraic  equations.  For  all  practical  purposes  It  Is  sufficient,  however, 
to  retain  the  coefficients  up  to  F and  they  have  been  solved  in  explicit  form;  numerical  checks  show 
that  such  a solution  can  satisfy  the  Integro-dlfference  equation  quite  closely  for  the  frequency  para- 
meters of  practical  Interest.  As  an  example  of  such  approximate  solution.  A,  of  equation  (4)  can  be 


written  out  as 
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This  new  solution  reveals  nonuni  fortuity  associated  with  the  previous  scheme  of  frequency  ex- 
pansion discussed  In  Section  5.  If  we  once  again  look  at  the  equation  (3),  we  note  that  the  low 
frequency  expression  possesses  a factor  1 - i'u  as  Its  denominator  and  hence  becomes  singular  at 
zero  Interblade  phase  angle  (w  ■ 0).  As  pointed  out  In  [1],  from  this  singular  behavior  the  unique- 
incidence  effect  can  be  recovered  In  the  steady  limit.  However,  the  existence  of  the  singularity  at 
aqy  frequency  Implies  that  even  In  the  unsteady  case,  the  pressure  distribution  would  become  Infini- 
tely large  at  this  zero  Interblade  phase  angle.  The  reason  for  this  anomaly  of  the  previous,  low- 
frequency  solution  can  readily  be  observed  If  we  Inspect  the  expression  of  the  denominator  of  (5). 
Namely,  as  long  as  the  Interblade  phase  angle,  u»  Is  sufficiently  far  from  the  value  of  zero,  one  can 
legitimately  expand  the  expression  In  the  power  series  of  the  frequency  parameter,  k.  In  such  an 
expression,  the  denominator  becomes  the  leading  term,  2(1  - e1(J)*»  and  one  can  show  that  the  expres- 
sion In  fact  agrees  with  the  corresponding  one  obtained  by  the  low  frequency  analysis.  However,  the 
leading  term  (and  also  the  next  term)  of  the  denominator  of  (5)  vanishes  at  zero  Interblade  phase 
angle  and  consequently  the  formal  expansion  In  the  frequency  parameter  breaks  down  at  this  point;  In 
other  words,  the  regular  perturbation  scheme  In  the  frequency  1$  not  uniformly  valid  (though  In 
practical  flutter  calculation,  this  particular  situation  of  zero  Interblade  phase  angle  Is  of  little 
significance).  The  present  Improved  solution  such  as  (S)  removes  this  nonuniformity.  In  addition  to 
extending  Its  range  of  applicability  to  the  practical  frequency.  The  expression  of  (5)  becomes 
singular  only  If  both  p and  k simultaneously  become  zero  and  from  this,  one  can  again  obtain  the 
unique-incidence  effect  as  the  steady  limit  In  exactly  the  same  manner  as  described  In  Section  4 of 
(lj.  Except  for  this  condition,  there  Is  no  other  singular  point  In  (5)  (see  [12]). 


The  wake  velocity  can  be  determined  likewise  by  expanding  It  In  a Taylor  series  of  x about  the 
trailing  edge  and  we  obtained  a closed-form  expression  for  the  pressure  distribution  now  valid  up  to 
and  Including  the  frequency  parameter  of  practical  Interest. 


Figure  (4)  shows  an  example  of  the  unsteady  pressure  distribution  for  a cascade  arrangement 
corresponding  to  Figure  (3);  the  airfoils  (N  - 1.3,  solidity  » 1.0,  stagger  angle  - 60°)  are  exe- 
cuting torsional  motion  at  Its  mid-chord  with  Interblade  phase  angle  of  90°  (a)  and  180°  (b).  The 
solid  lines  correspond  to  the  solution  obtained  by  the  present  scheme  and  the  dashed  lines  to  the 
low-frequency  solution,  equation  (3).  It  Is  of  Interest  to  observe  that,  at  these  Interblade  phase 
angles,  the  low-frequency  solution  holds  up  remarkably  well  In  comparison  with  the  more  exact  present 
solution.  The  stability  analysis  shows  that  the  region  of  torsional  Instability  tends  to  shrink  as 
the  value  of  the  frequency  parameter  Increases,  but  even  at  the  frequency  parameter  of  unity  where 
the  Isolated  airfoil  becomes  completely  stable,  there  still  remains  a considerable  region  of  In- 
stability for  cascaded  airfoils.  As  a matter  of  considerable  Interest,  the  flutter  analysis  for  the 
bending  motion  shows  that  although  at  low  frequency  the  bending  motion  Is  always  predicted  to  be 
stable,  at  the  frequency  parameter  of  practical  Interest,  'Islands'  of  bending  instability  emerge. 

Me  feel  that  the  results  are  encouraging,  because,  as  mentioned  In  Section  2,  the  bending  mode  In- 
stability as  well  as  the  torsional  one  was  observed  for  the  actual  fans  In  supersonic  unstalled 
flutter. 


7.  CUMUUTIVB.  NONLINEAR  EFFECT  OF  AIRFOIL  SHAPE 

Me  here  recall  a remark  made  earlier  In  Section  2 on  the  first-order  effect  of  airfoil  shape 
upon  supersonic  unstalled  flutter.  As  mentioned  In  Section  3,  the  acoustic  theory  does  not  take 
Into  account  the  nonuniform  fluid  and  acoustic  velocity  produced  by  the  very  presence  of  the  airfoil 
with  thickness.  Consequently  one  has  to  depart  from  the  conventional  acoustic  theory  In  order  to 
capture  the  quintessential  feature  of  the  effects. 


Prior  to  embarking  on  the  description  of  the  analysis,  we  first  consider  what  one  can  anticipate, 
from  the  physical  reasoning,  as  the  consequence  of  nonuniform  flow  field.  Let  us  first  assume  that 
only  a single  point  on  the  alrfbll  Is  sinusoidally  oscillating.  Mhen  one  plots  at  a given  point  In 
flow  the  time-trace  of  the  disturbance  emitted,  the  departure  of  the  nonuniform  acoustic  and  con- 
vective velocity  from  the  uniform  ones  (acoustic  theory)  would  be  graphically  revealed,  mainly,  as 
the  phase  difference  between  the  actual  signal  and  the  one  predicted  by  the  acoustic  theory.  The 


ag  would  Increase.  Suppose  now  that  the  alrfbll,  as  a whole.  Is  oscillating.  Then  the  above  phase 


— 61  — 


for  an  Individual  disturbance,  which  differs  from  one  signal  to  another,  and  that  alone  — to  say 
nothing  of  the  modification  In  the  amplitude  of  each  signal  — could  Introduce,  when  signals  are 
vectorial ly  added,  a significant  correction  to  both  the  amplitude  and  phase  of  the  unsteady  flow  In 
the  far  field.  Thus  the  nonuniform  flow  field  associated  with  the  airfoil  shape  would  Introduce  a 
change  In  the  very  substance  of  the  fluctuating  pressure  (In  addition  to  the  usual  alteration  to  be 
made  to  the  direction  of  the  characteristics).  The  modification  Induced  to  the  far  field  signal  has 
the  following  Implication,  which  appears  to  warrant  sufficient  emphasis:  contrary  to  the  situation 
In  the  near  field,  the  unsteady  signal  at  a large  distance  — even  to  the  first  order  of  small  per- 
turbation — can  by  no  means  be  separated  from  such  effect,  as  the  airfoil  shape,  camber  and  angle 
of  attack,  which  cause  the  properties  of  propagation  nonuniform.  If  this  Is  the  case,  no  doubt  this 
appears  to  explain  the  significance  of  the  airfoil  shape  In  the  flutter  boundary  of  cascaded  airfoils. 


Me  Intend  to  confirm  these  physical  expectations  and  for  that  purpose  we  first  obtained  the  un- 
steady flow  off  the  surface  of  an  Isolated  oscillating  airfoil  with  thickness,  [13].  The  upper  and 
lower  surfaces  of  the  airfoil  are  assumed  to  consist  of  parabolic  arcs.  (Since  the  bow  shock  Is 
assumed  to  be  attached,  the  supersonic  flow  above  the  upper  surface  of  the  airfoil  Is  Independent  of 
one  below  and  conseguently  the  analysis  can  treat  the  effect  of  camber  as  well  as  thickness.)  The 
governing  equation  Includes  the  second-order  terms  which  amount.  In  the  far  field,  to  first-order 
unsteady  term;  this  can  be  accomplished  by  the  use  of  the  strained  coordinate  technique  of  Llghthlll, 
Hhltham  and  Lin.  We  relegate  the  details  of  the  analysis  to  [13]  and  at  present  It  suffices  to  point 
out  that  In  the  final  result,  the  unsteady  velocity  potential  can  be  given,  to  the  first-order  of 
small  perturbation,  by  the  following  Integral  representation: 
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where  Is  the  steady  velocity  potential,  the  shape  of  the  airfoil  (In  the  mean  position  of 
oscillation)  Is  given  by 

e(>j  a x2  + Bx) 

and  where  N » (y  + l)Mf/2(m)2 , s and  p are  characteristic  coordinates  based  on  the  steady  velocity, 
and  M Is  the  confluent  hypergeometric  function.  Contrary  to  the  acoustic  solution  which  breaks  down 
In  the  far  field,  the  present  solution  Is  uniformly  valid  In  the  entire  flow.  When  reduced  to  various 
limits.  It  agrees  with  such  known  solution  as  the  Whltham's  rule  at  the  steady  limit,  oscillating 
flat  plate  and  the  wedge  solution  of  Carrier  and  Van  Dyke.  More  Importantly,  the  above  solution  for 
a parabolic  airfoil  reveals  many  physical  features  relevant  to  the  propagation  of  unsteady  distur- 
bance through  nonuniform  flow. 

The  numerical  results  are  presented  In  Figure  (5),  where  the  unsteady  pressure  distribution  for 
a parabolic  airfoil  (thickness  to  chord  ratio  ■ 2.5*),  computed  from  equation  (6),  Is  shown  In  com- 
parison with  the  result  for  a flat-plate  airfoil  at  two  different  frequencies  of  oscillation  — 
mc/U  - 0.1  In  Figure  (5a)  and  u>c/U  ■ 1 In  Figure  (5b).  There,  both  the  amplitude  and  phase  are 
plotted  as  functions  of  s,  or  the  distance  of  the  root  of  a straight  Mach  wave  and  at  three  different 
locations  of  y.  We  observed  that  - though  for  u>c/U  *0.1  the  effect  of  the  airfoil  shape  does  not 
become  prominent  at  these  locations  — for  wc/U,,,  * I,  except  for  the  close  vicinity  of  the  leading 
edge.  It  Indeed  alters  the  pressure  distribution  significantly,  as  anticipated  by  the  physical 
considerations. 

This  result  for  an  Isolated  airfoil  appears  to  explain  the  Importance  of  airfoil  shape  in  super- 
sonic unstalled  flutter  and  currently  efforts  are  under  way  to  include  the  above  results  Into  cascade 
analysis.  Though  evident  from  the  foregoing  analysis.  It  appears  worthwhile  to  emphasize  that.  In 
all  phases  of  the  present  study,  efforts  are  focused  on  deriving  the  solution  In  closed-fOrm  or  at 
least  retaining  Its  analytical  structure  In  order  to  facilitate  the  examination  of  the  roles  played 
by  the  various  parameters  and  also  to  minimize  the  computational  time  expended  for  repeated  flutter 
calculations. 
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On  the  Issue  of  Resonance  in  an 
Unsteady  Supersonic  Cascade 

M.  Kurosaka* 

General  Electric  Company.  Schenectady,  N.  Y. 

IN  a recent  paper,  Verdon  and  McCune1  presented  a 
linearized  analysis  of  an  unsteady  supersonic  cascade  with 
subsonic  axial  velocity.  It  is  an  interesting  advancement  on 
Verdon's  previous  paper.-’  To  achieve  its  objective  of  com- 
puting the  pressure  distribution,  the  analysis  of  Ref.  I starts 
out  to  add  the  contributions  from  unsteady  disturbances 
generated  at  all  the  oscillating  airfoils  below  the  reference  air- 
foil in  the  cascade.  The  partial  sums  of  the  series  were  found 
to  oscillate  about  apparent  limiting  values  in  general  but,  un- 
fortunately, convergence  was  not  proved.  Aside  from  the 
problem  of  convergence.  Ref.  I also  reports  that  under  certain 
circumstances  the  numerical  scheme  broke  down  rather  inex- 
plicably. Thus  it  states  in  the  concluding  remarks  that  further 
work  is  needed  to  resolve  these  questions. 

Several  years  ago,  the  author  of  the  present  Note  en- 
countered essentially  the  same  series  and  found  that  it 
diverges  at  a certain  number  of  discrete  points,  although  these 
results  were  never  published.  The  appearance  of  Ref.  I. 
therefore,  seems  a Titling  opportunity  to  point  out  the 
divergence  of  the  series,  to  offer  it  as  a possible  explanation 
for  the  aforementioned  breakdown  of  the  numerical  scheme 
and  to  discuss  its  physical  implications  in  regard  to  resonance 
and  other  salient  points. 

Consider,  for  example,  the  following  kernel  function  Kix), 
which  appears  in  the  integral  equation  (23)  of  Ref.  I : 


derivative  of  7„  and  obtain 


Or 

Six)  = £ e 

n - I 


M'U.y)  =«r(  (jr  + nx, )-' iy  + nyA  )>)'  (3) 

We  will  show  that  Six)  is  divergent  for  certain  combinations 
of  parameters.  For  this,  rewrite  tYix.y)  as 


= U*  <n) -2z(rt)/(r.y)cos  <t>ix,y)  +Z1  (jr,y)  J 1 ’ (4) 

where 

Zix.y)  =*(•»' -mV  >'•'  (5a) 

<M.v,y)  = — r|  tanh  _ / (py/jr) 

-tanh  ‘iyyAlxA)\  -t  (5b) 

This  transformation.  Eq.  (5)  may  be  conveniently  achieved  by 
introducing  a set  of  auxiliary  variables  defined  by  x = p cosh 
v,  p y = p sinh  i>,  x,  = r>A  cosh  •>,,  and  py,  =pA  sinh  v, . Now 
from  Neumann’s  addition  formula,'  /„  can  be  expressed  as 


7i>(K)  = E «n,7„l(Z)7„.(z)cos  md 


Kix)  = - (//p)S(.v) 


where  Six)  is  given  by 
- / 

S(jr)=-  E kit!ny4e'" 


x J' l*l<* - n*  i • ; — ( M*.v , ) -T’J 

| (sr-rur,  )•’  - (p/ty , j ’■  1 

where  jr,  and  yt  are  (he  spatial  distances  between  the  ad- 
jacent airfoils  in  the  cascade  defined  in  Ref.  I,  p = (W-  — I ) ’’ 
where  M is  the  Mach  number,  jr .,  -py,  2 0 (subsonic  axial 
velocity).  * is  a frequency  parameter,  and  tl  is  related  to  the 
interblade  phase  lag  o.  The  nth  term  in  the  sum  represents  the 
influence  of  disturbances  generated  at  the  nth  (n<0)  below 
the  first  airfoil.  We  replace  n by  - n,  rewrite  7,  in  terms  of  the 
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, = / for  m = 0 


<„,  =2  for  m = 1,2,... 


for  any  complex  values  of  Z.  z,  and  <t>.  Substituting  Eq.  (6)  in- 
to Eq.  (2)  and  assuming  the  validity  of  the  interchange  of  the 
order  of  differentiation  and  summation,  we  obtain 
00 

SU>  = E »w 


d[J„\Z{x,y)  jeos  md(jr,y)  | 
dr 


where  is  given  by 

Ob 

«„■  = E e ""'7„.  | nk  ( x-\  - pV, ) ' ’ I (8) 

» - 1 

In  deriving  Eq.  (7),  we  have  made  explicit  use  of  the  fact  that 
z(n)=n*(x^-pVi)''  which  appears  in  Eq.  (8).  is  in- 
dependent of  x and  y.  while  7.  and  « are  independent  of  n.  The 
series  of  a„  is  called  a Schlomlich  series.  Here  we  want  to 
single  out  a„  and  examine  its  real  part,  which  is  given  by 

OB 

Re|o„|=  E cos  (nn>70|»*(ari-pV4)'’  1 (9) 
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Now  from  Ref.  4,  for  any  r and  Jr  beiween  0 and  x,  ihe  series 

OO 

2 cos  (n/)7„(ru r) 

H = I 

is  divergenl  at  /= jr  and  therefore  Eq.  (9)  diverges  at 

It  is  easy  to  observe  that  Eq.  (9)  has  another  singularity  at 
0=  -AOrJ-^yJ)1’ 

or  more  in  general  has  an  infinite  number  of  singularities  at 

0=  ±k{x2A- y!yA)'‘  + 2xn  n = 0.  ±7,  ± 2 .... 

If  we  express  0 by  the  interblade  phase  angle  a,  this  becomes 

a + kMxA  = ±*(jr  2A-aiyA)u 

+ 2xn  n=0,  ±1,  ±2,...  (10) 

At  these  points  the  series  S(jr)  diverges.  We  take  special  note 
that  the  divergent  condition  for  Six)  is  independent  of  the 
value  of  rr. 

One  can  check  the  divergent  condition.  Eq.  (10).  more 
directly  by  the  numerical  evaluation  of  the  original  series  Eq. 
(I)  in  the  vicinity  of  divergence.  Let  the  departure  from  the 
divergence  be  A.  i.e.. 

o + kMxA  ±k^x2^-|l^y^4),’  -2xn  = A 

where  A = 0 corresponds  to  the  divergent  condition.  Eq.  (10). 
The  partial  sum  of  Eq.  (1)  can  then  be  written,  after  replacing 
n by  -A,  as 

v 

Ss(x)  = ky2  5^  cxpl  -/7i|  ±x(jr'j-M;Vj>'’  + -l|  l«>'i 

* -- 1 

J,[kUx  + iixA)2 -IgAy ., ) ' I' ' I 
| ( jr  + njr  j ) •’  - ( ) ; 1 1 ! 

A number  of  numerical  checks  have  been  performed  and  they 
in  fact  confirm  the  divergence  of  the  series  at  A = 0.  For  exam- 
ple, for  the  cascade  A of  Ref.  I , the  real  and  imaginary  part  of 
Sv(jr)  (divided  by  ky!)  at  x = 0.5  are  shown  in  Fig.  I as  func- 
tions of  A for  the  various  values  of  A ( where  the  positive  sign 
in  the  argument  of  the  exponential  function  of  Eq.  (II)  is 
chosen  | . As  can  be  observed  immediately,  the  series  tends  to 
diverge  at  <5  = 0.  In  addition,  we  note  the  following  important 
points; 

a)  The  divergence  does  not  appear  when  A 
is  5 or  20.  It  starts  to  emerge  at  A = 200  and  becomes 
prominent  at  A = 2000.  In  other  words,  Ihe  cause  of  the 
divergence  is  not  the  effect  of  the  nearby  airfoils  but  the 
cumulative  influence  of  those  airfoils  located  far  from  the 
reference  airfoil. 

b)  As  A increases,  the  effect  of  jr  on  the  r)th  term  in  the 
series,  Eq.  (II),  becomes  insignificant. t When  we  combine 
this  with  a),  it  becomes  evident  why  the  divergent  condition  of 
Eq.  (10)  forS(jr)  is  independent  ofx. 

c)  For  a given  value  of  A away  from  the  divergence  con- 
dition, the  partial  sum  5v(Jr)  appears  to  oscillate  as  A in- 
creases. This  behavior  is  in  agreement  with  what  was  reported 
in  Ref.  I. 

Other  series  in  Ref.  I can  also  be  shown  to  diverge.  For  in- 
stance, Eq.  (20)  of  Ref.  I contains  a series,  which,  in  the 

♦ This  enables  one  to  check  the  divergent  condition  directly  by  using 
the  asymptotic  formula  of  J t in  Eq.  ( 1 1 ). 


region  downstream  of  the  Mach  wave  emanating  from  the 
leading  edge  of  the  first  airfoil,  becomes 

T(x.ri,y)=  £ e ‘"aJ„[k\ {x-ii  + nxA  )> 

H = I 

-n!{y  + nyA)!)':  I (12) 

By  a method  exactly  identical  to  that  used  for  S(x),  we  can 
readily  show  that  T{x,r),y)  again  becomes  divergent  at  the 
condition,  Eq.  (10). 

For  the  cascade  A of  Ref.  I,  the  linear  relationship  beiween 
o and  k as  indicated  in  the  divergent  condition,  Eq.  (10)  is 
plotted  in  Fig.  2 for  several  values  of  n.  As  mentioned,  break- 
down of  the  numerical  scheme  was  reported  in  Ref.  1 for  cer- 
tain combinations  of  parameters;  in  Fig.  7 of  Ref.  I,  for  a 
given  value  of  k,  o was  continuously  increased  starling  from  0 
and  the  failure  occurred  at  a certain  value  of  a.  When  we  add 
such  combinations  of  a and  k closest  to  the  breakdown 
situation  to  Ihe  present  Fig.  2,  where  they  are  designated  by 
circled  points,  breakdown  is  observed  to  occur  near  the 
divergent  condition. 

The  present  divergent  criteria,  Eq.  (10),  are  identical  with 
those  obtained  by  Samoilovich.!  Fie  obtained  it  under  the 
restricted  condition  corresponding  to  x = T)=y  = 0 in  Eq.  (12). 
In  such  a case,  the  argument  of  J0  becomes  proportional  to  n 
and  the  use  of  the  Poisson  summation  formula  readily  enables 
the  series  to  be  transformed  into  another  series  represen- 
tation, which  possesses  the  singular  points  at  Eq.  (10).  In  the 
present  derivation,  the  divergent  condition  has  been  obtained 
for  the  more  general  case  of  any  finite  values  of  jr,  ij,  and  y. 

It  is  a matter  of  considerable  interest  to  observe  that  Eq. 
(10)  is  formally  the  same  as  the  resonance  condition  in  a sub- 
sonic cascade,*  and  therefore  it  may  be  given  the  physical 
meaning  similar  to  that  discussed  in  Ref.  7.  Thus,  for  the 
present  divergence  of  the  series  associated  with  the  supersonic 
cascade,  Samoilovich  also  gave  Ihe  physical  interpretation  of 
"resonance."  Contrary  to  the  subsonic  cascade,  it  is, 
however,  highly  unlikely  that  resonance  at  these  conditions 
could  indeed  occur  for  a supersonic  cascade.  The  reason  is  as 
follows.  We  have  found  that  the  divergence  is  the  direct  con- 
sequence of  Ihe  cumulative  contribution  of  those  airfoils 
located  far  from  the  reference  airfoil;  in  Fig.  I,  the  peak  does 
not  appear  when  Ihe  number  of  the  preceding  airfoils  A is  5 or 
20,  and  it  starts  to  emerge  only  at  A = 200.  Needless  to  say, 
such  an  effect  is  computed  within  the  framework  of  linearized 
formulation.  It  is  well  known,  however,  that  the  linearized 
treatment  of  a supersonic  flow  breaks  down  in  the  field  far 
from  an  airfoil.  According  to  the  linear  theory,  disturbances 
created  by  the  airfoil  would  propagate  unattenuated,  even  to 
infinity,  for  both  the  low-  and  high-frequency  limits  (for  high- 
frequency  behavior,  see  Ref.  8.)  Although  in  the  near  field  Ihe 
linearized  theory  is  a good  approximation,  the  contribution  of 
the  nonlinear  terms  become  no  longer  negligible  in  the  far 
field  and  there  it  encroaches  on  and  modifies  the  effect  predic- 
ted by  the  linear  theory.  Physically,  this  follows  from  the  fact 
that,  by  the  time  disturbance  reaches  the  far  field,  two 
nonlinear  effects  ignored  in  the  acoustic  theory— the  con- 
vection of  disturbance  by  the  local  and  instantaneous  fluid 
velocity  and  its  propagation  at  the  nonuniform  speed  of 
sound— have  cumulatively  taken  their  toll  and  distorted  the 
shape  of  the  wavelet  given  by  the  linearized  theory.  Thus,  as 
Lighlhill9  puts  it:  ".  ..  the  failure  of  linearized  theory  ...  is  ex- 
plained by  the  fact  that  ...  while  yielding  adequate  results  in  a 
limited  region,  may  yield  a worse  and  worse  approximation  to 
the  solution  farther  and  farther  from  where  the  boundary 
conditions  determining  the  solution  were  applied."  Con- 
sequently, only  the  influence  of  a limited  numbers  of  airfoils 
(A9  5,  say,  for  a typical  cascade)  in  the  neighborhood  of  the 
reference  airfoils  can  accurately  be  predicted  by  the  linearized 
theory.  As  the  distance  from  the  reference  airfoils  increases, 
the  nonlinear  effect  would  rapidly  alter  the  linearized  con- 
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tribution  whose  pro  forma  sum  over  a laige  number  of  air- 
foils is  responsible  for  “resonance.”  These  considerations 
therefore  cast  a serious  doubt  on  the  prospect  of  resonance  at 
the  present  divergent  condition  in  an  actual  supersonic 
cascade,  unless  evidence  to  the  contrary  is  found.  (When  the 
Mach  number  is  increased  in  the  present  subsonic  leading  edge 
problem,  the  solution,  including  the  behavior  at  ‘resonance,- 
should  approach  the  lowest  limit  of  the  supersonic  leading 
edge  problem  of  Ref.  10.  Since  the  latter  does  not  possess  any 
such  ‘resonance,’  this  appears  to  provide  additional  evidence 
to  the  contrary.) 

Even  when  the  sum  is  convergent,  the  aforementioned 
failure  of  the  supersonic  linearized  theory  in  the  far  field 
raises  the  obvious  question  on  the  very  scheme  of  the  sum- 
mation over  an  infinite  array  of  cascaded  airfoils  within  the 
framework  of  linearized  analysis.  We  hasten  to  add,  however, 
that  this  does  not  mean  that  one  cannot  solve  (he  supersonic 
cascade  problem  by  linearized  analysis.  It  is  indeed  feasible  to 
accomplish  this  by  using  a different  formulation  adopted  in 
Ref.  II.  The  formulation  may  be  called  a “passage  ap- 
proach” and  attention  is  focused  on  the  reference  passage  be- 
tween two  blades.  In  place  of  the  infinite  cascade 
arrangement,  we  impost  an  equivalent  periodicity  re- 
quirement that  the  flow  at  any  given  point  in  the  reference 
passage  be  the  same  as  the  flow  at  the  corresponding  point  in 
the  adjacent  passage  (with  the  exception  of  an  interblade 
phase  lag).  Since  the  problem  is  set  up  within  the  blade 
passage  and,  furthermore,  because  of  the  limited  domain  of 
dependence  of  a point  in  a supersonic  flow,  it  turns  out  that 
the  formulation  involves  only  the  near  field  where  the 
linearized  theory  is  valid.  This  approach  leads  to  a closed- 
form  expression  for  the  pressure  distribution,  correct  to  the 
linear  order  of  the  frequency.  The  analysis  yields  the  so-called 
unique  incidence  effect  as  the  steady  limit,  agrees  with  the 
sonic  limit  of  the  supersonic  leading  edge  analysts,  and 
correctly  reflects  the  effect  of  the  back  pressure.  The  ex- 
pression does  not  suffer  breakdown  ai  Eq.  (10)  and  this  is  one 


more  reason  to  preclude  resonance  under  these  conditions!  in 
a supersonic  cascade. 
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In  this  paper  we  examine  how  the  unsteady  flow  field  radiated  from  an  oscillating 
body  is  altered  from  the  result  of  acoustic  theory  as  the  direct  consequence  of  dis- 
turbances propagating  through  the  non-uniform  flow  produced  by  the  presence  of 
the  body.  Taking  the  specific  example  of  an  oscillating  airfoil  placed  in  supersonic 
flow  and  having  the  contour  of  a parabolic  arc,  we  derive  a closed-form  representation 
for  the  unsteady  flow  field  in  terms  of  the  confluent  hypergeometric  function.  The 
analytical  expression  reveals  explicitly  that,  though  the  body  shape  has  a negligible 
effect  in  the  near  field,  it  inextricably  affects  the  unsteady  flow  at  a large  distance, 
both  in  its  amplitude  and  phase,  and  substantially  modifies  the  results  of  acoustic 
theory.  In  addition,  we  display  the  relation  of  this  solution  to  the  ‘fundamental 
solution  ’ and  the  other  salient  physical  features  connected  with  disturbances  propa- 
gating through  non-uniform  flow.  The  present  results  recover  Whitham’s  rule  in  the 
limit  of  zero  frequency  of  oscillation  and  also  include,  as  another  special  case,  the 
unsteady  solution  for  a wedge  obtained  by  Carrier  and  Van  Dyke. 


1.  Introduction 

As  is  well  known,  acoustic  theory  in  a moving  medium  is  based  on  two  major 
assumptions*  that  a disturbance  propagates  at  a uniform  acoustic  velocity  and  is 
swept  downstream  at  a constant  free-stream  speed.  Although  this  approximation  is 
sufficiently  accurate  in  the  vicinity  of  the  body,  the  acoustic  theory  for  a supersonic 
flow  is  manifestly  unfit  for  the  description  of  the  far  field;  it  fails,  for  example,  to 
reproduce  the  fanning  out  or  coalescence  of  Mach  waves.  The  reasons  for  the  break- 
down have  long  been  understood  (e.g.  Ligh thill  1954):  as  a wavelet  spreads  out,  two 
nonlinear  effects  ignored  in  the  acoustic  theory,  i.e.  the  non-uniform  acoustic  and 
flow  velocities,  which  vary  with  both  position  and  time,  emerge  and  exert  an  influence 
over  a large  distance.  The  nonlinear  effects  are  locally  small  everywhere,  including 
the  far  field.  However,  not  only  is  the  disturbance  at  a given  point  influenced  by  the 
slightly  perturbed  flow  properties  at  that  location  but  it  has  been  undergoing  a con- 
tinual distortion  while  propagating  through  a non-uniform  flow  field.  It  is  this  cumu- 
lative distortion  or  ‘ memory  ’ content  of  the  signal  which  encroaches  upon  the  result 
of  acoustic  theory  and  eventually  alters  it  in  the  far  field. 

For  a steady  flow,  the  task  of  surmounting  the  shortcomings  of  acoustic  theory  has 
drawn  the  attention  of  Friedrichs  (1948),  Lighthill  (1949)  and  Whitham  (1950, 1952), 
to  mention  only  a few.  These  efforts  culminated  in  the  following  celebrated  rule  due 

t Present  address:  University  of  Tennessee  Space  Institute,  Tullahoma,  Tennessee  87388. 
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to  Whitham  f (1952):  to  a good  approximation,  the  result  of  acoustic  theory  can  be 
amended  if  one  replaces  the  linearized  Mach  wave  by  one  revised  using  linearized 
velocities  but,  along  this  improved  Mach  wave,  retains  the  values  of  the  fluid  proper- 
ties predicted  by  acoustic  theory.  Crudely  speaking,  then,  the  only  visible  consequence 
of  the  nonlinearity  is  the  directional  change  in  the  Mach  waves ; the  fluid  velocities 
remain  essentially  unchanged.  We  reiterate  here  that  the  flow  is  steady  in  the  frame 
of  reference  fixed  to  the  body. 

In  contrast  to  the  above  steady  flow  situation,  relatively  less  attention  appears  to 
have  been  paid  to  problems  where  the  flow  is  unsteady,  again  with  respect  to  the 
co-ordinate  system  fixed  to  the  body.  To  be  sure,  related  studies  have  been  published 
but  they  seem  mostly  to  be  restricted  to  a one-dimensional  problem  and  its  diverse 
variants  (e.g.  Lesser  1970;  Romanova  1970;  Nayfeh  1975).  There  have  been  very  few 
attempts,  if  any,  to  obtain,  in  the  spirit  of  the  above  steady  problems,  a complete  and 
uniformly  valid  solution  and  then  display  the  global  behaviour  of  the  unsteady  flow 
field  in  either  two-  or  three-dimensional  space.  Yet  there  are  many  important  practical 
problems,  like  the  unsteady  aerodynamic  interference  between  a multitude  of  oscil- 
lating bodies  in  a flow,  e.g.  flutter  of  cascaded  airfoils,  and  other  similar  phenomena, 
where  such  an  improved  prediction  of  the  unsteady  flow  valid  even  in  the  far  field  is 
critically  needed.  Prompted  by  this,  we  address  here  the  problem  of  obtaining  a first- 
order  correction  to  the  acoustic  field  radiated  from  an  oscillating  body,  accounting  for 
the  interaction  with  the  non-uniform  flow  created  by  the  body  itself. 

In  the  case  of  unsteady  flow,  the  nonlinearity  will  have  additional  consequences, 
as  one  can  anticipate  from  the  following  physical  reasoning.  Let  us  first  assume  that 
only  a single  point  on  the  body  is  oscillating  sinusoidally.  When  one  plots  at  a given 
point  in  the  flow  the  time  trace  of  the  disturbance  emitted,  the  departure  of  the  non- 
uniform  acoustic  and  convective  velocities  from  the  uniform  ones  (from  acoustic 
theory)  will  be  graphically  revealed,  mainly,  as  a phase  difference  between  the  actual 
signal  and  the  one  predicted  by  the  acoustic  theory.  The  phase  lag  depends  on  the 
position  and,  the  more  one  moves  away  from  the  source,  the  more  the  phase  lag  will 
increase.  Suppose  now  that  the  whole  body  is  oscillating.  Then  the  above  phase  lag 
for  an  individual  disturbance,  which  differs  from  one  signal  to  another,  and  that  alone 
(to  say  nothing  of  the  modification  in  the  amplitude  of  each  signal)  could  introduce, 
when  signals  are  vectorially  added,  a pronounced  correction  to  both  the  amplitude  and 
the  phase  of  the  unsteady  flow  in  the  far  field.  Thus  the  nonlinearity  would  cause,  in 
addition  to  the  alteration  to  be  made  to  the  direction  of  the  characteristics,  a change 
in  the  fluctuating  pressure  itself.  The  modification  induced  in  the  far-field  signal  has 
the  following  implication,  which  appears  to  warrant  emphasis:  contrary  to  the  situ- 
ation in  the  near  field,  the  unsteady  Bignal  at  a large  distance,  even  to  first  order  in  the 
small  perturbation,  can  by  no  meanB  be  separated  from  such  effects  as  the  body  shape, 
camber  and  angle  of  attack,  which  cause  the  properties  of  propagation  to  be  non- 
uniform.  The  effect  of  thickness,  for  example,  would  be  inextricably  embedded  in  the 
far-field  unsteady  signal. 

Our  present  aim  is  to  confirm  these  expectations  and  we  shall  do  so  by  investigating 

t In  early  literature  this  was  referred  to  as  Whitham's  hypothesis.  Now  that  it  has  become 
well  established,  it  appears  more  appropriate  to  call  it  a rule  instead.  This  rule  should  not,  of 
course,  be  confused  with  another  rule,  due  also  to  Whitham,  relevant  to  the  propagation  of  a 
shook  through  a region  of  varying  cross-sectional  area  (e.g.  Whitham  1974). 
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the  effects  of  a non-planar  body,  whose  presence  creates  a non-uniform  surrounding 
environment,  upon  the  unsteady  flow  field.  We  shall  expressly  limit  our  investigations 
to  the  case  of  a two-dimensional  slender  body  whose  upper  and  lower  surfaces  consist 
of  parabolic  convex  arcs  and  which  is  oscillating  sinusoidally  in  a supersonic  flow 
(figure  1).  Though,  by  confining  our  attention  to  this  particular  shape,  we  Bhall 
inevitably  forfeit  formal  generality,  the  present  approach  will  give  a closed-form 
solution  which  is  amenable  to  detailed  study ; from  this  we  hope  to  glean  the  essential 
features  of  the  flow  non-uniformities.  With  regard  to  figure  i again,  the  thickness  of 
the  body  is  characterized  by  a parameter  e and  the  amplitude  of  oscillation  by  60. 
We  shall  examine  the  cumulative  effects  of  the  second-order  terms,  which  ascend  in 
the  far  field  to  a first-order  unsteady  term  0(6 9).  There  are  three  second-order  terms, 
0(e&9),  0(6 5)  and  0(e*),  of  which  only  the  first  two  are  relevant  for  the  present  un- 
steady problem.  If  one  assumes  e > ffg,  one  can  discard  the  term  0(0J),  whose  presence 
would  cause  undesirable  higher  harmonics.  With  this  assumption,  we  are  now  in  a 
position  to  focus  attention  on  the  remaining,  0(6 9 e)  term,  which  represents  the  genuine 
coupling  effect  of  present  interest.  It  should  be  remembered,  however,  that,  as  pointed 
out  by  Hayes  (1954)  for  steady  flow,  only  a few  selected  second-order  terms  contribute 
cumulatively  to  the  first-order  effects.  Hence  we  shall  pick  out,  by  the  use  of  the 
strained  co-ordinate  technique,  those  terms  O(60e)  whose  cumulative  effects  amount 
to  O(60)  in  the  far  field.  Thus  our  aim  is  clearly  different  from  Van  Dyke’s  (1963q) 
second-order  theory  for  an  oscillating  airfoil  including  the  effect  of  thickness.  There, 
because  of  his  interest  in  the  flow  on  the  airfoil  surface,  combined  with  a situation 
involving  only  slow  oscillations,  he  used  a regular  perturbation  scheme  in  60  and  e 
and  obtained  a solution  to  O(60e);  consequently,  the  hierarchical  ascent  of  terms 
of  second  order  to  first  order  in  the  far  field  was  neither  expected  to  take  place  nor  was 
his  concern.  On  the  contrary,  our  interest  centres  on  just  such  an  evolutionary, 
ascending  process. 

In  the  next  section  we  shall  begin  with  the  governing  equation  and  simplify  it  in 
$3  by  employing  the  strained  co-ordinate  technique.  In  $$4  and  5,  we  shall  describe 
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the  procedure  for  solving  this  simplified  equation.  We  set  out  to  obtain  the  corre- 
sponding Riemann  function  appropriate  for  a parabolic  airfoil ; the  Riemann  function 
can  be  constructed  explicitly  and  exactly  in  terms  of  the  confluent  hypergeometric 
function.  With  the  Riemann  function  thus  obtained,  the  solution,  equation  (5.4), 
follows  from  it  without  much  difficulty.  In  §6,  before  embarking  on  a physical 
interpretation  of  the  solution,  we  pause  and  confirm  that  the  present  results  can  be 
reduced,  through  the  limiting  properties  of  the  confluent  hypergeometric  function, 
to  some  known  results.  In  the  limit  of  zero  frequency  of  oscillation,  we  shall  recover 
Whitham’s  rule;  for  an  oscillating  wedge  with  small  apex  angle,  the  present  result 
will  embrace,  as  a special  case,  the  solution  obtained  by  Carrier  (1949)  and  Van  Dyke 
(1953b).  We  shall  resume  the  discussion  of  the  curved  airfoil  in  §7,  where  we  observe 
that  Tricomi’s  (1949)  expansion  formula  for  the  confluent  hypergeometric  function 
is  ideally  suited  to  the  extraction  of  a physical  interpretation ; the  gradual  ascent  of 
second-order  terms  to  alter  the  acoustic  signal  in  the  far  field  will  become  effortlessly 
visible ; and  there  the  effect  of  body  shape  will  be  found  to  be  tenaciously  inseparable 
from  the  unsteady  flow  field.  This  will  be  followed  in  §8  by  further  description  of 
salient  physical  features  related  to  the  disturbances  propagating  through  the  non- 
uniform  medium. 

2.  Problem  formulation 

The  governing  equation  for  the  perturbed  velocity  potential  is  given  to  second 
order,  to  which  order  the  flow  can  still  be  regarded  to  be  irrotational,  by 

uao  uco 

= ^°{(y-D  (®x +<£  ot)  <**, + •„) + 2<J>,  + 2 

+^-(<I>,  <&*  + <!>,  <>„)},  (2.1) 

where  the  perturbed  velocity  components  («',  v')  are  related  to  0 by 

<DX  - u',  - v', 

Um  is  the  free-stream  velocity,  a,  the  speed  of  sound  in  the  free  stream,  Mm  — f/®/o«, 
m - (Jf*,- 1)1  and  y is  the  adiabatic  exponent  of  the  gas.  We  express,  according  to 
Van  Dyke  (1953a),  the  co-ordinate  of  the  moving  upper  surface  os 

y - ef(x)  - 99eMg{x),  (2.2) 

where  ef(x)  (e  < *'  ^signates  the  shape  of  the  body  in  its  mean  position  of  oscillation 
and  the  second  t.  represents  its  harmonic  motion  with  frequency  ut.  The  two  small 
non-dimensional  parameters  e and  9t  characterize  the  slenderness  of  the  body  and  the 
amplitude  of  motion,  respectively.  As  long  as  the  shock  remains  attached,  we  need 
consider  only  the  flow  above  the  upper  surface.  The  boundary  condition  on  the 
surface  of  the  airfoil,  as  given  by  Van  Dyke  (1953a)  to  second  order,  is 

%-(Um  + QtHcf'-O9ei**g')-iw0oe<«g-{ef-0oeMg)<lf„  at  y - 0.  (2.3) 

Also,  <t>  vanishes  upstream  of  the  bow  shock,  whose  position  moves  in  time.  Sinoe  the 
flow  variables  are  discontinuous  at  the  shock  and,  strictly  speaking,  do  not  possess 
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Figure  2.  Smoothing  technique  of  Van  Dyke. 

derivatives  there,  the  governing  equation  is  not  formally  satisfied.  Hence  in  principle 
jump  conditions  across  the  shock,  which  is  moving  and  whose  temporal  position  is 
unknown  a priori,  must  be  imposed  to  ensure  the  conservation  of  mass,  momentum 
and  energy  there ; this  would  introduce  complications.  However,  this  knotty  problem 
can  be  completely  circumvented  by  the  smoothing  technique,  which  was  first  devised 
by  Courant  & Friedrichs  (1948,  p.  365)  for  steady  flow  and  later  extended  to  the 
unsteady  case  by  Van  Dyke  (1953a).  We  first  imagine  that  an  extension  has  been 
added  to  the  leading  edge  of  the  actual  airfoil:  a sufficiently  smooth  and  flexible  tip 
of  such  a shape  and  moving  in  such  a way  as  to  prevent  the  formation  of  the  shock 
in  the  flow  above  the  upper  surface  (figure  2).  We  then  regard  the  desired  solution  as 
the  limit  as  the  extension  shrinks.  Once  this  device  has  been  employed,  as  here,  the 
need  to  impose  jump  conditions  at  the  shock  can  be  eliminated  for  the  solutions  up  to 
second  order.  (Also,  whenever  necessary,  we  shall  hereafter  regard  the  discontinuity 
in  the  flow  variables  at  the  shock  in  the  sense  of  the  above  limiting  process.)  The 
smoothing  technique  provides,  in  effect,  a formal  justification  for  the  following  point 
of  view:  the  global  behaviour  of  the  unsteady  flow  downstream  of  the  bow  shock  can, 
to  a good  approximation,  be  determined  essentially  independently  of  the  presence  of 
the  shock  and  various  complications  arising  from  its  motion  (except  in  the  close 
vicinity  of  the  shock,  where  such  a solution  fails) ; the  situation  is  akin  to  the  familiar 
steady  problem  (Whit ham  1952). 

Following  Van  Dyke  a little  further,  we  separate  the  perturbed  velocity  potential 
into  a steady  and  time-dependent  part  by  writing 

<t>  = e<f>(x,  y:  e)  + 0oexp  [i(eof  - kx)]  i]r(x,  y : e,  0O),  (2.4) 

where  k = The  first  term  represents  the  steady  base  flow  and  ijf  in  the 

second  term  corresponds  to  the  unsteady  flow ; our  interest  is  focused  on  ijr.  We  sub- 
stitute (2.4)  into  (2.1)  and  (2.3)  and  assume  e p 0O,  as  stated  in  the  introduction.  We 
thus  obtain  the  following  two  sets  of  equations:  for  <j> 

e [ - mVx*  + 4„]  = e*(MyU„)[m*(N -!)?)*  + #]x,  (2.5a) 

with  the  boundary  condition 

e0„  = e£k/'+«W'-/0w)  at  y = °.  (2.56) 

and  for  \jr 

0,[  - m*  \jrXI  + ir„-  (km/M)*  f]  = 2(edtJUv){M\>[m\N  - 1)(4X^X  + ^„^]X 

-ik[(2N- i)mt4>x$x  + NmtM„  + 4>x$w]-N{kmlMa,)t4>,f},  (2.0o) 

where  N «=  j(y  + l)Af*0/m,l  with  the  boundary  condition 
0O^„  - -09[Umg'  + (ikUmm»IM*,)g]e^ 

+€0of/'(i>1  - ik\jf)  + ( - 4>xg’ +g4>„Wkx-N,n\  at  y«0.  (2.66) 
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In  obtaining  the  above  equation,  some  simplification  on  the  right-hand  side  has  been 
achieved  by  using  the  expressions  for  the  left-hand  side  and  neglecting  terms  of  higher 
than  the  second  order.  It  should  be  noticed  that,  although  (2.5a)  is  nonlinear  in  <j>, 
equation  (2.6a),  the  basis  of  this  paper,  is  a linear^  function  of  >jr  involving  variable 
coefficients.  Both  <f>  and  \jr  vanish  upstream  of  the  bow  shock. 


3.  Application  of  strained  co-ordinate  technique 

The  right-hand  sides  of  (2.5a)  and  (2.6a)  are  of  higher  order  than  the  left-hand  sides. 
Consequently,  if  one  uses  a regular  perturbation  scheme,  they  successively  yield  the 
first-  and  second-order  equations  of  Van  Dyke  with  the  right-hand  sides  either  zero  or 
expressible  in  terms  of  the  first-order  velocities,  respectively  ; the  first-order  equation 
for  ifr,  in  particular,  is  the  (reduced)  acoustic  equation  and  it  obviously  precludes  the 
ascent  of  the  terms  on  the  right-hand  side  to  first  order.  In  order  to  achieve  our  stated 
objective  of  examining  such  an  evolutionary  process,  we  shall  employ  the  strained 
co-ordinate  technique  instead:  this  is  the  point  of  departure  of  the  present  analysis. 
Although  the  original  strained  co-ordinate  technique  developed  by  Lighthill  and 
Whitham  involves  only  a single  family  of  characteristics,  the  present  unsteady  prob- 
lem requires  two  families  of  characteristics  for  adequate  description  of  the  flow  field. 
It  is  therefore  convenient  to  use  Lin’s  (1954)  extension  of  the  strained  co-ordinate 
technique  (see  also  Oswatitsch  1962)  or  the  analytic  method  of  characteristics,  which 
enables  one  to  treat  the  case  of  two  families  of  characteristics.  According  to  this 
method,  the  independent  variables  (x,  y)  as  well  as  the  dependent  variables  are  to  be 
expanded,  with  the  characteristic  parameters  » and  p regarded  as  new  independent 
variables: 


x = afm(a,p)  + &fl){s,p)+  ..., 

(3.1a) 

y = y(0)(*,p)+ey(1>(8,p)+  , 

(3.16) 

etj>  = e<f>a)(8,p)+e*<f>a)($,p)+  ... , 

(3.2  a) 

Orf  = 0o (*.  P)  + e0o {s,p)  + 0J ^ (a, p)  + ... , 

(3.26) 

where  s and  p are  constant  along  the  corresponding  characteristic  curves,  respectively. 
With  respect  to  the  characteristic  curve,  we  first  observe  that,  comparing  (2.5a) 
and  (2.6a),  all  the  coefficients  of  the  second  derivatives  in  the  steady  equation  are 
the  same  as  the  corresponding  ones  in  the  unsteady  part.  This  dictates,  then,  that 
the  characteristic  curves  for  both  the  steady  and  the  unsteady  equations  are  identi- 
cally the  same  and  given  by 


dy 

dx 


•-conit 


dy 

d^-  P -conit 


if.  in ,«  ., . w i J i 


(3.3a) 


(3.36) 


f Besides  the  usual  shock  emanating  from  the  leading  edge  (and  the  one  at  the  trailing  edge, 
which  does  not  matter  for  the  flow  field  upstream  of  it),  no  additional  shook  is  created  owing  to 
the  motion  of  the  airfoil;  consequently,  the  entire  unsteady  flow  can  be  uniformly  described 
by  the  linearised  equation. 
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Into  these  we  substitute  (3.1a, b)  end  (3.2a)  and  equate  the  coefficients  of  equal  powers 
of  6.  We  then  determine  successively,  using  the  boundary  condition  (2.6b),  the  terms 
in  the  series  expansion;  while  the  zeroth-order  terms  in  (3.1),  z*0*  and  y*®,  give  the 
expression  for  *‘ie  characteristic  parameters  corresponding  to  acoustic  theory,  the 
first-order  tei  *<M  and  give  the  desired  nonlinear  correction.  We  direct  atten- 
tion towards  the  fact  that  the  procedure  is  dependent  wholly  on  the  steady  flow 
and  excludes  the  unsteady  part  (3.2b).  This  process  of  co-ordinate  stretching  in 
steady  flow  being  familiar,  it  suffices  here  to  write  down  the  following  results; 


4F  - - H(s)(VJm)f(s ), 
s *=  x-my-  emyN(M\,/Uv)d^m(s)lds, 
p - x + my-  (e/2Ux)(N  - 2)  - <t>a)(p)l 


(3.6  a) 
(3.6  b) 


where  H(s)  is  a unit  step  function.  The  above  expressions  for  s and  p have  been  put 
in  the  present  form  by  rewriting  the  results  corresponding  to  (3.1  a,  b).  Geometrically, 
s represents,  as  shown  in  figure  1,  the  root  of  the  straight  Mach  wave  passing  through 
a given  point  (x,  y)  and  along  this  s remains  constant  (Van  Dyke  1975) ; likewise,  p re- 
presents the  root  of  the  cross  Mach  wave,  along  which  p remains  constant.  (As  a matter 
of  fact,  the  constants  of  integration  in  (3.3),  their  choice  being  at  our  disposal,  are  so 
adjusted  that,  at  y - 0,x  = « = p.)  Equation  (3.4)  indicates  that  the  steady,  first-order 
velocity  potential  is  dependent  on  s only  and  it  obviously  embodies  Whitham’s  rule. 

Having  thus  specified  a and  p,  we  then  substitute  the  expansion  for  the  unsteady 
part  (3.2b)  into  (2.6).  In  obtaining  the  equation  for  the  leading  term  i^u,  we  proceed 
with  caution  and  retain  the  terms  associated  with  k on  the  right-hand  side  because, 
for  sufficiently  high  frequencies,  they  could  become  comparable  with  the  terms  on  the 
left-hand  side;  the  terms  not  associated  with  k can  be  neglected.  One  thus  obtains 

0#{-4l Wi+Jf  Nek^^  + ~ (N-l)ekpMW 

where  <f>'a)  designates  the  derivative  of  with  respect  to  s;  in  differentiating  we 
recall  and  envisage  the  smoothing  process  described  in  § 2 and  discard  the  term  asso- 
ciated with  the  delta  function.  (When  obtaining  (3.6a),  the  term  ( k/M, *,)*  in  the  braces 
initially  appears  as  {kjMw)*[\  - 2{NjUa,)e<f>'w]  but  the  second  term  in  the  square 
brackets  is  neglected.)  The  boundary  condition  (2.6b)  becomes 


d0{iekf'(8)  \j^l)  - mtjr^  + m\Jr^>  + V (a)  eik‘)  = 0 at  8 « p, 

where 

V(x)  = Uwg'(x)  + (ikUmm*/Ml)g(x). 

Also,  the  upstream  condition  becomes 

*0  for  s < 0. 

It  is  convenient  at  this  point  to  introduce  the  function  F defined  by 
y^>  - exp[t(efe/t/J  N(p  - s)f  <»(«)]*’. 


(3.6b) 


(3.6c) 


l 
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Then  (3.6a)  becomes,  to  the  order  consistent  with  the  present  approximation, 

<?.  ( v + ET I [■ - + (p  ■ - «W'(1)(‘)]  FP 

+[(£),+a£  <3Sa) 

with  the  boundary  conditions 

^{«fc[-i/'W-^^'»(e)]j,+Ji-J,-^K(#)^j-0  at  s = p,  (3.8 ft) 

and 

F = 0 for  a < 0.  (3.8c) 

In  (3.8a,  6),  the  factor  d0  is  retained  as  a reminder  that  the  equations  are  valid  to  order 
60,  the  higher-order  terms  such  as  those  0(eda)  in  (3.26)  being  neglected.  Our  aim  is  to 
obtain  the  explicit  solution  for  F and  we  shall  do  so  for  an  airfoil  whose  shape  consists 
of  parabolic  arcs. 


1 


I 


4.  Construction  of  the  Riemann  function 

If  0'«>(«)  were  either  zero  or  a constant,  (3.8a)  would  be  reduced  to  the  telegraph 
equation.  In  the  present  case  of  a parabolic-arc  airfoil,  f(x)  in  (2.2)  is  quadratic  in  x 
and  from  (3.4)  the  derivative  is  linear  in  s.  Thus  (3.8  a)  is  a second-order  linear 
hyperbolic  equation  whose  coefficients  are  variable  (and  linear  in  s).  It  is  well  known 
that  the  solution  of  any  second-order  linear  hyperbolic  equation  can  be  expressed  in 
the  form  of  an  integral  representation,  once  the  corresponding  Riemann  function  has 
been  obtained  (e.g.  Courant  & Hilbert  1962,  p.  449).  If,  in  general,  u satisfies 

lF[u]  = + aux  + buy  + cu  = 0, 

where  a,  b and  c are  given  functions  of  x and  y,  then  it  can  be  represented  by  an  integral 
along  the  boundary  (where  Cauchy  data  are  assumed  to  be  prescribed)  whose  integrand 
involves  the  Riemann  function  R of  the  operator  SF.  R does  not  satisfy  the  operator 
equation  SF(R)  = 0 but  rather  satisfies  the  adjoint  operator  equation 

X*[R]  s Rxy- ( aR)x - (bR)v +cR  = 0. 

For  our  purpose,  it  is  convenient  to  derive  first,  instead  of  R,  the  Riemann  function 
R*  of  the  adjoint  operator  which  satisfies  the  operator  equation  for  SF itself;  then  we 
obtain  R through  the  symmetry  property  of  the  Riemann  functions.  For  the  present 
equation  (3.8a),  the  Riemann  function  of  the  adjoint  operator  R*(£,y,8,p)  satisfies 
the  following  three  conditions  (Courant  & Hilbert,  ibid.): 


(a)  JF^R*]  - Rt„+jj-jk[-(2N- l)^(€)  + (V-i)If^»(i)]R* 

+ 

(b)  along  ^4C  in  figure  3, 


(4.16) 


l dR*  A 
Ri^m0  °n 


] 

J 


v'ji.  - J- 
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Figure  3.  Integration  contour. 

and  along  BC, 

on  r,  = p;  (4.1c) 

(c)  B*(s,p,s,p)  - 1.  (4.1  rf) 

Intagrating  (4.1e)  and  determining  the  constant  of  integration  from  (4. 1 d),  we  obtain 

X0(€,p;s,p)  = exp/i*.  (4.2) 

where 

M*  - ( i/Um)ek{(N - 1)0U>(£) - #»(»)}  - N[p(4>'«'(l)  - *'<«(«)) 

If  we  write 

- exp(n*)M(x),  (4.8) 

where 

* - - ( ilUm)NeHi  - sHn-p)*™.  (4.4) 

then  for  a parabolic-arc  airfoil,  for  which  la  a constant,  (4.1a)  is  reduced  to  the 
following  ordinary  differential  equation: 

uMm+(l-s)M'-aM  » 0, 

a . } + Umk{H«N . 


where 


(4.5) 
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This  is  known  as  Kummer’s  equation  and  its  only  solution  which  satisfies  (4. Id)  is 
the  following  confluent  hypergeometric  function  (e.g.  Slater  1960,  p.  2): 


defined  by 
where 

«»  = 

and 

Hence  (4.3)  becomes 


M = M(a,  1 ,2), 


Jf(«16,«)-l+f+g|f+...+?^  + .... 


0.2" 


bHn\ 


o(a  + i)(o  + 2)...(o  + n-l),  for  n = 1,  2, 
a0  = 1. 

R*(i,V>s<P)  =exp(ji*)Jf(a,  1,2). 


(4.6) 


(4.7) 


Along  £ - 8,  R*(8,i)  \8,p)  = 1 and  this  obviously  satisfies  the  remaining  requirement 
(4.1b)  for  the  Riemann  function. 

The  Riemann  function  R(£,ij;a,p)  may  be  immediately  derived  from  R*  through 
the  symmetry  property  of  the  Riemann  function  (Courant  & Hilbert  1962,  p.  464)  by 
replacing  £ and  tj  with  a and  p,  respectively.  Thus  we  obtain 

R(£,  ir,8,p)=  ex  p {/*)  M(a,  1 ,2),  (4.8) 

whore 

p = (ilUJek{(N-  1 ) [<j>™  (8)  - (£)]  - N[V(<t>™  (8)  - t'V  (£)) 

-(8<t>'"H8)-E<f,'vm}, 

a = l + Uxk(4kNMl<j>'V)-'1 

and 

2 = -(ilUJekNis-nip-v)#™. 


5.  An  Integral  representation  of  the  solution 

Once  the  Riemann  function  has  been  thus  derived,  one  is  in  a position  to  employ 
Riemann’s  formula  (Courant  & Hilbert,  ibid.)  to  obtain  the  integral  representation  of 
F in  (3.8a),  provided  that  Cauchy  data  are  prescribed  on  the  boundary.  Unfortunately, 
the  present  boundary  condition  (3.8b),  which  applies  along  the  segment  OA  of  figure  3 
(this  corresponds  to  the  x axis  of  figure  1),  is  not  Cauchy  data.  Rather,  it  expresses  a 
linear  relationship  between  the  function  F and  its  derivatives ; this  induces  some  com- 
plication. If  one  applies  Riemann’s  formula  to  the  contour  around  the  shaded  region  of 
figure  3 (OACBO'),  although  the  contributions  from  the  line  segments  AC,  CB,  BO'  and 
00'  vanish  identically,  one  ends  up  with  an  integral  along  OA ; since  it  turns  out  that 
the  integral  involves  the  value  of  F,  which  is  unknown  as  yet,  one  has  to  solve  a 
complicated  integral  equation  to  determine  it. 

The  difficulty  is  by  no  means  unique  to  the  parabolic  airfoil,  and  in  fact  the  same 
complication  arises  even  in  the  more  simplified  situation  of  a flat-plate  airfoil,  where 
<j>'m  is  zero.  In  such  a case,  (3.8a)  is  reduced  to  the  telegraph  equation  and  the  corre- 
sponding Riemann  function  is  a Bessel  function  (e.g.  Courant  ft  Hilbert,  ibid.): 

R(£,v;».P)  - ^{(i/if.ias-fHp-?)]*}.  (8.1) 

To  construct  the  flat-plate  solution,  Temple  ft  Jahn  (1945)  used  this  and  applied 
Riemann’s  formula  for  a closed  curve ; the  contour  around  the  shaded  region  of  figure  3 

* 
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is  none  other  than  their  path  of  integration.  Their  final  result  for  F at  a general  point 
in  the  flow  was  left  in  a somewhat  awkward  form  involving,  inside  the  integral  along 
the  segment  OA,  the  unknown  values  of  F to  be  evaluated  there,  although  the  un- 
desirable term  vanishes  from  the  integral  for  a point  on  the  surface  of  airfoil,  i.e.  on 
OA . It  turns  out,  however,  that  one  can  advance  a step  further  and  eliminate  the  term 
entirely.  More  specifically,  by  substituting  the  expression  derived  for  F on  the  surface 
of  the  airfoil  into  the  integral  representation  for  an  arbitrary  point  and  noting  an 
identity  involving  a product  of  Bessel  functions,  F can  be  written  exactly  as  the 
following  integral  of  the  Riemann  function: 


F(s,p)  = ~^jQ  V(r)eikrR{i  = r,i)  ® T,s,p)dr 


(6.2a) 


= Jo V{T)  e<kT  J°  (i [(tf  ~T)(p~ T)]*) dT • (8-26) 

This  expression  is,  of  course,  the  well-known  flat-plate  solution  obtainable  by  a 
number  of  other  methods  (e.g.  Miles  1959,  p.  50). 

Motivated  by  (5.2a),  in  the  present  case  of  a parabolic  airfoil  we  try 

F(»,p)  = ^|o*  F(t)c‘*R(£  - r,V  - T,s,p)dr,  (5.3) 

where  R is  now  given  by  (4.8)  and  this  can  be  directly  verified  to  satisfy  the  governing 
equation  (3.8a).  Also,  substituting  this  into  the  boundary  condition  (3.8b)  and  re- 
calling that  it  is  valid  to  0(dQ),  it  can  be  shown  by  using  some  of  the  results  obtained 
by  the  present  author  (1974)  that  the  boundary  condition  is  indeed  satisfied  to  the 
same  order,  the  details  being  given  in  appendix  A.  From  (5.3),  F obviously  vanishes 
for  s < 0.  Hence  (3.8c)  is  satisfied  and  (5.3)  is  in  fact  the  solution  sought.  Before  we 
write  down  the  final  solution  explicitly,  we  restore,  in  order  to  obtain  in  (3.7),  the 
exponential  factor,  which  may  be  written  to  the  present  order  of  approximation  as 

ex?  [*"  JT  Ar(j3-*)0'a,(«)]  - ®xp  [i^  iV2roy0'<»(s)j . 

When  we  collect  all  the  results  obtained  so  far,  we  have  the  following:  if  the  airfoil 
shape  in  the  mean  position  is  given  by 

ef(x)  = e(\ax*+Pz), 

where  a < 0 (a  convex  surface),  and  the  co-ordinate  of  the  moving  upper  surface  is 
given  by 

y = ef(x)-d9eMg(x), 

where  the  amplitude  of  the  motion  g(x)  is  an  arbitrary  function  of  x,  then  the  leading 
term  of  the  unsteady  part  of  the  velocity  potential,  ^<u  in  (3.2b),  becomes 


^(s.p) 


xexp{»£r[(tf-l)  (*»  (•)-*»  (r)) + N(s  - r)  (•)]} 
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in  \ tt  //  v , ikUmmx  , w *r  y+l/if.'V* 

FM-tu-M+Tf r-«M.  "-  V(-S-)  • 


0(1)(«)  = (-C^/m)(J<M* +/?«), 

* - "sy » - *+-»-ifc 

|(e/m)(ow  + /S)|  « 1. 

This  integral  representation  is  the  solution  we  have  been  seeking. f (The  last  inequality 
is  a restriction  due  to  the  assumption  of  a small  perturbation.)  Before  attempting  to 
extract  physical  meanings,  we  pause  in  the  next  section  to  observe  that  the  present 
solution  embraces  the  various  known  results  as  special  limiting  cases. 


6.  Limiting  cases 


6.1.  Steady  limit 

In  the  limit  <o  ->  0 or  k-*0,  from  the  limiting  form  M(a,  6,0)  = 1 of  the  confluent  hyper- 
geometric  function  (e.g.  Abramowitz  & Stegun  1964,  p.  108),  (6.4)  is  immediately 
reduced  to 


H(s)  . 
~ m 


(6.1) 


from  (3.6 b).  This  is  Whitham’s  rule  for  steady  flow  and  becomes  identical  to  (3.4)  if 
we  replace  / by  -g.  We  wish  to  emphasize  that  g is  an  arbitrary  function  and  that 
we  have  recovered  the  above  as  the  limit  for  zero  frequency  of  oscillation. 


6.2.  Oscillating  flat-plate  air/oil 

In  the  limit  e->0,  when  we  note  that  (Abramowitz  & Stegun  1964,  p.  606) 

lim  M(a,  1,  -z/o)  = ./0(22*), 

a—*  oo 

(6.4)  becomes  at  once 

= H{xrmy]  V{T)e<*J,  {JL  [(*-T)«-m*y*]i)  dr,  (6.2) 

which  is  precisely  the  well-known  flat-plate  solution ; the  physical  meaning  of  this 
integral  representation  was  given  by  the  present  author  (1974). 

t In  this  connexion,  it  is  of  interest  to  note  that  Goldstein  ft  Rioe  (1973)  found  a solution 
for  sound  propagating  through  a uniform  shear  flow  in  terms  of  the  parabolic  cylinder  function, 
which  is  intimately  connected  with  the  confluent  hypergeometrio  function. 
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6.3.  Oscillating  wedge 

The  third  case  which  invites  comparison  with  the  present  result  is  that  of  an  oscillating 
wedge.  In  the  limit  a -►  0,  with  the  aid  of  the  limiting  formula  cited  in  § 6.2  we  obtain 

Mi)  — f y (T)  eikr  eXp  ( _ 2iNeSlcy)  exp  [ — — e/3k(2N  — 1 ) (s  - r )|  (6.3) 

m Jo  [ m J 

* J°  (f*  ^ “ T)  ^ " T)^)  dr' 

where  s = x — my  + 3f*,  Ne/iy, 

p = x + my—(N—2)M%e/Jy 

and  efi  is  the  semi-vertex  angle  of  the  wedge. 

In  order  to  confirm  the  agreement  of  this  formula  with  that  obtained  by  previous 
workers,  we  first  restore  the  factor  e~ikx  to  (6.3).  It  is  convenient  to  rotate  the  co- 
ordinate system  from  (x,  y)  to  (xt,  y2),  where  xt  is  parallel  to  the  upper  surface  of  the 
wedge  and  yt  normal  to  it.  At  the  same  time,  we  refer  the  flow  properties  to  the  mean 
steady  flow  behind  the  shock  instead  of  those  upstream  of  the  shock  and  designate 
them  by  a subscript  2.  Furthermore,  we  change  the  integration  variable  from  t to 
y = r(l  +mieft).  All  this  transforms  the  right-hand  side  of  (6.3),  upon  discarding 
negligible  quantities,  into  the  following  expression: 

I Czt-mJIt 

60e~ikx)/r<l)  ~ 0a—  H{xt - mtyt)  V(y) exp  ( - iktxt) exp (ikty) 

mt  Jo 

xJ°  (6-4) 

This  is  identical  to  the  flat-plate  solution  (6.2)  if  the  latter  is  expressed  in  terms  of  the 
(xvyt)  co-ordinate  system  and  the  flow  properties  downstream  of  the  shock.  This 
result  is  not  unexpected,  since  it  is  known  that,  if  one  takes  the  second-order  equation 
for  the  unsteady  component  of  the  velocity  potential  to  O(e0o)  and  expresses  it  in  terms 
of  these  co-ordinate  systems  and  flow  variables,  then  for  a wedge  it  exactly  reduces 
to  the  acoustic  equation.  (The  reason  why  the  relationship  (6.4)  is  approximate  rather 
than  exact  is  obviously  due  to  the  fact  that,  in  the  course  of  applying  the  strained 
co-ordinate  technique,  some  non-essential  second-order  terms  have  been  discarded.) 

Carrier  (1949)  obtained  a solution  for  a wedge  oscillating  at  its  apex;  the  solution 
was  derived  in  a more  generalized  way  by  including  the  rippling  motion  of  the  shock 
and,  in  addition  to  the  irrotational  component  of  the  flow,  rotational  flow  behind  the 
shock.  His  solution  was  later  generalized  to  include  the  case  of  a moving  vertex  by 
Van  Dyke  (19636),  who  also  corrected  typographical  errors  in  Carrier’s  paper.  The 
solution  was  expressed  in  the  form  of  a series  involving  Bessel  functions.  In  order  to 
facilitate  direct  comparison,  we  recast  the  present  solution  (6.4)  in  the  following 
alternative  form: 

- -0,0^2  (xS-(m1y,)t]i|  exp (-»*,*,),  (6.6a) 
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tanhd  = mtyjxt, 

K = (iMtvlktmt)  [«*+(-  <)"'] + 6 cos  (c#  [r  - ( - *)—].  (6.66) 

t = i(Mt  + mt) 

and  where  V (xt)  = Ut  + it o(xt  - ft  cos  e/9),  ft  being  the  pivotal  position  of  the  oscillating 
wedge  measured  from  the  apex.  The  above  identity  is  given  in  appendix  B.  Now 
Carrier’s  solution  for  the  irrotational  component  of  the  flow  becomes,  in  the  present 
notation, 

90at  2 [<*„  cosh  vd+b„ sinh  vd]  J„  [xj  - (to, y,)*]ij  exp ( - ihtxt).  (6.6) 

(The  expression  for  b,  given  in  (6.56)  is  the  corrected  one  given  by  Van  Dyke  1953 b.) 
Carrier  showed  that  as  long  as  the  shock  is  sufficiently  weak 

a,  * -M 

and  in  such  a case  (6.6)  is  indeed  identical  to  the  right-hand  side  of  (6.6a).  This  agree- 
ment naturally  endorses  the  present  viewpoint  that  the  global  behaviour  of  the 
unsteady  flow  downstream  of  the  weak  bow  shock  can  be  determined  essentially 
independently  of  the  presence  and  movement  of  the  shock. 


7.  Alternative  representation  of  the  solution  and  interpretation 

Returning  now  to  the  immediate  subject  of  a parabolically  curved  airfoil,  the 
solution  as  given  in  (6.4)  is  not  appropriate  for  extracting  its  physical  significance. 
Such  an  interpretation  will,  however,  be  obvious  once  we  recast  (6.4)  in  a more  reveal- 
ing form  by  making  use  of  the  following  Tricomi  (1949)  expansion  formula  for  the 
confluent  hypergeometric  function  in  a series  of  Bessel  functions: 

Jlf(a,6,x)  = r(6)(Ax)i<1-ft>exp(Jx)^^B^j*  Jn+ft_j[2(Ax)*]  for  Reft  > 0,  (7.1) 

where  A is  the  Whittaker  parameter,  given  by  A = J - a,  and. 

A0  = 1,  Ax  = 0,  At  = 


(n  + 2)A„+l  = (n  + I)i4„-2A^4B_1. 
When  we  insert  this  into  (5.4),  we  obtain 


^ f V(t)  e***  exp  ( ickfr) 
m jo 

x{.'.(3£[(*-r,(3>-r)Jl) 

+ 2 An\- 

n— 1 L 

(NeaMA*  „ .1*- 

-(„)(.  r )lr  ,)j 

| dr, 

(7.8a) 

t Van  Qyke  (19536,  also  private  communication)  proved  that  for  a small  wedge  angle 

aj 6,  a -l+SfW.efo  + Ote*). 

As  for  the  rotational  component  of  the  flow,  the  first  term  of  its  series  representation  may  be 
shown  to  be  O(0tA). 


: 

if 


I 


i 

i 


*&**--•■ ' 


where 
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At~\,  ^s  = -!A.  it4-»,  (7.26) 

(»+  1M,,+1  = nAH_1~2AAH_t,  A = mkiiiM^Nea)-1, 

<r  = -2 Nmy±  ^'<»(s)-(N- 1)^-  [#«(«) -^a)(r)]-^(e-r)  i-  0'<»(s) 

aN  * 
--{s-r){P-r) 

(AH  three  of  the  limiting  cases  of  the  preceding  section  are  now  directly  derivable  from 
the  present  form:  for  example,  when  k is  set  equal  to  zero,  the  result  (6.1)  follows  at 
once.)  Equation  (7.2a)  immediately  surrenders  itself  to  the  following  physical  inter- 
pretation. Let  us  first  examine  the  flow  field  near  the  leading  edge,  where  both  y and  s 
are  small.  Then  (7.2a)  becomes,  approximately, 

Vl)  = H^nm—  Jo  y V(T) eikr  Jo  (-g-  [(* - t)* - (my)*]ij  dr. 

This  is  the  flat-plate  solution  (6.2),  and  in  this  region  the  effect  of  the  body  shape  is 
indiscernible  as  yet;  the  unsteady  flow  field  is  completely  separated  from  the  non- 
uniform,  steady  flow.  Physically  the  decoupling  occurs  because  the  unsteady  dis- 
turbance, having  travelled  only  a short  distance  from  the  leading  edge,  has  suffered 
little  distortion. 

We  now  move  away  from  the  leading  edge  by  increasing  the  value  of  y while  keeping 
the  value  of  s constant  (along  the  straight  Mach  wave)  or  penetrate  downstream  by 
increasing  the  value  of  s while  keeping  y constant.  In  either  case,  if  we  look  at  the 
integrand  of  (7.2a)  or  the  signal  emitted  at  a point  t on  the  airfoil,  the  complex 
exponential  term,  which  can  be  written  as 

exp  ( - iek<r)  = exp  [Nmy{$'a)  (s)  + (r))  + (2JV  - 1 ) (0<»  (s)  - 0<l>  (t))]  j , 

immediately  discloses  the  following  key  aspect:  no  matter  how  slender  (e  1)  the 
airfoil  may  be,  this  phase  shift  (induced  by  the  presence  of  the  body)  will  eventually 
amount  to  an  increasing  delay  at  a large  distance  y or  s.  Moreover,  it  is  also  crucial 
to  recognize  here  that  the  phase  lag  of  the  signal  received  at  a position  s depends  not 
only  on  the  local  flow  at  that  point,  but  also,  through  the  very  difference  in  the  Bteady 
velocity  potential,  i.e.  $a)(s)-<friV(r),  upon  the  entire  flow  field  which  the  signal  has 
traversed;  the  disturbance  ‘remembers’  its  past.  Thus  we  might  call  this  exponential 
factor  the  phase  memory,  a term  commonly  used  in  connexion  with  the  propagation 
of  a radio  wave  through  a stratified  ionosphere  (e.g.  Budden  1961).  As  stated  in  the 
introduction,  the  existence  of  phase  memory,  which  differs  from  one  signal  to  another, 
is  by  itself  quite  sufficient  to  induce,  upon  superposition,  a change  in  the  amplitude 
of  the  unsteady  flow  field.  The  change  is,  however,  further  enhanced  because  the  shape 
of  the  airfoil  alters  even  the  amplitude  of  the  individual  signal  in  the  far  field  when  the 
contributions  from  the  higher-order  terms  of  (7.2a)  in  the  series  of  Bessel  functions 
begin  to  surface.  Thus,  in  the  far  field  the  airfoil  shape,  in  its  effect  of  causing  non- 
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R (parabolic)//?  (flat  plate) 


Figure  4.  Amplitude  and  phase  of  unsteady  pressure ; — p'(pv  U%  e'"1)-1  = Re1*.  The  ordinate  for 
the  left  figure  of  each  pair  is  the  ratio  of  the  amplitude  R for  a parabolic-arc  airfoil  to  that  for  a 
flat  plate.  The  ordinate  on  the  right  is  the  difference  in  phase  <f>\  = 1-3,  y = 1-4,  e = 0*1, 

a = —1,  /?  = 0-5  (max  y/c  = 0-0125)  and  the  pivot  axis  is  at  the  leading  edge,  y/c  = 

e[t<x (z/e)*+/S(x/e)].  (a)  wc!Um=  0-1  (kc  = 0-245) : ,y/c  = 0 [ky  = 0); — ,y/c  = 0-82(ty  = 0-2) : 

, y/e  = 2-04  [ky  = 0-5).  (6)  wc/C/„  = 0-3  (fcc  = 0-735) : , y/c  = 0 (Jfcy  = 0); — , y/c  = 0-82 

(iy  = 0-8);  , y/c  = 2-04  (ky  = 1-5).  (c)  htclUm  = 1 (kc  = 2-45):  , y/c  = 0 (ty  = 0) : 

— .y/c  = 0 82  [ky  = 2);  ---,y/c  = 2-04  (fcy  = 5). 

uniform  surrounding  flow,  is  inextricable  from  the  unsteady  flow  field  and  deeply 
affects  both  its  phase  and  amplitude,  as  well  as  the  directional  change  in  the  charac- 
teristic curves. 

This  point  is  illustrated  in  figure  4,  where  the  unsteady  pressure  distribution  for  a 
parabolic  airfoil  (max  y/c  = 0-0125),  computed  from  (7.2),f  is  compared  with  the 
result  for  a flat-plate  airfoil  at  three  different  frequencies  of  oscillation:  wc/t/*  — 0-1 
in  figure  4(a),  (ocjUm  = 0-3  in  figure  4(6)  and  ok/U^  = 1 in  figure  4(c).  There,  both  the 
amplitude  R and  phase  4>  are  plotted  as  functions  of  «,  i.e.  the  distance  between  the 
root  of  a straight  Mach  wave  and  the  leading  edge,  and  at  three  different  values  of  y. 
(If  the  flow  were  steady  then,  regardless  of  y,  the  amplitude  would  remain  the  same 
along  the  characteristics  » = constant.)  We  observe  that,  though  for  (oc/U^  — 0-1  the 
effect  of  the  airfoil  shape  does  not  become  prominent  at  these  values  of  y,  it  begins 
f For  numerical  computations,  (7.2)  is  also  more  convenient  than  (5.4). 
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to  emerge  at  (oc/U^  = 0*3 ; and  for  atc/U^  = 1 , except  for  the  close  vicinity  of  the  lead- 
ing edge,  it  indeed  alters  the  pressure  distribution  substantially. 

The  profound  modification  of  the  unsteady  linear  theory  displayed  here  raises  an 
obviously  disquieting  thought  on  the  upshot  of  the  acoustic  theory  when  multi-body 
aerodynamic  interference  is  involved  and  deepens  concern  expressed  (Kurosaka  1975) 
with  regard  to  some  of  the  consequences  arising  from  a pro  forma  sum  of  linearized 
unsteady  up  washes. 

8.  Further  interpretation 

Pursuing  the  physical  interpretation  further,  we  seek  the  connexion  between  (7.2a) 
and  the  ‘fundamental’  solution.  We  shall  not,  however,  merely  reconstruct  (7.2a)  by 
the  superposition  of  the  fundamental  solution.  Rather,  we  shall  reverse  the  usual 
process  and  obtain  the  fundamental  solution  from  (7.2a):  that  is  to  say,  we  regard 
(7.2)  as  the  spectrum  at  frequency  <o  or  the  Fourier  transform  and  take  its  inverse 
transform  so  as  to  derive  the  transient  response  to  an  arbitrary  time-dependent  motion 
of  the  airfoil.  The  ‘ fundamental  ’ solution  will  arise  naturally  in  the  course  of  obtaining 
the  transient  response  (Miles  1959,  p.  53).  Let  us  go  back  to  (2.4)  and  rewrite  the 
unsteady  part  in  a more  general  way  as 

<H  = efl+0o£l(z,y:t). 

Then  the  Fourier  transform  fi(w)  of  O (its  leading  part)  is  equal  to  e~ikxi/^l\  \j^l)  being 
given  by  (7.2a),  provided  that  V is  regarded  as  the  Fourier  transform  P of  itself,  i.e. 


An(<o)cnJH  [(«— r)(p— T)]i|J  dr, 

where  c,  - [ - (^>f°°)  (»  - r)  (p  - r)j  ** . 

Taking  the  inverse  transform 

we  obtain,  by  convolution, 

£l(z,y,t)  - 1JJ  drj^  V(T,t-l)F(i)dg.  (8.1a) 

midy}^  = V(z,t) 


) 

! 

j 

i 

i 
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and  where 
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a = (x-r  + eoh  b = -±—  [(«-T)(p-r)]*, 

axm*  awmM 

cos  6 = a/b,  r = (6*  - a*)i. 

We  note  that  Cl(x,y,t)  contains,  through  Fn  in  (8.16),  the  term  1/r,  which  can  be 
written  as 

^ = a<Bmji^-t{8-T)(p-T)-^aKmi-^  (x-r  + ea)^  j*.  (8.2) 

The  meaning  will  become  immediately  recognizable  if  we  note  that  at  e = 0 the 
denominator  of  (8.2)  may  be  reduced,  after  some  algebra,  to 

{(a„Z)*-(x-T-U„£)*-y*}k. 

This  represents,  when  set  equal  to  zero,  a circular  wave  front  of  a disturbance  which 
was  emitted  at  a source  point  (r,  0)  and  is  propagating  through  uniform  flow  after  a 
time  £.  Thus  the  denominator  of  (8.2),  when  put  equal  to  zero,  i.e. 

1 r M 1* 

— (s-T)(p-T)-\a„m£--2(x-r  + e<T)\  =0,  (8.3) 

m1  L m J 

now  describes  the  distorted  wave  front  propagating  in  a non-uniform  flow  field.  In  fact, 
we  can  directly  show  that  the  expression  for  £ obtained  from  (8.3)  does  satisfy,  within 
the  approximation  consistent  with  the  present  analysis,  the  appropriate  eikonal 
equation  at  large  distances;  (8.2)  is  indeed  the  fundamental  solution.  In  general,  for 
a given  point  (x,  y)  in  flow  and  for  a given  source  point  (r,  0),  there  are  two  values  of  £ 
satisfying  (8.3):  one  corresponds  to  the  time  when  the  disturbance  first  arrives  at 
(x,  y)  and  the  other  to  the  time  when  it  departs  from  (ar,  y).  In  the  particular  case  when 
the  point  (x,  y)  is  located  such  that  either 

s = r or  p = T, 

there  is  only  one  such  moment  for  £,  which  implies  that  the  wave  front  is  tangential 
to  either  s = t or  p = r.  a = r corresponds  to  the  straight  Mach  wave,  whose  root  is 
located  at  (t,  0) ; p = r is  the  cross  Mach  wave  passing  through  the  same  point.  Hence, 
as  expected,  two  families  of  Mach  waves  passing  the  source  point  form  envelopes  for 
the  disturbance  emitted  from  the  source.  In  particular,  the  time  required  for  the  signal 
to  arrive  at  a point  on  the  straight  Mach  line  s = r is  given  by 

i,A) 

It  is  of  interest  to  note  that  this  can  be  obtained  in  the  following,  more  physical  way. 
The  wave-front  velocity  c is  in  general  the  vectorial  sum  of  the  local  acoustic  speed 
in  the  direction  of  the  normal  n to  the  front  and  the  convective  fluid  velocity,  i.e. 
c «■  an  + u.  However,  along  the  enveloping  Mach  waves,  which  are  tangential  to 
the  wave  front,  the  acoustic  speed  does  not  contribute  to  the  component  of  the  wave- 
front  velocity  in  the  direction  parallel  to  the  Mach  wave ; only  the  fluid  velocity  con- 
tributes. In  particular,  along  the  straight  Mach  line  the  component  of  the  fluid 
velocity  or  the  wave-front  velocity  remains  constant.  If  we  divide  the  distance  from 
the  source  (r,  0)  to  the  point  (x,  y)  by  the  component  of  the  flow  velocity  in  the 
direction  of  the  straight  Mach  wave,  we  can  directly  derive  (8.4),  as  the  time  elapsed. 
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9.  Concluding  remarks 

It  has  been  our  aim  to  find  a uniformly  valid  solution  for  the  unsteady  flow  field  and 
examine  it  in  detail.  We  have  shown,  through  an  explicit  solution  obtained  for  the 
specific  case  of  a parabolic-arc  airfoil  oscillating  in  supersonic  flow,  that  the  prediction 
of  the  unsteady  signal  in  the  far  field  demands  the  detailed  description  of  the  contour 
of  the  moving  boundary.  The  non-uniform  surrounding  flow  produced  by  the  very 
presence  of  the  body,  no  matter  how  slender  it  may  be,  cumulatively  and  inextricably 
affects  both  the  amplitude  and  the  phase  of  the  unsteady  disturbance  at  a large 
distance  from  the  leading  edge. 

As  a further  related  effort,  it  would  appear  to  be  worth  while  to  pursue  a study  for 
other  airfoil  shapes  so  as  to  enlarge  our  stock  of  particular  solutions.  With  regard  to 
the  question  of  similar  cumulative,  first-order  effects  of  nonlinearity  in  subsonic  flow, 
we  still  remain  uncertain.  It  is  intriguing,  however,  to  note  that  in  a very  recent  paper 
of  Goldstein  & Atassi  (1976),  where  an  exact  second-order  solution  is  obtained  for  an 
airfoil  subject  to  a convected  gust,  the  incoming  gust,  in  its  nonlinear  interaction  with 
the  steady  non-uniform  flow  field,  is  found  to  suffer  distortion  in  wavelength  in  a 
manner  akin  to  the  present  supersonic  result  though  the  flow  treated  there  is 
incompressible. 

The  author  would  like  to  extend  his  gratitude  to  Mr  C.  E.  Danforth  for  calling  his 
attention  to  the  problem,  to  Mr  I.  H.  Edelfelt  for  assisting  in  numerical  calculations, 
to  Dr  S.  D.  Savkar  for  his  many  helpful  suggestions  and  to  Dr  M.  E.  Goldstein,  Dr 
L.  J.  Slater  and  Professor  M.Van  Dyke  for  answering  his  queries.  He  is  especially 
indebted  to  Professor  Van  Dyke  for  generously  making  available  details  of  his  earlier 
notes.  The  work  was  supported  by  the  Air  Force  Office  of  Scientific  Research  under 
Contract  No.  F44820-74-C-0040. 

Appendix  A 

In  this  appendix  we  shall  show  that  the  expression  for  F given  by  (5.3)  does  satisfy 
the  boundary  condition  (3.86)  to  order  60.  We  denote  the  left-hand  side  of  (3.86)  by 

1{F)  *<?0{iefc  ["/'<«)-£  **'<»(.)]  F + F.-F,-!  F(s)e**),  (Al) 

and  we  shall  prove  that  along  a - p this  vanishes,  to  order  0O.  Substituting  (5.3)  into 
the  above,  one  obtains  for  a > 0 

x{«  [-;/'<•>- ErV+Er  (•’.-*',>]  Jf+Mt-Jy)*,  <AJ) 

where 
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We  observe  that  along  s -p 


v,  = (2N-l)<fi’(»(s)-N(T-s)t'<»,  pp  = 0,  Mt  = Mp.  (A3) 

Furthermore,  by  Tricomi’s  expansion  formula  cited  in  §7,  the  confluent  hyper- 
geometric function  M can  be  expressed  as 

M = exp  |-2^-  — t)V*<1,J 

<A4) 

where 

A0  = 1,  Aj  = 0,  At  = \ 

and  the  other,  higher-order  An  are  the  same  as  those  given  in  (7.26).  The  leading  term 
of  31  is  given  by 

M ~ exp  [ - iekN(s  - r)*  0*»>]  J0  [A-  (s  - r)] . (A  6) 

From  (A  3)  and  (A  5),  (A  2)  becomes 

1(F)  = 6gek  Jo  V*(t ) exp  (i/fr)  exp  rj  exp  J^-~  iV(s-r)*9i'(l,j 

*J*[wJ8-T)]dT)'  (A6) 

where 

V*(r)  = - F(r)»  U~JN(r-s)<f"a>. 

Equation  (A  6)  can  contribute  to  O(0o)  only  when  ek  is  such  that,  if  properly  non- 
dimensionalized,  0(ek)  = 1 or  k = 0(  1/e).  For  such  large  values  of  k,  we  apply  the 
following  method  of  obtaining  an  asymptotic  expansion  (Kurosaka  1974):  we  first 
express  J0  in  terms  of  an  integral  involving  an  exponential  and  use  the  stationary- 
phase  method  repeatedly.  This  yields 


and  (A  6)  becomes 


1(F)  - 0(0oe), 

which  is  of  higher  order  than  O(90) ; the  other  terms  of  (A  4)  may  similarly  be  shown  to 
be  of  higher  order.  Hence  to  O(0o),  1(F)  «=  0. 


Appendix  B 

In  this  appendix  we  shall  prove  the  identity  (6.5a): 

If  exp  ®xp  ( “ *****>  (~Mt  “ *)*  ~ (”*1?*)*]*)  drl 

“-«»*•  2 b,  t~*J,  (|l  [xj  - (mjy,)*)*!  exp ( - »*,*,),  (B 1 ) 
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xt  > m,y„  tanh  6 = b,  = *7^  [<*’+(-<)']  + 6 cos  (e/J)  [«'-(-  <)-'], 

X|  Km  7llj 

1 = + m,),  F(x,)  = t/,  + ko(x,-6cose/?).  (B2) 

First  we  write 

V(if)exp(ikt7i)  = ^ j^Jf,-^cose/jj  [yexp(*'i:,y)j  + ^»  [y*exp(i*,y)]. 

As  suggested  by  Carrier  (1949),  we  expand  7exp(t£,7)  and  i)*exp(iktii)  in  series  of 
Bessel  functions,  through  the  generating  function  of  the  Bessel  functions,  and  obtain 

F(7)exp(ti,7)  = 7) • 

Substitution  of  this  into  the  left-hand  side  of  (B  1 ) yields 


-at00  2 b,vexp(-iktxt)F(xt), 


where 


► -1 


(B  3) 


,F(xt)  = [|jr  [(^»-7),-(»»tyl),]*]j^  Jr  (B4) 

If  we  take  the  Laplace  transform  ^*of  F(xt),  defined  by 

exp  ( - sx,)  F(x,)  rfx„ 


then,  by  convolution,  we  obtain 
& 


Inverting  this  gives  (e.g.  Erd&yi  et  at.  1954,  p.  250) 

^(*1)  = l M ~ (roiFi)1]*)  (B  5) 


for  xt  > TO,y,.  By  substituting  (B  5)  into  (B  3) , one  may  establish  the  required  identity 

(Bl). 
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"The  Effect  of  Airfoil  Contour  upon  the  Unsteady 
Aerodynamics  of  Supersonic  Cascades " 


The  Effect  of  Airfoil  Contour  upon  the  Unsteady 
Aerodynamics  of  Supersonic  Cascades 


M.  Kurosaka * 
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Tullahoma , Tennessee  37388 

and 

I.  H.  Edelfelt' 
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1.  Introduction. 

It  is  now  an  indisputable  fact  that  a subtle  change  in  airfoil 
contour  induces  significant  modification  upon  the  unsteady  aerodynamics 
of  supersonic  cascades  and  consequently  results  in  substantial  shift  in 
its  flutter  boundaries.  Its  phenomological  description  is  briefly  given 
in  Ref.  1 along  with  physical  explanation  based  upon  the  concept  of 
'phase  memory';  for  an  isolated  airfoil  of  parabolic  contour,  a detailed 
analysis  substantiating  this  concept  is  presented  in  Ref.  2. , where  the 
profound  effect  of  airfoil  contour  upon  the  unsteady  far  field  around  an 
oscillating  airfoil  is  in  fact  quantified. 
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The  objective  of  the  present  paper  is  to  extend  the  isolated 
airfoil  analysis  of  Ref.  2 to  cascade  arrangement;  in  doing  so,  we  shall  ‘ 

rely  heavily  on  the  formulation  and  results  of  Ref.  2,  particularly  in 

1 

the  usage  of  Riemann  function.  More  specifically,  in  analyzing  the  so-  j 

3 

called  pre-interference  zone,  which  is  invariably  the  most  difficult  portion  j 

| | 

of  the  problem  on  hand,  we  shall  make  exclusive  use  of  Riemann's  integral  j 

4 j i 

representation  of  a solution  for  hyperbolic  equations  ; applying  it  to  an  j 

appropriate  contour  near  the  entrance  section  of  cascades  adopting  the  ' j 

2 5 

passage  approach  and  with  the  aid  of  cascade  periodicity  condition,  we  j 

j 

shall  be  able  to  determine  the  entire  flow  field  in  the  pre-interference  j 

1 

•r 

zone  . 

Based  upon  these  analysis,  our  results  indeed  confirm  the  significant 
effect  of  airfoil  contour  upon  unsteady  pressure  distribution  acting 
on  the  surface  of  cascaded  airfoils. 


3 


2.  Formulation. 

Consider  oscillating  airfoils  in  a cascade  where  the  coordinate  of 
moving  surface  of  the  n-th  airfoil  is  expressed  by 

y » ef(x)  - eQ  eiUlt  einyg(x)  , (1) 

and  the  coordinate  system  (x,y)  is  shown  in  Figure  1;  in  the  above 
£f(x)  (£  « 1)  designates  the  shape  of  the  airfoil  in  its  mean  position 
of  oscillation  and  the  second  term  represents  its  harmonic  motion  with 
frequency,  w and  interblade  phase  angle,  U ; here  we  shall  restrict 
ourselves  to  parabolically  shaped  airfoil  where  f (x)  is  a quadratic  of  x. 

The  two  small  non-dimensional  parameters,  e and  0^  characterize  the  slender- 
ness of  the  airfoil  and  the  amplitude  of  motion,  respectively.  Far 
upstream  of  the  cascade,  the  flow  velocity  is  equal  to  and  the  Mach 
number,  Following  Ref.  2,  we  separate  the  perturbed  velocity  potential 
into  a steady  and  time-dependent  part  by  writing 

$ - e(|)(1)  + 0Q  exp  [ i (cot  - kx)]  ^(1),  (2) 

2 2 

where  k ■ w/m  U^;  the  first  term  represents  the  steady,  base  flow 
and  che  second  term  corresponds  to  the  unsteady  flow;  our  present  interest 
is,  of  course,  focused  on  tj/.  It  turns  out  to  be  convenient  to  introduce 
the  characteristic  coordinate  (s,p)  defined  by 


4 


1 

1 2 2 2 2 
where  N ■ j (T  + D ^ /«  . and  m = (M^  - 1)  , y is  the  adiabatic 

exponent  of  the  gas;  instead  of  the  Cartesian  coordinate  system,  we 

shall  hereafter  use  the  above  characteristic  coordinates  as  the  independent 

variables.  The  steady  part  of  the  velocity  potential,  <t>^\  is  given  by 

(1)  U» 

- -H  (s)  — f (s)  , (4) 

where  H is  the  step  function.  With  regard  to  the  unsteady  part 
upon  introduction  of  the  function  F defined  by, 

¥(1)  - exp  [i(£k/UJ  N (p  - s)  $,(1)(s)  F],  (5) 

it  is  found  to  satisfy  the  following: 

F + ^ t-  <2N  - 1)  <J>'(1)(s)  + (p  - S)  N 4>"(1)(s)]  F 

« oo  * 

2 

+ [fcg£>  + N <j>' ' (1)  (s)  ] F * 0,  (6) 

with  the  boundary  conditions  on  the  airfoil  surface 

iek  l"f'(s)  - f N $’(1)(s)]F  + F - F - ^ V(s)  eikB  - 0 , (7) 

m U s p m 

00  r 

where 

V(x)  « [U^g'Cx)  + (ik  U^/M^2)  g(x)]  einU.  (8) 


The  initial  condition  is  that  at  any  given  point,  say  Q^,  of  Figure  1 
just  upstream  of  the  shock  emanating  from  the  leading  edge  of  the  reference 
airfoil,  the  unsteady  velocity  potential  there  must  be  equal  to  the  one 
at  the  corresponding  point,  Qj,  upstream  of  the  leading  edge  shock  of 
the  following  airfoil,  with  the  exception  of  the  phase  lag,  i.e.. 


LA 
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3.  Solution. 

(2) 

Our  present  aim  Is  to  determine  F , which  satisfies  the  governing 
equation  (6)  along  with  the  following  boundary  condition 


iek  [-if(s)  - -f  N 

00 


(1) 


(•)]  F 


(2) 


+ F 


(2) 


(2) 


(12a) 


and  the  initial  condition 


[e  ikx  f(2)], 


[e"ikx(F(1)  + F(2))], 


-ip 


(12b) 


The  governing  equation  (6)  is  hyperbolic;  it  is  well  known  that  the  solution 

of  any  second-order  linear  hyperbolic  equation  can  be  expressed  in  the  form 

of  integral  representation,  once  the  corresponding  Riemann  function  R is 
4 

obtained  . More  specifically,  for  any  point C of  Figure  1,  whose  character- 
istic coordinates  are  given  to  be  (s,p),  the  Riemann's  integral  representa- 
tion can  be  expressed  as 


F(2)(s,p)  - | [F(2)(A)  R (A)  + F(2)  (B)  R (B) ] 


[-Fn(2)R  + F^2^]  dr,  + [F_(2* 


+ 2bF 


((2)  R + f(2)  r^]  dn 


[F  (2)R  -F(2)R  + 2bF(2^R]tfC 


OA 


(13a) 


where,  as  shown  in  Figure  1,  A and  B are  intercepts  of  the  characteristics 
passing  the  point  C and  the  airfoil  and  the  bow  shock,  respectively,  and  0 


7 


is  the  origin;  the  line  integrals  should  be  carried  out  in  the  sense 
indicated  above,  that  is,  in  the  first  integral  from  point  B to  point  0 
and  in  the  second  integral,  from  point  0 to  point  A.  Furthermore,  in  (13) 
b is  equal  to 


b = ~r  t-  (2N  - l)  4>' (1) (£)  + (n  - C)  N <t(1)"(?)i. 


(13b) 


and  R,  the  Riemann  function,  is  given  by 


R(5»  n.s,  p)  = expv'  M (a,l,z) , 


(13c) 


where 


= iekTj-  ( (N  - 1)  (<t(1)(s)  - <j>(1)(0)  " N [n  (0,(1)(8)-  <t>,(1)(  )) 


-(s4>(1)'(s)  - C4>(1),(0)]  } , 


(13d) 


and  where 


3 = 2 + 4ie  NM  2 ’ 


ekN  (s  - O (p  - n)  <t>’ 

Uoo 


(13e) 


(13f) 


We  rewrite  the  second  integral  appearing  in  (13a)  with  the  aid  of  (12a) 
and  (13b)  and,  noting  that  along  the  path  OA,  C ■ n , (13a)  becomes 


F(2)(s,p)  - i [F(2)(A)  R(A)  + F(2)(B)  R(B)] 


+ j j ]^-Fn(2)R  + P(2)Rn]  dn  + [F^(2)R  - F(2)R?  + 2bF(2)R]  dcj 
BO 


(14) 


f(2) 


We  then  replace,  from  (12b),  those  values  of  Fv  ' to  be  evaluated  along  the 

path  OB  located  just  upstream  of  the  first  shock  by  the  ones  at  the 

corresponding  points  just  upstream  of  the  second  shock.  By  this,  all 

the  values  of  F^  appearing  both  sides  of  the  equation  become  those 

(2) 

downstream  of  the  first  shock.  Consequently  F involved  there  does  not 
contain  any  discontinuity  and  regarding  the  resulting  expression  as  an 
integral  equation,  one  solves  it  by  the  collation  method.  This  determines 
F'2^  and  hence  we  obtain  the  flow  field  in  the  pre-interference  zone. 
Pursuant  to  this,  we  obtain  the  flow  field  between  the  blade  passage  by 
the  enperposition  principle. 


i 


9 


4.  Discussion  on  the  Results. 

The  computed  pressure  distributions  are  shown  in  Figure  2-5. 

In  all  of  the  figures,  the  shape  of  the  airfoil.  f(x)  of  (1)  is  as  follows: 

f(x)  - -1.7633  (x^  - xc)  , 

where  c is  the  airfoil  chord;  the  airfoil  is  executing  torsional  motion 
at  its  mid  chord,  i.e.  g(x)  in  (1),  is  given  by 

g(x)  = x - c/2. 

In  addition, 

H = 1. AO 
00 

Stagger  angle  - 65°, 

Solidity  - 1.589 
M, interblade  phase  angle  ■ 180° 

The  trajectory  of  the  shock  locus  at  its  steady  position,  which  is  needed 

in  solving  (14),  is  computed  according  to  the  routine  used  in  Ref.  5. 

Figure  2 shows  both  the  real  and  imaginary  part  of  the  pressure 

1 2 

coefficients,  defined  by  C ■ (p  - p^J/y  Poo^co  » an(*  co®Pare  the  present 
results  obtained  for  e ■ 0 (the  limits  of  zero  thickness  airfoil)  with  the 
ones  previously  computed  for  flat  plate  airfoils^;  as  expected,  both  results 
agree  completely.  Figures  3 to  5 compare  the  present  results  for  e ■ 0.1 
with  the  ones  for  flat  plate  airfoils  at  three  different  values  of  k 
and,  as  might  be  anticipated,  the  effect  of  airfoil  contour  indeed  induces 
significant  modification. 
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I.  Introduction. 

One  of  the  treacherous  flow-induced  vibration  problems  in  high-speed 

aircraft  engines  is  the  so-called  supersonic  flutter,  which,  if  allowed 

to  persist,  may  be  capable  of  inflicting  excessive  damage  to  their  structural 

Integrity.  For  more  detailed  information  on  this,  we  relegate  to  elsewhere*". 

Here  we  only  remark  that  flutter  takes  place  when  the  relative  velocity  at 

the  tip  exceeds  sonic  speed  and  at  the  time  it  occurs  near  the  operating 

line,  the  flow  remains  unstalled.  Because  of  its  critical  importance, 

there  have  in  recent  years  an  accelerated  activity  in  our  effort  to  under- 

2 

stand  and  quantify  the  phenomena  . However,  due  to  the  very  complexity 
of  the  problem  on  hand,  past  emphasis,  in  both  theoretical  and  experimental 
research,  has  been  laid  upon  two-dimensional,  cascade  problem.  Highly 
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useful  these  approaches  have  been  In  contributing  to  our  understanding 

of  the  phenomena,  there  always  been  a nagging  question  upon  the  validity 

of  the  two-dimensional  approximation  to  what  is  in  reality  three-dimensional 

flow  through  turbomachines.  The  question  is  certainly  a legitimate  one, 

since  even  in  the  steady  situation  where  bladings  do  not  vibrate,  it  has 

long  been  recognized  that  the  three-dimensionality  of  flow  induces  sub- 

3 4 5 

stantial  modification  to  cascade  representation  * ’ . If  once  the 

compressor  bladings  start  to  oscillate,  two  extra  factors  introduce 

additional  deviation,  both  attributing  to  influence  the  propagation  of 

disturbances  emanating  from  the  airfoils  oscillating  with  their  natural 
* 

frequency.  Firstly,  the  encasing  walls  serve  to  reflect  the  Incident 
walls  and  thus  affect  the  unsteady  flow  field;  secondly,  strong  radial 
gradient  in  steady,  base  flow  causes  refraction  of  transmitting  waves. 

The  full  detailed  analysis  of  three-dimensional  flow  through  turbomachinery 
being  dishearteningly  difficult  even  for  steady  flow,  the  satisfactory 

treatment  incorporating  the  above  features  in  addition  to  retaining 

other  necessary  constituents  such  as  contour  of  individual  airfoil  and 

geometry  associated  with  the  arrangement  of  bladings  appears,  in 

spite  of  some  recent  attempt^,  to  be  formidable.  Consequently,  instead 
of  attacking  the  subject  as  such,  we  turn  our  present  attention  toward 


Parentically,  if  the  natural  frequency  of  the  airfoil  Is  u and 
the  fundamental  blade  passing  frequency  is  it  can  be  easily  n 
shown  that  the  frequency  perceived  by  the  stationary  observer  is  equal  to 

«n  + (n  - u/2tt), 

where  n is  an  integer  and  y is  the  Interblade  phase  lag. 


the  following  problem,  which  seems  to  be  more  crucial  in  practical 
applications:  suppose  one  replaces  outer  casing  of  compressors  with 
material  which  is  capable  of  absorbing  unsteady  fluctuations,  what 
would  be  its  effect  on  compressor  flutter  boundary?  The  installation 
of  acoustic  lining  on  casing  walls  of  aircraft  engines  is,  of  course, 
widely  used  practice  to  alleviate  the  aeroacoustic  noise  but  its  potential 
benefit  upon  the  flutter  suppression  does  not  appear  to  be  exploited. 
However,  its  obvious  advantage  appears  to  be  worthy  of  exploration,  since 
owing  to  the  very  capacity  of  sound  absorbing  material  to  relieve  the 
flow  fluctuation,  the  unsteady  pressure  acting  upon  the  bladings  would 
directly  be  reduced.  Moreover,  in  supersonic  flutter  the  tip  portion 
of  the  bladings  is  the  most  critical  region  immersed  in  high  supersonic 
flow;  therefore,  the  provision  of  acoustic  liners  on  the  surface  of  outer 


casing  appears  — due  to  its  closest  vicinity  to  bladings  — to  produce 
an  inmediate  beneficial  effect  upon  the  tip.  Furthermore,  since  the 
Installation  of  lining  material  can  be  carried  out  entirely  independent  of 
the  aerodynamic  performance  of  bladings,  it  does  not  interfere  with  or 
compromise  the  other  various  design  consideration  of  compressors. 

To  appraise  this  concept,  in  the  present  paper  we  formulate  and 
analyze  a simplified  model  problem  where  an  isolated  oscillating  airfoil 
is  placed  in  a supersonic  duct  whose  upper  wall  is  lined  with  sound 


absorbent  material.  By  evaluating  the  influence  of  wall  liners  upon 
unsteady  pressure  received  at  points  off  the  surface  of  the  isolated 
airfoil,  we  assess  its  effect  in  turbomachines,  whose  individual  blading 
is,  of  course,  subject  to  far-field  upwash  generated  by  the  other 
members  of  airfoils. 
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II.  Model  Problem  Formulation. 

Consider  an  isolated  airfoil  placed  in  a two-dimensional  duct 
(Figure  1)  where  the  base,  steady  flow  moves  in  the  x direction  with 
supersonic  velocity,  U;  U is  assumed  to  be  uniform  across  the  duct  and 
there  is  no  variation  in  other,  steady-state  flow  variables.  The  air- 
foil executes  small-amplitude  harmonic  motion  in  the  transverse  z direction; 
one  surface  wall  of  the  duct  located  at  y • h is  provided  with  sound 
absorbent  material  having  specified  Impedance  while  the  other  wall  of 
the  duct  is  untreated.  If  one  represents  the  unsteady  quantities  by 
primes,  the  linearized  governing  equations  become 


(1) 


where  c is  the  unperturbed  speed  of  sound,  p the  ambient  density; 
u',  v',  w'  are  perturbed  velocities  in  the  x,  y and  z directions, 
respectively;  p'  is  the  perturbed  pressure  and  t is  the  time.  From 
all  of  these  above  perturbed  quantities,  we  separate  the  time- dependent 
factor  by  writing  p'  ■ p'e^U)t,  u'  ■ u,e*wt,  etc.  and  (1)  becomes 
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(2) 


The  boundary  condition  on  the  airfoil  is  that  the  z-component  of  fluid 
velocity  is  prescribed  and  equal  to,  say,  W (x,y) 


i 

i 

I 


w'  (x,y,  z « 0)  - W (x,y)  . 


(3) 


The  condition  on  the  hard  wall  is 


v'  (x,  y - 0,  z)  - 0 , 
and  on  the  soft-wall 


(A) 


(5) 


where  A is  the  specific  acoustic  admittance  and  is  assumed  to  be  constant. 
The  initial  conditions  are  such  that  upstream  of  the  leading  edge  Mach 
cone,  the  flow  remains  unperturbed  or  equivalently 


p'  ■ u'  ■ v’  ■ 


w'  - 0 


(6) 
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Because  of  the  supersonic  nature  of  the  flow,  we  focus  our  attention 
to  the  flow  for  z > 0 (above  airfoil),  which  is  uncoupled  from  the  flow 
for  z < 0 (below  airfoil) . 


7 

III.  Analutioal  Solution. 

In  order  Co  derive  an  analytical  solution,  we  take  the  Laplace 
transform  of  all  the  perturbed  quantities  in  the  x direction,  defined, 
for  example,  by 

00 

*'•  I «6X  — 

p - / e p*  (x.y.z)  dx, 

0 

and  obtain  the  following: 

-4jr  id)  p + Us  p + p (s  u + -|^  + “ 0 » (7a> 
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“ " p 
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w 
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3p 
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with  the  following  boundary  conditions: 
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0 , 


Pc  ~ 

" 411,1  ar 


y - h 


(8a) 

(8b) 

(8c) 
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Elimination  of  u,  v,  and  w from  the  above  equations  at  once  yields  the 
following  expression  for  p: 


32p  . 32o  ,22.  2iaUs  u)2  ~ n 

— 21  p " 0 ’ 

c c 


2 2 
3y  3z 


3 £_  r ‘ - a. 
-c—  - [ms  + 


with  the  boundary  conditions, 


3. 

3z 


-p  (iio  + Us)  w (y) , 


z ” 0 


y - 0 


y * h 


= 0 , 


c 3p 

A(io)  + Us)  3y 


y = h 


The  solution  for  p is  immediately  obtained  to  be 


w 

V’ 


n n 


P - > p (iu)  + Us)  e~anz  • 6,  cos  B„y, 

L- — )_  n 

n - 0 


where  Bn  is  given  as  solutions  of 


B h • tan  B h - - ■*-&?-+  Usi>h  , 


n n 


and  a is  related  to  0 by 
n n 7 


<*n2  - en2  - »2  t®2  + - 0, 

c c 


(9a) 


(9b) 


(9c) 


(9d) 


(10a) 


(10b) 


(10c) 
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where  S are  the  coefficients  associated  with  the  expansion  of  W(y)  in 
a 

the  series  of  cos  (Jny,  i.e.. 


W(y) 


V 


£ cos  6 y 
n n 


(10d) 


n • 0 


or 


h 

Sn  * [ J « 


(y)  cos  ($ny)  dy) /[  / (cos  B^y)  dy] 


(10e) 


If  we  denote  the  inverse  Laplace  transform  by  L , the  inverted  form  of 
(10a)  may  be  expressed  as 


P' 


p f ±ct^p1  + Up2)  , 


where 

P,  ■ L 1 t(~  e-ctnz)  (6  cos  B y)  ] , 
i a n n 

n 

P2  ■ L [(s  e anZ)  (£n  cos  6 n y) I - 
n 


(11a) 


(lib) 


(11c) 


In  order  to  evaluate  (lib)  through  convolution,  we  have  to  invert  the 
following  two  transforms,  respectively: 

e^n*  and  6 cos  8 y . 
a n n 

n 


10 


(The  other  transform  s/a  e 0nZ  appearing  in  (11c)  will  be  promptly 


n 


obtained,  once  1/a^  e anz  is  inverted,  as  described  shortly.)  The 


inversion  of  the  above  obviously  demands  the  explicit  expression  of 


3 and  in  order  to  obtain  it,  we  assume  that  the  specific  acoustic 
n 


admittance.  A,  is  small  and  then  (10b)  yields  approximately 


(12a) 


3 - ~ [1  - ±3—^-1 , n - 1,  2 

n h (2nir)2c 


(12b) 


A.  Inversion  of  — e~(XnZ . 


Using  these  expressions  and  with  the  aid  of  the  following  well-known 


inversion  formulae  , 


L1  { (s  + a)2  (s  + 3)  2 exp  [-b  (s  + a)2  (s  + 3)2] 


- H (x  - b)  exp  | (a  + 3)x]  JQ  [|  (a  - 3)  (x2  - b2)2]  , 
where  Jq  is  the  Bessel  function,  1/  an  e0*02  is  readily  inverted  to  be 


C1  [ft  “»'] 
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— H (x  - mz)  exp 
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. AU  . 

+ »]  x 

2hcnt 


x J 


{r  iuU 

l-  2 2 

c . 

{r_w_  i 1mA,.,  /_2  _2 _2\2  1 

[— + — J-12  (x  - ■ z ) / , 
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(13a) 
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and 
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1 „ , . / r icoU  . AU  , 

- — H (x  - m2)  exp  \ [ 2“ 2 + 2 x 

c in  hem 


„ 1 , (ji)^  , 2iwA  . /2nTrN^1^  2 2 2V  , 

x Jo  ^ C~TT  + ~4  + 1 (x  " m 2 > /» 

c m hem 


(13b) 


where  n ■ 1,  2,  . 

If  we  denote  the  right  hand  side  of  (13)  as  (x),  i.e., 


-1  1 - a„z 

L e n ] =»  6 (x)  , n « 0,  1, 

n 

_ a z 

then  the  Inversion  of  s/a  e appeared  in  (lie)  is  immediately  given  to 
be8 

-1  _ a z 

L [~  e n ] - 0'(x)  . (13c) 

n 


B.  Invention  of  bn  008 

In  order  to  invert  b cosg  y of  (lib),  we  recognize  that  if  the 

n n 

normal  fluid  velocity  W of  (3)  is  independent  of  y,  and  is  a function 
of  x only,  for  instance,  equal  to  q (x),  i.e., 

W(x,y,  * - 0)  - q (x), 
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then 

W (y)  - q, 

where  q is  a constant,  (10e)  yields,  to  the  present  order  of 
appr oximat ion , 

bQ  - q , (14a) 

£ * [iw  q + Us  q]  , (14b) 

n 2(nn)  c 

where  n * 1 , 2 . Making  use  of  these  and  noting 

-1 

L (s  q)  = q'(x)  + q(0)  6(x), 

where  5(x)  is  the  delta  function,  we  obtain  the  inverse  of  £ cosS  y as 

n n 

-1  2 

L [bQ  cos  eQy]  - q(x)  + U’W  + q (x)]>  (15a) 

and 

L ^ [b  cos3  y]  “ — ^ v — cos  (2QTTy) 
n n 2(nn)2c  h 

x [iu)q(x)  + U (q'(x)  + q(0)  5 (x))],  (15b) 

where  n ■ 1,  2,  ; the  derivation  of  (15b)  is  straightforward,  while 

that  of  (15a)  requires  some  consideration,  and  is  therefore  given  in 
Appendix. 
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C.  Final  Solution. 


With  the  aid  of  (13)  and  (15),  (11a)  becomes  as  follows: 


x mz 

m — f (mz)  B (x  - mz)  + — / g«  (x  - t)  B(t)  dt 
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(16) 
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where  fQ(x)  ■ U exp  (YqX)  Jq  [6q  (x  - m z ) , 
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In  the  limist  of  zero  acoustic  admittance,  one  can 

the  above  expression  will  be  reduced  to  a solution 

9 

oscillating  airfoil  immersed  in  supersonic  flow  . 


inanediately  show  that 
for  two-dimensional, 

1 


i 
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IV.  Discussion. 

Figure  2 show  the  comparison  of  the  unsteady  pressure  distribution 
between  treated  and  untreated  walls;  the  ordinate  Is  the  absolute  magnitude 
of  unsteady  pressure  and  the  abscissa  is  the  distance  measured  from  the 
bow  shock,  x - mz,  plotted  at  varous  values  of  z.  The  additional  parameters 
are  as  follows: 

airfoil  chord  = 1.0, 

(Mach  number)  = 1.3, 

h - 1.5, 

k - w/U  • M = 0.5, 

A = -1.140  + 0.5471  (for  treated  wall) 

=*  0 (for  untreated  wall). 

The  airfoil  is  executing  bending  motion  and  its  instantaneous  position 
is  given  by  z » e*Ut;  the  unsteady  pressure  in  the  figure  is  evaluated 
on  the  top  surface,  i.e.,  y * h. 

It  is  clearly  evident  that  the  provision  of  acoustic  liners 
significantly  attenuates  the  unsteady  pressure  fluctuation  and  tends 
to  suppress  the  occurrence  of  flutter. 
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APPENDIX 


Here  we  shall  derive  the  inversion  of  (15a):  we  note  that 
cos  eQy  becomes,  from  (12a) 


COS  80y  = COS  Ai^±Uslj2 

Taking  the  inverse  transform,  we  obtain 


-1 


L * (cos  60y]  = exp  (-i  ^)  L_1  ( cos  [y(-  %2  /^] 


ch 


= exp  (-i  !f)  L_1  [1  + ^ s + 0 (A2)  + _ j 


exp  (-i  ^)  (6(x)  + $'(x)  + 0 (A2)  + — ]. 


From  (14a)  and  with  the  aid  of  convolution,  one  obtains  (15a). 


18 


REFERENCES 


^Kurosaka,  M. , and  Edelfelt,  I.  H. , "Some  Recent  Development  in  Unsteady 
Aerodynamics  of  a Supersonic  Cascade,"  Proceedings  of  Symposium 
IUTAM  sur  L'Aeroelasticite  dans  les  Turbomachines,  Revue  Francaise 
de  Mecanique,  Numero  Special,  1976. 

2 

Platzer,  M.  F.,  "Unsteady  Flows  in  Turbomachines  - A Review  of  Current 
Developments,"  AGARD  Symposium  on  Unsteady  Aerodynamics,  September, 
1977. 

^McCune,  J.  E.,  "A  Three-Dimensional  Theory  of  Axial  Compressor  Blade 
Rows-Appllcatlon  in  Subsonic  and  Supersonic  Flows,"  Journal  of 
the  Aerospace  Sciences,  Vol.  25,  No.  9,  1958,  pp.  544-560. 

McCune , J.  E.,  "The  Transonic  Flow  Field  of  an  Axial  Compressor 
Blade  Row,"  Journal  of  the  Aerospace  Sciences,  Vol.  25, 

No.  10,  1958,  pp.  616,  626. 

^Okurounmu,  0.,  and  McCune,  J.  E. , "Three-Dimensional  Vortex  Theory 

of  Axial  Compressor  Blade  Rows  at  Subsonic  and  Transonic  Speeds," 
AIAA  Journal,  Vol.  8,  No.  7,  1970,  pp.  1275-1283. 

^Salaiui,  P. , "Instabilite  de  Flottement  dans  Une  Grille  Annulalre," 
Proceedings  of  Symposium  IUTAM  sur  L'Aeroelasticite  dans  les 
Turbomachines,  Revue  Francaise  de  Mecanique,  Numero  Special,  1976. 

^Erdelyl,  A. , et  al.,  "Tables  of  Integral  Transforms,"  Vol.  1,  1954, 
McGraw-Hill . 

g 

Kurosaka,  M. , "On  the  Unsteady  Supersonic  Cascade  with  a Subsonic 
Leading  Edge  - An  Exact  First  Order  Theory  - Part  1 and  2.", 

ASME  Transactions,  Series  A,  Journal  of  Engineering  for  Power, 

Vol.  96,  January,  1974,  pp.  13-31. 

9 

Miles,  J.  W. , "The  Potential  Theory  of  Unsteady  Supersonic  Flow," 
Cambridge  University  Press,  1959.  — 


Treated  wall 


h 


Untreated  wall 


APPENDIX  6 


" Towards  Simpler  Representation  of  Wave  Propagation 
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Towards  Simpler  Representation  of  Wave  Propagation 
Through  Non-Uniform  Media 


M.  Kurosaha 

The  University  of  Tennessee  Space  Institute 
Tullahomat  Tennessee  37388 

1.  Introduction. 

Pressing  problems  beset  with  unsteady  flow  phenomena  in  practical 
applications  appear  to  involve  almost  invariably  the  propagation 
of  disturbance  through  non-uniform  media  — such  as  flow  with  velocity 
and  temperature  gradients.  Notable  among  them  are:  the  flutter  problems 
in  turbomachinery  bladings  where  the  imparted  swirling  motion  induces 
substantial  radial  gradient  in  steady  base  flow;  the  aeroacoustic  noises 
propagating  through  aircraft  engine  ducts  whose  contoured  form  and  the 
boundary  layer  growth  over  its  inner  surface  result  in  pronounced 
variation  of  velocity;  and  high-intensity  sound  generated  in  heat 
exchangers,  gas-cooled  reactors  and  the  like  whose  complicated  internals 
inescapably  produce  considerable  spatial  distribution,  both  in  velocity 
and  temperature. 

Having  been  spurred  by  these  compelling  incentives,  the  last 
decade  has  witnessed  intensified  activity  in  the  study  of  acoustic 
propagation  through  non-uniform  flow  field  ; for  extensive  bibliography 
in  this  area,  the  reader  is  referred  to  the  recent  survey  articles  by 
Nayfeh,  Kaiser  and  Telionis  (1975)  and  Goldstein  (1976).  However, 
but  for  few  exceptions  (e.g.  Pridnore-Brown  (1958),  Shankar  (1971), 
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Goldstein  and  Rice  (1973)),  the  bulk  of  the  work  has  leaned  heavily  on 
entirely  computational  method.  Though,  needless  to  ;ay,  extremely  useful 
as  these  numerical  methods  are,  they  may  still  tend  to  fall  short  of 
achieving  the  desired  goal  of  revealing  the  explicit  functional  inter- 
relationship between  the  assorted  parameters  and  displaying  the  underlying 
physical  features.  Even  in  the  situations  which  are  amenable  to  analytical 
solutions,  they  appear  to  suffer  from  various  rigid  constraint  imposed  by 
the  particular  choice  of  velocity  profile  or  the  geometry  studied;  for 
example,  the  analytical  expressions  accrued  for  unbounded  media  may 
not  readily  be  applied  to  eigenvalue  problems  associated  with  ducted 
flow  and  often  they  are  altogether  unwieldy  for  physical  interpretations. 

It  therefore  appears  that  what  is  genuinely  needed  is  a simple,  albeit 
approximate,  technique  which  enables  one  to  cope  effectively  with  the  general 
wave  propagating  problems.  In  this  paper,  we  shall  describe  such  a technique 
which  yields  an  analytical  solution  of  disturbances  propagating  through  arbitrary 
velocity  profile  in  simple,  closed  form;  and  we  shall  present  it  from  a 
unified  viewpoint  where,  once  the  free-space  solution  is  derived,  a slight 
modification  of  the  results  will  promptly  furnish  both  a solution  and  physical 
interpretation  even  for  wave  guide  problems.  The  essential  idea  that  we  shall 
promote  is  not  unknown,  and  it  is  in  fact  a variant  of  'slowly  varying'  method, 
but  somehow  it  does  not  appear  to  have  been  pursued  in  the  way  it  will  be 
exploited  here;  with  the  assumption  of  small  velocity  (or  temperature)  gradient, 
the  technique  to  be  employed  here  is  strikingly  simple  yet  appears  to  be 
capable  of  exposing  physical  meaning  in  unmistakably  transparent  manner. 

In  the  problem  of  sound  propagation  through  non-uniform  media, 
the  original  governing  equations  are  invariably  reduced  into  a second- 
order  ordinary  differential  equation  with  varying  coefficients,  the 
latter  reflecting  the  non-uniformity  of  velocity  (and  temperature) 
profile;  the  equation  is  sometimes  called  Prldmore-Brown's  equation. 
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According  to  the  method  suggested  perhaps  first  by  Jeffreys  and  Jeffreys 
(1956,  p.  522),  we  recast  the  second-order  differential  equation,  via  ex- 
ponential transformation,  into  a first-order  equation.  The  equation  is  of 
Rlcatti's  type,  nonlinear  and  inhomogeneous;  the  transformation  is,  of 
course,  the  one  sometimes  used  as  the  first  step  in  deriving  the  WKB  for- 
mula. In  the  standard  WKB  method,  one  then  proceeds  to  solve  the  nonlin- 
ear first-order  equation  by  iterative  process.  Here,  instead,  we  seize 
upon  the  explicit  advantage  that  under  the  assumption  of  small  velocity 
gradient,  the  first  order  equation  turns  into  the  linear  one  with  constant 
coefficients;  a closed  form  solution  follows  immediately  from  this,  provid- 
ing a free-space  solution  for  sound  propagation  through  non-uniform  media. 

In  compact  form,  our  solution  (2.11)  gives  the  expression  for  wave  trans- 
mission through  any  arbitrary  shear  profile;  as  a special  case,  we  shall 
prove  that  the  solution,  when  applied  to  linearly  sheared  flow,  recovers 
the  result  corresponding  to  the  exact  solution  of  Goldstein  and  Rice  (ibid. ) , 
which  they  obtained  in  terms  of  Weber's  parabolic  cylinder  function.  The 
present  solution  also  includes,  for  another  specific  shear  profile,  a solu- 
tion of  Miles  (1957)  and  Ribner  (1957)  for  sound  transmission  through  sudden 
velocity  discontinuity.  These  shall  be  discussed  in  the  following  section,  2. 

Most  significantly,  from  the  wave-like  behaviour  of  the  solution 
explicitly  embodied  in  the  exponential  transformation  and  furthermore,  ow- 
ing to  the  simplification  arising  from  the  assumption  of  slightly  sheared 
flow,  the  present  solution  will  conspicuously  reveal  the  physical  features 
of  the  wave  propagation  in  non-uniform  media,  which  will  be  described  in 
section  3.  For  example,  we  shall  readily  recognize  that  for 
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i 

disturbances  propagating  slightly  but  arbitrarily  sheared  profile,  * , 

the  wave  front  advancing  past  a given  point  can  simply  be  constructed  | 

by  replacing  the  role  of  convective  velocity  in  uniform  flow  by  the  mean  ! 

i 

velocity  cumulatively  averaged  up  to  that  point  (as  might  be  antici- 
pated from  'slowly  varying'  concept);  this  appears  to  provide  the  * 

approximate  but  simpler  means  than  the  ones  where  the  wavefront  has 
to  be  successively  constructed  by  chasing  an  Instantaneous  front  surface 
after  surface.  After  these  descriptions  of  physical  implication,  we 
shall  compare  the  present  solution  with  the  ones  derived  from  WKB 
formula;  to  obtain  this,  we  transform  the  Pridmo re-Brown  equation 
into  one-dimensional  Schrodlger's  equation  and  apply  the  standard  WKB 
method;  we  shall  observe  that  the  present  solution,  when  rendered  into 
alternative  form  of  infinite  series  representation,  is  tantamount  to  providing 
all  the  higher-order  terms  lacking  in  the  solution  based  upon  the  standard 
WKB  method.  In  the  final  section  4,  we  shall  show  how  these  solutions  ob- 
tained for  unbounded  media  can  be  utilized,  with  slight  modification,  to 
waveguide  problems,  leading  to  closed-form  solution  (4.13).  In  addition, 
we  shall  observe  there  that,  with  the  aid  of  free-space  solution,  a certain 
interesting  phenomena  associated  with  some  class  of  velocity  distribution 
will  become  physically  interpretable. 

Our  undeviating  aim  in  this  paper  is  to  obtain  a solution  for  waves 
propagating  through  non-uniform  media  in  a form  as  uncomplicated  as  possible 
and  in  order  to  strive  for  this,  it  will  inevitably  be  subject  to  certain  formal 
limitations;  nevertheless  we  shall  find  that  the  present  approach  appears  to  be 
highly  effective,  both  in  providing  a compact  expression,  from  which  one  can 
easily  detect  various  trend,  and  in  exposing  physical  features. 
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2.  General  Solution  For  Unbounded  Medium. 


Consider  acoustic  waves  propagating  in  a two-dimensional  and 
flowing  medium;  the  base,  steady  flow  is  unidirectional,  with  velocity 
U in  the  positive  x direction  and  it  varies  in  the  transverse  y direction 
as  shown  in  Figure  1,  e.g.  U * U(y).  For  the  sake  of  simplicity,  the 
ambient  density  and  temperature  are  assumed  to  be  constant,  although,  if 
necessary,  the  present  technique  can  easily  be  employed  for  stratified 
medium  as  well;  the  flow  will  be  treated  as  inviscid.  In  this  section, 
we  consider  the  medium  to  be  unbounded  and  in  order  to  fix  our  ideas, 
we  shall  be  interested  in  acoustic  waves  travelling  from  y * -»and 
propagating  through  sheared  flow  in  the  positive  y direction;  at  y * -00, 
U is  assumed  to  be  uniform.  The  governing  equations  are  given  in  the 
linearized  form  as 


1 /9p'  . / x 3p*  v . _ ,3u'  , 3v’  x 

1 (aT~ + u <y>  + p (—  + — > “ 


3x 


3x 


3y 


3u* 

3t 


+ U 


(y) 


v.  Mil 

dy 


1 1e1 

p 3x  ’ 


> 


(2.1) 


3v' 

3t 


+ u <y>  Ir 


I i£l 

p 3y  ’ 


where  p and  c are  unperturbed  density  and  the  speed  of  sound,  respectively, 
which  are  constant;  u'  and  v'  are  the  fluctuating  part  of  the  velocity 
component  in  the  x and  y direction,  respectively;  p'  is  the  perturbed 
pressure  and  t is  time.  Eliminating  u'  and  v'  from  (2.1),  one  obtains 
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3x  3y  c c 3t 


2 dU(y)  3_£l. 

dy  3x3y 


(2.2) 


where  M (y)  « U (y)/c  is  the  usual  Mach  number.  If  one  considers  a 
sinusoidal  wave  with  frequency  ui  and  travelling  with  wave  constant  a 
in  the  x direction,  the  fluctuating  pressure  takes  the  following  form: 


P (y)  e 


i(otx  - wt) 


(2.3) 


Parenthetically,  this  embodies  a klnematical  condition  that  the  wave 
number  in  the  x direction,  ot/27r  , remains  always  constant  even  though  the 
waves  suffer  refraction  due  to  shear  flow  in  the  y direction.  Into 
(2.2),  we  substitute  (2.3)  and  obtain 


P"  " p1  + t(aM(y)  " k)2-a2l  p - 0.  (2.4) 

where  k » tu/c  and  the  primes  denote  the  differentiation  with  respect  to 
y;  the  above  ordinary  equation  involving  variable  coefficients  is 
sometimes  called  Pridmore-Brown's  equation  (Pridmore-Brown,  1958). 

For  the  extremely  simple  situation  of  M (y)  ■ constant  * M^,  say» 
the  solution  is,  of  course,  given  by 

P * exp  [+iy  ✓(oMq  - k)2  - a2]  . 


(2.5) 
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i 

i 

i 
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Now  consider  slightly  sheared  flow  with  arbitrary  profile 


M (y)  - MQ  + e g (y) , 


(2.6) 


where  e is  a small  parameter  and  the  shear  distribution  g(y)  is  any 
given  function  of  y;  we  take  Mq  to  be  the  Mach  number  at  y a - “ 
(hence  g (-  °°)  * 0).  Motivated  by  the  form  of  (2.5)  for  uniform  flow, 
we  write  in  the  present  case 


p = exp  [+  i y 


/(aM„ 


. ,2  2 

- k)  - a + e 


F (O  d5) 


(2.7) 


Substitution  of  (2.7)  into  (2.4)  at  once  yields  the  following  equation, 
to  the  order  of  e,  for  F: 

^ ± 2 i y F + [2  (oM0  - k)  a g (y)  + ] = 0,  (2.8a) 

where 

y - /(aM0  - k)2  - a2  . (2.8b) 

It  is  crucial  to  recognize  here  that  the  above  first  order  equation  is 
linear  and  the  coefficients  are  constant;  were  the  perturbation  scheme 
not  explicitly  introduced  in  (2.7),  the  exponential  transformation  would 
have  led  to  the  nonlinear  Ricatti's  equation.  In  the  above  (2.8a),  we 
also  call  our  attention  to  the  fact  that  shear  flow  profile  g (y),  which 
is  a function  of  y,  now  appears  only  within  the  bracket  or  as  non-homogeneous 
term.  With  regard  to  the  solution  of  (2.8a),  the  homogeneous  solutions 
is  obviously  equal  to  exp  (+  2 i y y) . But  since  this  does  not 
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contain  any  effect  of  shear,  the  solution  Is  physically  meaningless 
and,  therefore,  is  discarded;  only  the  particular  solution  of  (2.8a) 
is  relevant.  In  order  to  obtain  the  latter,  we  take  the  Fourier  transform 
of  (2.8a)  and  obtain 


iaF  + 2 iyF  = - 2 (aHQ  - k)  a g + ^p'_gk  . (2.9) 


where  F and  g are  Fourier  transform  of  F and  g,  respectively,  defined  by 


1 

/2tt 


F eiCTy  dy,  I 


/2rf 


iOy 

g e dy. 


Solving  for  F and  inverting,  we  obtain 
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Thus  we  have  already  attained  the  desired  goal  of  obtaining  an  expression 
for  waves  propagating  through  arbitrarily  sheared  flow  in  the  positive  y 
direction.  (For  waves  propagation  in  the  negative  y direction,  one  only 
has  to  change  the  lower  limits  of  the  integrals  from  - 00  to  +<*».)  If 
the  flow  is  uniform  up  to  certain  height,  say,  y ■ y^,  i.e., 

g (y)  - 0 for  y < yQ  , 

(2.11a)  becomes 


Although  the  case  of  y ■ 0 or  zero  incidence  angle  is  physically  trivial 
for  the  present  free-space  transmission  problem,  we  record  here,  for 
subsequent  reference  to  be  used  in  ducted  wave  analysis,  that  at  y • 0 
(2.11b)  becomes 


A 


t 

l 
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For  given  shear  flow,  the  propagating  wave  at  any  point  can  be  determined  from  i 

either  of  the  above,  once  the  integration  involved  is  carried  out.  A 
note  of  caution  is,  however,  in  order  here.  When  the  integration  over  g 
is  carried  out  for  specified  shear  profile,  sometimes  the  terms  in  the  form 
of  exp  (+  2iyy)  appear.  Since  they  are  of  the  same  form  as  the  afore- 
mentioned homogeneous  solution  of  (2.8a),  these  terms  are  parasitic  in  their  na- 
ture and  therefore  they  should  always  be  discarded  for  the  reason  already  stated. 

In  what  is  to  follow,  we  shall  work  out  examples  for  three  shear 
distribution:  (a)  shear  profile  with  exponentially  decaying  transition 
(b)  linear  shear  and  (c)  sudden  velocity  discontinuity  (see  Figure  2). 

In  (b)  and  (c),  we  shall  observe  that  the  present  results  agree 
satisfactorily  with  the  ones  based  on  the  available  solution. 


Example  1.  Shear  profile  with  exponential  transition: 
g(y)  ■ 1 - e for  y >.  0 

and  g(x)  - 0,  for  y <_  0 . 

Upon  substitution  of  g(y)  into  (2.11b)  where  we  take  y^  « 0 and  discarding 
the  terms  of  exp  (+  2iyy)  after  carrying  out  the  integral,  one  obtains, 
when  y t 0 


p’  - A e1(ctx  " wt)  e±  iYy 


This  represents  the  transmitted  wave  in  the  region  above  y **  0 when  the 


incident  wave  for  y < 0 is  prescribed  to  be  p'  ■ e1^01*  ~ Ut^e^  We 

particularly  emphasize  here  that  the  above  expression  provides  even  the 
correct  amount  of  the  amplitude  of  the  transmitted  wave  corresponding 
to  the  amplitude  of  Incident  waves  specified  to  be  unity.  (If  other 
than  unity,  the  amplitude  of  the  transmitted  wave  should,  of  course,  be 
adjusted  proportionately) . This  is  the  reason  why  the  constant  A is 
deliberately  retained  in  the  above  and  there  is  no  need  to  multiply 
it  by  some  other  additional  factor.  One  can  directly  confirm  this  point 
by  applying,  at  the  interface  of  y ■ 0,  the  continuity  of  pressure  and 
the  particle  displacement  condition  between  the  expression  of  (2.12), 
valid  for  y _>0,  and  the  expression  for  y <_  0;  the  latter  is  given  by  the 
following  combination  of  plane  incident  and  reflecting  waves 

p*  - ei(ax  - «t)  e±  lYy  + R ei(ax  " wt)eT  iyy,  for  y < 0 
where  R is  the  amplitude  of  the  reflected  wave  to  be  determined  in  the  process. 
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Example  2.  Linear  Shear  : g(y)  ■ y for  y >_  0 and  g(y)  - 0 

for  y <.  0. 

In  thla  case,  where  at  y ■ 0 Che  base  flow  changes  from  uniform 
distribution  to  linearly  varying  one,  aquation  (2.11b)  becomes  as 
follows: 


p»  . eK«x  - wt)^  lyy 


(2.13) 


x exp  < + 1 


ea(aH0  - k)  2 a(u*Q  " 


2y 


y + e l- 


2 Y 2 


aM0  - k 


J y > • 


In  obtaining  this,  the  term  In  the  form  of  exp  (+  2 lyy)  , which  has  arisen 
in  the  course  of  Integration,  Is  again  discarded.  As  before,  (2.13)  re- 
presents the  wave  transmitted  in  the  region  above  y ■ 0 when  the  incident 
wave  impinging  on  y « 0 is  specified  to  be  p’  « e1^0®  “ . 

For  linearly  sheared  profile,  Goldstein  and  Rice  (1973)  obtained 
the  exact  solution  of  propagating  waves  in  terms  of  the  parabolic  cylinder 
function.  In  the  present  notation,  their  solution  can  be  expressed  as 


i 

I 

! 

l 


] 

, i 

' \ 

i 

i 

: i 


j 


i 


i 

I 


p' 


ei(dx  - wt) 


[U'(b,  +()+b(U(b,  +€)], 


(2.14) 


> 


where  U is  the  parabolic  cylinder  function  (e.g.  Abramowltz  and  Stegun 
(1970),  p.  686)  defined  by  the  following  integral  representation: 


When  the  gradient  o£  shear,  , Is  small,  both  b and  £ defined 

above  tend  to  become  large.  In  such  a case,  one  may  prove,  as  shown  In 

Appendix  A,  that,  with  the  aid  of  the  Darwin's  formula  (1949)  for  the 

parabolic  cylinder  function,  (2.14)  Is  reduced  to  the  present  result, 

equation  (2.13).  (In  lieu  of  Darwin's  formula.  It  might  appear  possible 

to  start  with  the  Integral  representation  of  (2.15a)  and  apply  the  Laplace's 

method  for  large  values  of  b and  £•  It  turns  out  that  this  approach 

calls  for  the  usage  of  the  generalized  Laplace's  method  and  so-called 

Faxen  integral,  both  of  which  are  discussed  In  Olver  (1974,  p.  331-332). 

The  method  is,  however,  apparently  somewhat  too  crude  to  serve  the  present 

2 

need,  for,  although  this  accurately  reproduces  y term  in  (2.13),  it 
fails  to  yield  the  correct  coefficients  of  y.) 


Example  3.  Sudden  Velocity  Discontinuity: 
g (y)  ■ AM  for  y > 0 
and  g (y)  ■ 0 for  y < 0. 

When  the  steady  velocity  increases  suddenly  by  AM  at  y ■ 0, 

(2.11b)  yields 

p'  - A ei(ax  ' wt)  e - iYy  exp  [+  iect  (cMQ  - k)  , (2.16a) 
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where 

A “ exp 


CAM 


a[2y2  - (oMQ  - k)2] 
2y2  (CM0  - k) 


(2.16b) 


Once  again,  this  Is  the  wave  transmitted  in  the  region  above  the 
interface  y ■ 0,  upon  which  the  incident  wave  impinges. 

The  exact  solution  for  sound  transmission  through  a flow  velocity 
discontinuity  was  obtained  by  Miles  (1957)  and,  independently,  by 
Ribner  (1957).  According  to  them,  the  transmitted  wave  is  given  by 


p'  - T ei(ax  " ut>  ei0ty  tan  ^2  ’ 


(2.17a) 


where 


2 sin  (2  *x) 

T " sin  (2  4>1)  + sin  (2  ^)  * 

In  the  above,  4^,  and  4>2  are  incident  and  refracted  wave  angle, 
respectively,  measured  from  the  positive  x direction.  is  given  by 


(2.17b) 


1 

a 


tan  4^  ; 


(2.17c) 


and  4>2  is  related  to  4>^  by  the  following  law  of  refraction  analogous 
to  the  Snell's  law: 


1 

cos  4^ 


1 

cos 


+ e AM. 


(2.17d) 


m 


l 


When  the  amount  of  velocity  jump  is  small,  one  can  readily  prove 
that  (2.17)  becomes  Identical  to  (2.16),  the  details  of  which  are  shown 
in  Appendix  B.  Additional  numerical  comparison  between  the  Miles-Ribner's 
exact  solution  and  the  present  approximate  solution  is  shown  in  Table  1, 
where  the  transmission  coefficient,  A of  (2.16b)  for  the  present  solution, 
are  compared  with  T of  (2.17b)  for  various  incident  angles,  they 

are  given  both  for  positive  and  negative  velocity  discontinuity.  It  is 
clear  that  the  present  formula  provides  sufficiently  accurate  values 
for  a wide  range  of  incident  angles,  except  for  the  vicinity  of  either 
critical  incidence  (corresponding  to  total  reflection  for  positive 
velocity  jump)  or  the  zero  incidence  in  the  case  of  negative  velocity  jump 
(The  local  discrepancy  near  these  two  incidence  angles  Is  not  quite 
unexpected,  since  both  correspond  to  the  turning  point  of  (2.4),  about 
which  we  shall  discuss  more  in  section  3;  for  AM  > 0,  the  turning  point 
occurs  in  the  region  where  y > 0 and  for  AM<  0,  in  y < 0.) 


3.  Alternative  Representation  and  Physical  Interpretation. 


We  now  return  to  the  general  solution,  (2.11),  representing  waves 
propagating  through  arbitrary  shear  profile;  though  convenient  for 
obtaining  compact  expression  of  such  waves,  it  is  not  appropriately 
suitable  for  exposing  physical  features.  Therefore,  in  order  to  pave  the  way 
for  extracting  physical  interpretation,  we  recast  it  in  alternative  form. 
Taking  the  case  of  shear  distribution  where  the  flow  remains  uniform  up 
to  y ■ 0 and  whose  solution  is  given  by  (2.11b),  we  apply  the  integration 
by  parts  to  its  double  integral  appearing  in  the  argument  of  exponential 
and  obtain 

p,  _ ei(ax  - wt)^  iyy 


x exp  / +ie 


oKoMq  - k) 


g (n)  dn 


+ e 2a«*0-k>J 


a-^-5-  g(y)  -*-iTg,,(y)  +J-rg(4)(y)  — ] 


± U-^g'  (y) g"'(y)+ — g(5)(y)  — ]>}, 

(2  y)J  (2y  y (2y/ 

(3.1)  J 


where,  as  before,  the  terms  in  the  form  of  exp  (+  2 i yy)  originated  from 
the  lower  limits  of  the  integral  are  discarded.  The  above  expression,  (3.1), 
is  a series  expansion  of  (2.11b)  in  terms  of  y.  Though,  in  general,  it  is 
less  expedient  for  computational  purpose,  on  the  other  hand  it  readily 
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surrenders  to  the  physical  Interpretation,  whose  term-by-term 
description  will  be  given  below;  pursuant  to  this,  we  shall  compare 
(3.1)  with  the  other  series  solution,  which  we  derive  by  applying  the 
standard  WKB  method. 


3.1  Term  Associated  with  j g (V  dj. 


Before  we  consider  the  first  term  In  the  exponential 


+ 1 e 


a(oM0  - k) 
Y ~ 


J g(n)  <*n  , 

o 


(3.2) 


we  note  that  for  uniform  flow  ( e - 0) , the  surfaces  of  crests  and  troughs  or 
of  constant  phase  are  given  by 


ox  + yy  - constant. 


(3.3a) 


and  the  surfaces  move  as  time  changes.  For  the  present  case  of 
e l*  0,  the  term  (3.2),  which  is  the  leading  correction  term  for  large 
values  of  y or  at  higher  frequency,  modifies  the  above  Into  the  following: 


, a(aM0  - k) 

ax  + yy  + e 


g (n)  d n * constant 


(3.3b) 


The  Integral  form  indicates  that  phase  change  is  dependent  on  the 
entire  flow  field  over  which  the  wave  has  traversed;  hence,  the  term 
embodies  cumulative,  memory  content  and  might  be  called  "phase  memory".  (Though 
more  popularly  found  in  electromagnetic  propagation,  this  phase  memory  recently 
surfaced  in  the  context  of  unsteady  flow  field  radiated  by  an  oscillating 
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body  (Kurosaka,  1977).)  In  order  to  appreciate  the  meaning  contained 
in  (3.3b)  more  fully,  we  replace  the  uniform  velocity  appearing 
In  the  following  definition  of  y» 


.A 


Y - ✓ (oiM0  - k)2  - a2  , 


cy  the  mean  velocity  averaged  up  to  the  vertical  position,  y,  i.e.. 


m-  ^ M(n)dn  - mq  + ^ | g(n)  dn. 


- 2 2 1/2 

Then,  if  we  define  y to  be  [(aM  - k)  - a ] , this  becomes, 

to  the  present  order,  as 


7-  -Y+  £a(aMo  ~ k)  1 


g (n)  dn  , 


and  (3.3b)  becomes 


ax  + yy  - constant 


(3.4) 


Comparison  of  the  above  with  the  expression  for  uniform  flow,  (3.3a), 
immediately  reveals  the  following:  the  surfaces  of  constant  phase  for 
waves  propagating  through  shear  flow  passing  a vertical  point,  y,  may  be 
given  (as  a good  approximation  for  shortwave  length  in  the  y direction  or  at 
high  frequency)  by  replacing  the  role  of  uniform  velocity  by  the  mean 
one  cumulatively  averaged  up  to  that  point,  y.  From  the  intimate 
relationship  between  the  behaviour  of  the  waves  at  high  frequency  and  the 
wave  front  (e.g.  Whit ham,  1974,  p.  236),  one  is  naturally  led  to 
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expect  that  the  same  replacement  might  take  place  for  wave  fronts. 

In  order  to  obtain  the  direct  confirmation  of  this,  we  regard  the 
solution,  correct  to  the  present  order  and  expressed  in  terms  of  the 
above  y,  as  the  Fourier  transform,  F,  with  respect  to  a and  to:  that  is 


The  expression  for  the  wavefront  will  accrue  naturally  when  one 
takes  the  inverse  transform  so  as  to  derive  the  transient  response, 
(Miles,  1959,  p.  53;  Kurosaka,  1977).  Taking  the  inverse  transform 
of  (3.5a)  with  regard  to  a and  to  , 


F - 


00  oo 


1 gi(ax  - tot)e  + iyy 
/2n 


(3.6a) 


where  U is  the  averaged  shear  velocity  defined  by 
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and  H is  the  step  function;  the  derivation  of  this  is  relegated  to 
Appendix  C.  From  (3.6b),  it  is  clear  that  the  wave  front  is  given  by 

2 2 2 — 2 

c t - y + (x  - Ut)  . (3.8) 

2 2 2 2 

Comparison  with  the  one  for  uniform  flow,  c t = y + (x  - U^t)  , 
readily  shows  that  the  wave  front  propagating  through  non-uniform  flow 
can  easily  be  obtained  by  replacing  the  role  of  uniform  velocity  by  the 
cumulatively  averaged  one,  as  might  be  anticipated  from  'slowly-varying' 
notion.  We  emphasize  the  point  that  the  above  approximate 
but  analytical  expression  offers  a simpler  way  to  construct  wavefronts 
than  the  general  method  where  it  has  to  be  successively  constructed  from  the 
initial  surface  by  the  vectorial  addition  of  acoustic  velocity,  which 
is  in  the  direction  of  normal  to  the  surface,  and  fluid  velocity. 

3. 2 Term  Associated  with  g(y). 

With  respect  to  the  term  associated  with  g(y)  in  the  argument  of 
the  exponential  in  (3.1), 


£ t'arVv  - 2«  <<*0  ' k)]  —2  8(y)’ 

0 u (2y) 

we  observe  that  this  is  wholly  real  and  consequently  represents  the 
primary  change  of  wave  amplitude  due  to  the  presence  of  shear.  It 
may  be  of  interest  to  note  its  following  implication:  when  we 
express  the  amplitude  corresponding  to  the  above  by  A,  i.e. , 


1 


(3.9) 


- 2a  (aM  - k)]  — 

0 (2y) 


g(y)] 
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then 

^ A • g'(y)  . (3.10) 

Therefore,  to  the  first  order  of  short  wavelength  in  y direction  or  at 
high  frequency,  the  rate  of  the  change  of  the  amplitude  is  proportional 
to  the  product  of  velocity  gradient  at  that  position  and  the  amplitude 
itself.  (It  satisfies  the  equation  analogous  to  the  one  for  radioactive 
decay . ) 

3.3  Terms  Associated  with  derivatives  of  g(y) 

The  remainder  of  the  terms  in  the  argument  of  exponential  in  (3.1), 

- 2“ 

x ([-  — • g"  (y)  + — ^ g(4)(y)  — ] 

(2Yr  (2y) 


± i [ — H g'  (y) g'”  (y)  + g(5)y ]), 

(2y)^  (2y)  (2Y)' 

(3.11) 

obviously  provides  the  higher-order  correction  to  both  amplitude 
and  phase.  If  we  consider  a shear  profile  where  at  some  point  y “ y^, 
all  of  its  derivatives  vanish  (see  Figure  1 again),  i.e.,  g'(y^)  ■ 
g"(y1)  * — - * 0,  the  wave  transmitting  through  y does  not  perceive 

any  local  change  in  steady  velocity;  and  therefore,  both  the  amplitude 

i 

| and  phase  should  remain  unaltered  while  traversing  past  y^.  This 

| . fact  is  indeed  accounted  for  in  the  present  solution,  since  in  such 


t 
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an  instance,  (3.11)  vanishes  completely  while  (3.2)  and  (3.9)  retain 
the  same  value  just  before  and  just  after  leaving  y^. 


3.4  Co 


n-son  with  the  Series  Solution  Obtainable 


bu  the  Standard  WKB  Method. 


We  now  compare  the  alternative  series  representation,  (3.1), 
with  the  other  series  solution  derivable  from  the  standard  WKB  formula. 
For  this,  we  transform  the  Pridmore-Brown's  equation  (2.4),  into  the 
Schrodinger's  equation  for  one-dimensional  quantum-mechanical  motion 
where  the  standard  form  of  WKB  solution  is  readily  available.  First, 
we  remove,  in  (2.4),  the  term  in  first  derivative  by  the  following 
standard  transformation: 


p - zY,  (3.12a) 

where 

z ■ aM(y)  - k . (3.12b) 

This,  together  with  (2.6),  changes  (2.4)  into 


Y"  + [y2  - e U(y)]  Y - 0 , (3.13a) 


where 


i 


i 


(3.13b) 
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This  is  the  Schrodinger 's  equation  where  U (y)  corresponds  to  potential 
and  so  far  no  approximation  has  been  made.  The  standard  WKB  formula 
(e.g.  Morse  and  Feshbach,  1953,  p.  1092  et  seq.)  is  applied  here  and 
the  approximate  solution  of  (3.13a)  is  immediately  found  to  be 


(3.15) 


To  this,  we  restore  the  factor  z of  (3.12b),  which  we  recognize  to  be  equivalent, 
to  the  order  of  the  present  approximation,,  to 
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Multiplying  the  result  by  e^0*  we  obtain  p'  as 

; 

p.  . ei(ox  - “Ogl  ±yy 

i I +iea(oM0  - k)  f 

f X exp  / / g (n)  dn 

; L o 

! 2 

+ e ■ 2a  (OMo  ■ k)1  — « (y> 

“"o  k 0 (2Y)2 


I 

i 


We  take  a special  note  that  the  argument  of  the  above  exponential  function 
consists  of  finite  terms,  in  contrast  to  the  infinite  series  appeared  in 

(3.1) .  Further  scrutiny  between  the  above  WKB  solution,  (3.16),  and  our 

y 

series  representation,  (3.1),  reveals  that  the  terms  associated  with  l g (n)  dn 

0 

and  g (y)  agree;  the  other  terms  associated  with  g' (y)  and  g"(y)  are  not 
quite  the  same;  g' ' ' (y)  and  other  higher  derivatives  are  completely  missing 
in  the  WKB  solution.  One  can  show  that  whereas  the  present  series  solution 

(3.1) ,  or  more  precisely  its  y-dependent  portion  corresponding  to  p(y) 

of  (2.3),  satisfies  the  original  differential  equation  (2.4)  to  the  order 
> of  e,  the  WKB  solution  (3.16)  falls,  not  surprisingly,  to  do  so  even 

to  0 (c).  Indeed,  it  is  straightforward  to  prove  that  if  the  iteration 
procedure  described  by  Morse  and  Feshbach  (ibid.),  which  leads  to  the  standard 
form  of  WKB  solution  of  (3.14a),  had  been  repeatedly  continued,  the  resulting 
higher  order  terms  would,  to  the  order  of  £,  contribute  to  recover  the  entire 
terms  missing  in  (3.16).  Consequently,  the  present  infinite  series  representation. 


+ e^k  It-S*'  <y>  --4 

aMo  K (2yr  (2Y) 


j 

(3.16)  i 
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(3.1),  does  offer  the  terms  of  higher  order  which  would  be  lacking  in 
the  WKB  solution;  the  primary  expression  obtained  in  the  preceding  section 
2,  (2.11),  amounts  to  the  summation  of  these  terms  into  a concise  form 
convenient  for  arbitrary  shear  profile.  Observe  also  that  without  the 
assumption  of  e <<1,  the  solution  (3.14a)  directly  derived  from  the 
standard  WKB  formula  would  not  be  amenable  to  the  physical  Interpretation 
described  previously. 

Before  concluding  this  section,  we  wish  to  dissipate  a potential 
source  of  misconception  about  summation  formula,  (2.11):  on  the  face 
of  it, this  solution  might  appear  to  be  valid  even  in  the  neighborhood  of 
Y ■ 0,  which  is  a singular  point  in  the  series  representation  of  (3.1) 
or  WKB  solution  (3.16).  (The  turning  point  of  (3.14b)  is  q ■ / Y2  - e U(y)  - 0 

but  the  expansion  for  e « 1 has  transferred  the  singular  point  to  Y * 0, 
instead).  However,  except  for  the  case  where  y happens  to  be  exactly 
equal  to  zero,  this  is  not  so;  for  once  the  integration  for  (2.11)  is  carried 
out,  the  terms  like  1/y  reappears.  For  example,  one  can  immediately  recognize 
this  from  (2.12),  which  — although  of  compact  form  than  the  one  that  would 
be  obtained  from  the  infinite  series  representation  corresponding  to  (3.1)  — 
is  still  singular;  in  example  3 of  section  2,  we  also  have  observed  the 
similar  behaviour.  To  avoid  undue  complication,  we  defer  any  attempt  to  improve 
this  limitation  obviously  associated  with  the  regular  perturbation  scheme  to 

elsewhere,  but  in  the  next  section  we  shall  observe  that  even  in  the  present 
form,  the  expression  for  Y ■ 0,  (2.11c),  will  be  found  to  be  valid  and  of 
physical  significance  in  ducted  acoustic  waves. 
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4.  Application  to  Ducted  Shear  Flow  Problem. 

Let  us  turn  our  attention  toward  the  problem  of  acoustic  wave 
propagation  through  shear  flow  in  ducts;  here  we  shall  find  that  the  results 
of  the  previous  sections  obtained  for  unbounded  media  can  be  applied  to 
channel  flow  with  minor  modification  — in  particular,  without  any  change  in 
the  case  of  the  fundamental  mode. 

The  coordinate  system  is  shown  in  Figure  3 where  the  duct  walls  are 
placed  at  y * + £ . Assuming  that  the  walls  are  rigid  or  hard,  the  boundary 
conditions  on  the  walls  are  given  by 


3y 


at  y - + £ 


(4.1) 


In  order  to  ascribe  definite  meaning,  the  steady,  base  velocity,  M^,  of  (2.6),  is 
now  taken,  without  loss  of  generality,  to  be  the  one  averaged  across  the  duct,  i.e.. 


M0  - ii  j M(y)  dy* 

-i 


(4.2) 


and  accordingly 
£ 


g(y)  dy  - 0 


(4.3) 


-£ 


If  the  flow  inside  the  duct  is  uniform,  it  is  of  course  elementary 
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to  show  that  the  wave  constant  in  the  y direction,  y of  (2.8b),  takes 
the  following  eigenvalues  determined  from  the  wall  boundary  condition: 
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nir 


Y m ~s~  where  n ■ 0,  + 1,  + 2, 

n l — — 

and  correspondingly  Inverting  (2. 8b), a is  given  by 


(4.4) 


MQk  + /k2  - 


°n 


(1  - H02)  (f) 


Mq  -1 


for  n ■ 0,  + 1,  + 2,  , (4.5) 


where  we  assign  suffix  n to  a.  In  the  above,  the  positive  sign  corresponds 
to  the  waves  propagating  in  the  upstream  direction  and  the  negative  sign 
to  the  one  in  the  downstream  direction. 

In  the  particular  case  of  the  fundamental  mode  of  plane  wave 
where  n ■ 0,  y for  uniform  flow  becomes  zero  from  (4.4).  If  we  return  to 
the  present  subject  of  sheared  flow,  then  for  y ■ 0 the  solution  has 
previously  been  obtained  as  (2.11c)  and,  when  one  takes  the  lower  limit  of 
the  integrals  to  be  equal  to  -l , we  obtain  the  following: 


P-  - e1(ctx  " ut)exp[-2ea  (<*„  - k)  / dn  / g(0  d£)  . (4.6) 

-l  -l 

Upon  differentiation  with  respect  to  y,  it  is  immediately  apparent  from 
(4.3)  that  this  does  satisfy  the  boundary  condition,  (4.1);  consequently 
this  is  in  fact  the  lowest  mode  for  waves  propagating  in  duct;  the 
underlying  physical  reason  why  in  this  case  the  free-space  solution  also 
becomes  a ducted  wave  will  be  given  subsequently.  For  now,  from  the  expression  of 
a corresponding  to  n - 0 in  (4.5),  one  obtains 


b ‘ k)  j dn  Js 


_-iwt  x \ — . 2k' 

e exp  (+ik  - — — — r-)  exp  [+  e — 

M0  + 1 (1 


y/A  n 

SW-  f 


8(Od(f)U 


-l 


A 
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where  suffix  0 attached  to  p'  denotes  the  fundamental  mode  and  upper  signs 
inside  the  exponential  correspond  to  the  upstream  propagation  and  the  lower 
signs  to  the  downstream  propagation.  During  the  course  of  the  present  work, 
it  was  communicated  by  Dr.  Savkar  that  the  essentially  same  expression  was 
previously  obtained  by  him  in  hitherto  unpublished  memo  (1972),  which  is 
now  being  prepared  for  publication,  by  using  an  entirely  different  method 
where  a technique  akin  to  PLK  method  was  employed  along  with  certain  ortho- 
gonality condition  for  various  duct  modes;  on  the  other  hand,  as  we  have  just 
seen,  the  present  derivation  is  directly  obtained  from  the  free-space  sol- 
ution and  not  dependent  upon  the  latter  requirement. 

Of  interest  here  is  the  comparison  of  the  present  expression,  (4.7), 
with  Shankar's  solution  (1971)  for  initial  value  problem.  He  posed  and 
solved  a problem  where  a harmonic  plane  wave  disturbance  is  suddenly  switched 
on  at  x * 0 in  an  initially  quiescent  medium  with  slightly  sheared  flow.  For 
the  Initial  condition 

p'(0,y,t)  * A e iWtH(t) 

where  A is  constant  (plane  wave)  and  H is  the  step  function,  his  solution 
for  large  time  is  given,  according  to  the  present  definition  of  Mq,  to  be 
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x [ exp 


[ 


M.  2 2_2  - , 

* 2 ("  T + 1I7  - a ' Mo  >]  >' 

1 - c c b 


exp  ^ 1 7uVm0)  ? 1 + 0 (e2)  ’ 


(4.8) 


where  b - 2%  , y*  * y + £ and  aR  is  the  Fourier  coefficient  of  shear 
profile  g(y)  * (y*)  defined  by 


M(1)  (?)  cos  (^)  d C . 


We  combine  the  first  term  of  (4.8)  and  the  terms  of  the  same  exponential 
form  appearing  Inside  of  the  series,  i.e.. 


exp  ( - iu  (t  - 


c(l  + Hq)' 


E4  an  ^ , 2 . ,niTY*v  -iwt  iux 

-5 5-  -5-  b k cos  ( — f—)  e exp  r 1 . 

if2 (1  + Mn) 2 n2  b l c(l  + M0)  j 


n * 1 


(4.9) 


Upon  differentiating  twice  with  respect  to  y,  it  is  straightforward  to  show  that, 
from  the  Fourier  series  formula,  the  sum  is  Indeed  equal  to  the  present  result  of 
(4.7)  for  downstream  propagation,  correct  to  the  order  of  e.  ( The  other 
terms  in  the  series  of  (4.8)  obviously  correspond  to  higher  modes  associated 


with  the  particular  choice  of  the  initial  condition  which  is  assumed 
to  be  a plane  wave  in  Shankar's  paper.) 

The  present  expression  for  the  fundamental  modes,  (4.7),  reveals 
explicitly  various  features  of  the  shear  flow  effect.  For  example, 
for  a parabolic  shear  flow  satisfying  the  requirement  of  (4.3),  i.e. 
g(y)  = 1/3  - (y/f,)2 


p'o  * e"iwt  exp  (± 1 kM^7i)  exp  ( 1 e 6?H~v2  l(i)  ’ 21 


(4.10) 


where  the  upper  and  lower  signs  correspond  to  upstream  and  downstream 


propagation,  respectively.  Not  only  the  role  of  amount  of  shear  (e) , frequency 

(k)  and  mean  Mach  number  is  transparently  obvious,  but  also,  when  we  note 
2 

that  [(y/£)  - 2]  always  remains  negative  across  ducts,  it  is  plainly  visible 

that  shear  layers  refract  the  fundamental  mode  wave  toward  the  wall  for 

downatream  propagation  and  away  from  the  wall  for  upstream  propagation. 

The  effect  has  been  known  from  the  various  results  based  upon  numerically 

computational  scheme  (see  Nayfeh  et  al.  for  full  references.);  the  magnified 

effect  of  refraction  at  higher  frequency  and  greater  shear  has  likewise 

been  recognized.  We  feel,  however,  that  the  present  analytical  expression 

displays  these  effects,  perhaps  more  conspicuously  and  compactly.  The 

numerical  comparison  of  (4.7)  with  the  results  obtained  by  wholly  computational 

scheme  (e.g.  Mungur  and  Gladwell  (1969))  for  prescribed  shear  profile  such 

as  linear  distribution  shows  the  satisfactory  agreement, even  when  the  values 
2 2 

of  ek  l appearing  inside  the  exponential  of  (4.7)  become  moderate; 
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n 


therefore  in  this  paper,  in  the  interest  of  retaining  simplicity,  we 

shall  not  attempt  to  improve  the  , formal  limitation  at  large  values 
2 2 

of  ek  l through  the  obvious  application  of  singular  perturbation  procedure. 

For  higher  modes  in  ducted  flow,  the  free-space  solution  obtained 
previously,  requires  some  modification.  Since,  for  uniform  flow,  the 
wave  constants  in  the  x direction,  otQ  of  (4.5),  are  dependent  upon  the 
uniform  velocity,  Mq,  we  expect  that  for  sheared  flow,  it  is  dependent 
on  shear;  accordingly,  we  express 

a * - a + ex  (4.11) 

n n n 

where  is  the  correction  factor  due  to  shear.  In  general,  ^ has  to  be 
evaluated  from  the  boundary  condition.  However,  as  in  the  previous  case  of 
fundamental  mode  of  n ■ 0,  which  in  the  present  notation  corresponds  to 
Xq  • 0,  one  can  show  that  even  for  higher  mode  vanishes  for  a certain  shear 
profile;  that  is,  for  odd  component  of  velocity  distribution  in  the  y 
direction,  )^  ■ 0 or  the  wave  number  in  the  x direction  remains  unaffected  by 
the  presence  of  shear.  This  fact  appears  to  be  first  recognized  by  Savkar 
in  the  aforementioned  unpublished  memo  based  upon  the  duct-mode  analysis; 
in  what  is  to  follow,  we  shall  instead  show  this  from  physical  consideration 
based  upon  free-space  solution. 

The  sound  field  in  ducts  may  in  general  by  regarded  in  the  following, 
two  different  ways:  in  the  first  interpretation,  it  can  be  visualized  as 
the  superposition  of  two  sets  of  free-space  waves,  one  propagating  in  the 
positive  y direction  and  the  other  in  the  negative  y direction;  alternatively, 
it  can  be  viewed  that  these  waves  are  continually  reflected  back  and  forth 
between  two  rigid  walls.  Let  us  at  first  adopt  the  first  viewpoint  and 
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consider  these  two  waves  generated  by  the  initial  pressure  specified 
at  x - 0;  for  the  n-th  mode  in  uniform  flow,  the  wavelength  in  the 
y direction  is  equal  to  2£/n,  (shown  in  Figure  4(a)  for  n • 1 ) with  the 
pressure  oscillating  at  frequency, u . The  wave  constant,  ot , in  the  x 
direction  is  proportional  to  the  number  of  crests  (or  troughs)  in  unit 
length  in  that  direction.  Let  us  pose  the  following  question:  what  is 
the  change  in  the  number  of  crests  due  to  the  non-uniform  velocity 
distribution?  In  order  to  answer  this,  we  extend  the  shear  profile, 
originally  defined  inside  the  duct,  to  its  outside  in  periodic,  repeated 
manner  in  order  to  cover  the  entire  unbounded  space,  as  shown  in 
Figure  4(b) . (For  odd  component  of  velocity  distribution,  the  profile  in 

the  neighborhood  of  y ■ + & , + 2&,  may,  if  necessary,  be  envisioned 

to  be  slightly  modified  so  that  the  velocity  distribution  becomes  continu- 
ously smooth;  this  is  to  avoid  the  unnecessary  complication  arising  from  the 
refraction  which  would  otherwise  occur  at  the  interface  between  velocity 
discontinuity.)  Ue  note  that,  for  any  modes,  the  lines  of  constant  phase 
in  the  free  space  may  be  given  by  the  curves  along  which  the  imaginary 
part  of  (3.1)  remains  constant,  i.e., 

a(oMQ  - k) 
ox  ± Yy  ± e 

-CO 


♦ ~ 2 o(oM  - k)]  ^ g'(y)  - g’"(y)  + 

M0  k ° (2y)3  (2y)5  (2y)7 


* (y)  ] " constant. 


(4.12) 


We  observe  that  because  of  the  periodic  extension  of  the  shear  profile 
throughout  the  entire  space  and  with  the  aid  of  (4.3),  the  Integral 
Involved  may  be  reduced  to  the  following: 
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g(n)  dn  - I g(n)  dn 

-i 


We  at  once  recognize  that  at  the  points  corresponding  to  walls,  y = + £, 

the  integral  vanishes.  Furthermore,  for  odd  component  of  velocity 

distribution,  the  odd  derivatives  take  the  same  value  at  these  points: 

g'(£)  * g'(-&),  g' ' ' (£)  * g"'(-it)  and  so  on.  Consequently,  these  promptly 

lead  to  the  following:  in  (4.12),  all  the  terms  associated  with  e or 

shear  take  the  same  value  at  the  points  corresponding  to  walls.  Therefore 

lines  of  constant  phase  appear,  as  shown  in  Figure  (4c),  where, 

for  the  sake  of  illustration,  the  lines  are  adjusted  so  that  the  lines 

of  constant  phase  for  uniform  and  for  non-uniform  flow  cross  each  other 

initially  at  y * - H.  When  we  now  switch  to  the  second  viewpoint  of 

regarding  the  duct  wave  system  as  being  reflected  repeatedly  back  and 

forth  between  the  walls,  one  can  immediately  recognize  that  for  odd  component 

of  velocity  distribution,  the  number  of  wave  crests  remains  unaffected 

by  the  presence  of  shear; and  therefore,  in  such  a case  y of  (4.11)  vanishes. 

n 

(It  is  apparent  that  this  does  not  hold  for  even  component  of  velocity.) 

A similar  argument  can  also  be  used  to  provide  the  physical  explanation  of 
Xq  ■ 0 for  n * 0 or  the  lowest  mode,  for  in  the  corresponding  free-space 
solution  of  (2.11c),  the  shear  effect  appears  only  as  the  real  part  of 
the  argument  of  the  exponential  function,  leaving  the  phase  surface  unaltered. 
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For  general,  arbitrary  shear  profile,  the  solution  for  higher 
mode  may  be  obtained  by  enforcing  the  boundary  condition  (4.1),  upon 
the  linear  combination  of  the  solutions,  which  can  be  obtained,  in  a manner 
similar  to  the  one  used  for  free  space;  in  the  course  of  this  we  naturally  have 
to  keep  in  mind  the  correction  of  aR  as  expressed  in  (4.11).  The  final  solution 
for  the  n-th  mode  is  as  follows: 

Pn'  ■ y (exp  zn  + exp  zfl)  e^^n  + ^n^  x iu)t,  n = 1,  2 (4.13a) 


where  zq  is  the  complex  conjugate  of  z^  defined  by 

2a  y 

zn  - iYny  + I 8(n)  dn 

n 0 


+ e [ 


4a  Y 


n n 


- 2a  (a  M.  - k) ] 


a Mn  - k n v n 0 

n 0 


x [ / dn  / g(0  e"2iYn(n  “ Qd£ 
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and  where 


a y , 

n n 1 


*n  ‘ MM^oMq-  k)  - an]  IotnM0  - k " 2 an(anM0  _ k)  ] 


x r — [g», 


-5-  [g'a)  - g'  (-i)]  - ~ t [g'"w  - g,,,(-n)]+ 

(2YJ2  (2yJ4 


(4.13c) 


It  is  quite  apparent  that  for  odd  velocity  distribution,  X * 0,  as  we 

expected  before.  Needless  to  say,  in  the  above,  as  in  the  free  space 

+ 2iY  v 

solution,  the  parasitic  terms  of  e—  n , which  may  arise  in  the  course 
of  carrying  out  the  Integration  for  specific  velocity  distribution,  have 
again  to  be  discarded. 

Although  the  above  results  are  obtained  for  hard-walled  ducts, 
one  can  obtain,  in  similar  manner,  the  expressions  for  the  ones  with  soft 
or  treated  wall  as  well. 


> 


L 
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5.  Concluding  Remark. 


In  summary,  it  has  been  our  intent  to  demonstrate  that  the  present 
method  provides  a simple  and  effective  means  to  obtain  a solution  for 
waves  propagating  through  aribtrarily  sheared  flow.  Our  main  results, 
accrued  in  closed  form  for  both  unbounded  media  (2.11)  and  ducted  flow 
(4.13),  involve  only  quadratures  of  any  given  shear  profile;  as  special 
cases,  they  are  found  to  embrace  other  known  solutions  and  besides,  they 
can  be  rendered  into  forms  particularly  convenient  for  the  extraction 
of  various  physical  features.  The  method  appears  to  be  easily  applicable 
to  other  problems  involving  disturbance  propagation  through  non-uniform 
media.  (For  example,  the  method  is  now  being  found  to  be  particularly 
instrumental  in  assessing  the  effect  of  spanwise  flow  variation  upon 
flutter  in  turbomachinery  bladings.) 

The  author  owes  his  sincere  gratitude  to  his  former  colleagues, 

Drs.  Shankar  and  Savkar  for  arousing  his  early  interest  to  the  problem. 

The  work  was  in  part  supported  by  the  Air  Force  Office  of  Scientific 
Research  under  Contract  No.  F 44620- 74-C-0040. 
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APPENDIX  A 


We  shall  show  that,  under  the  assumption  of  small  velocity 
gradient,  (2.14)  expressed  in  terms  of  the  parabolic  cylinder  function 
becomes  equal  to  (2.13).  For  this  we  make  use  of  the  following  formula 
of  the  parabolic  cylinder  function,  due  to  Darwin  (1949): 


U(b,  + O * Y 

where 


: f 

exp  ( + i ($  + 

1 

— 2 2 

Y - (x  - 4a;  , 

1 _ 

<f>  = 4 x Y - a log 


+ -2—)  + 
-15'' 

Y 


(^  + 
x 


x + Y 
2 /a". 


(A- lb) 
(A-lc) 


a - -ib,  (A-ld) 

1 

x = i2  C,  (A-le) 


and  the  expression  for  g,  £,.  . . etc.  are  given  in  Darwin's  paper. 

Since  in  the  present  case  of  Example  2 of  Section  2,  g(y)  of  (2.6)  is  equal 
to  y,  it  follows 


dM 

dy 


e. 


and  from  (2.15b),  (A-ld)  and  (A-le)  become 


f 


(A-2a) 


! 
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i (^)2  (k  - cM0  - aey)  . 


(A-2b) 


Into  (A-lc),  we  substitute  the  above  and,  after  expanding  and  discarding 
a certain  constant  which  is  disposable,  we  obtain 


i 2 k “ gMq  2 

<J>  - yy  - T aey  — + 0 (e  ) , 


(A-3) 


where  Y is  defined  in  (2.8b).  Likewise 


22  aey<k  “ M0>  14221 

y = iy<£)2  n ^ + 


] 


(A-4) 


By  substituting  (A-3)  and  (A-4)  into  (A-la) , it  is  found  that  the 
terms  associated  with  g,  , . . . etc.  become  negligible  and  we  obtain 

— 1 2 k ~ °^o  i ~ k 

U(b,  + o * exp  (+  iyy)  exp  [+  j ictey  j aey ]. 

(A-5) 

Upon  substituting  the  above  into  (2.14)  and  noting  that  the  factor 
b£  associated  with  the  second  term  within  the  bracket  can  be  expressed  as 


b£  - exp  [log  (bC)] 


— (g)2  (k  - OM0)  exp  \ log  (1  - 


2lT  # <k  - **n>  S exP  -aSZ—  + 0 <e2>J  > * 


k - 0Mo 


it  follows  that,  due  to  this  factor,  the  second  term  dominates  for 
small  e;  from  this,  (2.14)  takes  the  same  form  as  (2.13). 
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APPENDIX  B 

Here,  we  shall  find  that  for  small  velocity  discontinuity, 
V-.17)  becomes  (2.16).  From  (2.17d),  we  obtain 


2a 
COS  <P1 

<(>,  - 4).  - eAM  ' — ; . 

2 1 sin  <p. 


(B-l) 


From  (2.17c)  and  the  definition  of  Y given  in  (2.8b),  a tan  <J>2  appearing 
in  (2.17a)  becomes 


AM 


a tan<f>2  - Y ± e a (aM0  - k)  — . 


(B-2) 


With  the  aid  of  (B-l),  T of  (2.17b)  becomes 


T - 1 + eAM 


a[2Y2  - (aM0  - k)2] 

2y2(oM0  - k) 


(B-3) 


which,  to  the  order  of  t , is  equal  to  (2.16b).  Thus,  (2.17)  becomes 
(2.16),  to  the  present  order  of  approximation. 


1 
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APPENDIX  C 


We  prove  that  (3.6a)  is  identical  to  (3.6b).  For  this,  we 
rewrite  (3.6a)  as 


/2rf 


ei0Cc  1(a)  da  , 


where 


I (a)  = — — f e-iut  e±iYy  dw  . 

/2i 


From  the  definition  of  Y given  by  (3.5b),  (C-lb)  becomes 


(C-la) 


(C-lb) 


1(a) 


/ 2TT 


-iait 

e exp 


+ i [(u  - ac(M  - l))(aj  - ac(M  + 1))]  > du. 

— c 


(C-2) 


Now  from  Campbell  and  Foster  (1961,  p.  109,  No.  860.0),  one  notes 


- f _ iut 

[(u  - c(M  + 1)X0J  - c(M  - 1))]  * e 


x exp  ( t i c “ ac(M  + 1))]*  ( du 


2iri  H(t  + c y)  exp  (-i  ca  Mt)  Jq  [ca  (t2  — y2)2  ] , 

c 


(C-3) 
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where  H Is  the  step  function  and  Jg  Is  the  Bessel  function*  Differentiating 
both  sides  with  respect  to  y,  we  obtain 


I (a)  - + /2tF  • c exp  (-icaMt) 


tt-  l H(t  + i y)  Jn  [ca  (t2  - -4  y2)  ] 


3y 


(C-4) 


Into  (C-la) , we  substitute  (C-4) ; with  the  aid  of  the  following  identity 


J0(b£)  eiCx  dt  - 2(b2  - x2)  2H  (b  - |x|), 

which  Is  found  in  Erdelyi,  et  al  (1954,  vol.  1,  p.  43),  and  with  the 
definition  of  U given  in  (3.7),  we  obtain  (3.6b). 
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TABLE  1 


COMPARISON  OF  TRANSMISSION  COEFFICIENT 


am  = 0.05 


incidence  angle 
^(degrees) 

transmission  coefficient 

T (Miles-Ribner) 

transmission  coefficient 

A (present) 

90 

1.000 

1.000 

80 

.996 

.996 

70 

.992 

.993 

60 

.992 

.992 

50 

.996 

.995 

40 

1.011 

1.008 

30 

1.059 

1.044 

20 

1.324 

1.166 

17.75* 

2.000 

1.232 

*critical  Incidence 


AM  - -0.05 


incidence  angle 
^(degrees) 

transmission  coefficient 

T (Miles-Ribner) 

transmission  coefficient 

A (present) 

90 

1.000 

1.000 

80 

1.004 

1.004 

70 

1.007 

1.007 

60 

1.008 

1.008 

50 

1.005 

1.005 

40 

.994 

.992 

30 

.967 

.958 

20 

.896 

.857 

10 

.691 

.464 

0 

0 

0 

I 
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Examples  of  Shear  Profile 
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Coordinate  In  Ducts. 
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-Dimensional,  Refraction  Effect  upon  Unsteady 
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Three-Dimensional,  Refraction  Effect  Upon  Unsteady 
Oscillating  Airfoils  in  Supersonic  Flow 

M.  Kurosaka* 

The  University  of  Tennessee  Space  Institute 
Tullahoma,  Tennessee  37388 

I.  Introduction. 

On  the  subject  o£  three-dimensional  effect  upon  the  supersonic 
flutter  in  aircraft  engines,  Ref.  1 identified  the  following,  two  key 
issues  which  could  significantly  modify  the  results  based  upon  two- 
dimensional  approximation:  (1)  reflection  of  waves  from  the  encasing 
walls  and  (2)  refraction  arising  from  the  propagation  of  disturbances 
through  steady  flow  with  strong  velocity  gradient  in  the  radial  direction 
existing  in  turbomachines.  Of  these  two.  Ref.  1 specifically  restricted 
itself  to  the  first  problem  or  its  variant  of  practical  significance  where 
the  sound  absorbent  material  is  Installed  upon  the  wall  surface. 

In  the  present  paper,  we  turn  our  attention  towards  the  second 
one  and  examine  to  what  extent  the  unsteady  pressure  on  the  airfoils  will 
be  affected  by  the  non-uniformity  of  surrounding,  steady  flow.  As  is 
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Contract  No.  F 44620-74-C-0040. 
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well  known,  the  sheared  velocity  distribution  always  refracts  the  sound 
waves;  for  example,  according  to  a geometrical  acoustics  picture  convenient 
for  high  frequencies,  the  rays  of  sound  will  be  bent  while  propagating 
through  any  stratified  media.  The  refraction  of  disturbances  traversing 
a flow  field  with  steep  radial  velocity  gradient  present  in  turbomachinery 
environment  may  considerably  alter  the  unsteady  pressure  signals,  when 
compared  to  analysis  based  upon  two-dimensional  approximations.  In  order 
to  focus  our  attention  solely  on  the  refractive  aspect  of  the  problem, 
we  shall  attempt  to  extricate  ourselves  from  the  other  unneccessary 
complications.  Therefore  we  pose  and  study  a model  problem,  which  seems 
to  capture  the  central  nature  of  the  problem,  and  consider  an  isolated 
oscillating  airfoil  which  is  placed  in  a supersonic  duct  with  non-uniform 
velocity.  The  steady,  sheared  velocity  is  in  the  postive  x direction  (Figure 
1)  and  is  function  of  y only,  i.e.  U(y).  The  walls  of  the  walls  located 
at  y • H and  y ” -A  are  rigid  and  the  airfoil  is  harmonically  oscillating 
with  small  amplitude  in  the  z direction. 

Our  aim  is  to  compare  the  effect  of  sheared  flow  upon  the  unsteady 
pressure  acting  on  the  surface  of  the  airfoil  with  the  uniform  flow.  In 
analyzing  the  problem,  we  rely  exclusively  on  the  results  of  Ref.  2,  where 
a simplified  representation  of  unsteady  disturbance  propagating  through 
non-uniform  flow  is  obtained. 


s: 

t 
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II.  Formulation. 


The  linearized  governing  equations  can  be  written 


+ 3r> 


+ p (• 


3u' 

3x 


hL 

3y 


*1)  - 
3z  ' 


3u' 

3t 


+ U(y) 


3u* 

3y 


+ v.  Mzl 

dy 


1 3p 1 
p 3x  * 


3v’ 

3t 


+ £r 


1 3£l 

P ay  ’ 


3w* 

3t 


+ u(y) 


3w* 

3x 


1 3p' 

p 3z 


Cl) 


where  p and  a are  constant  ambient  density  and  acoustic  speed,  respectively 
and  primes  denote  perturbed  quantities  with  conventionally  defined  meanings. 
As  usual,  in  the  above  we  write  all  the  perturbed  quantities  as  the  product 
of  two  factors,  one  for  the  spatially  dependent  amplitude  and  the  other 
for  time-dependent  factor  such  as  p'  ■ p'  eia)t  and  obtain 


+ U(y) 


iEl 

3x 


) + PC- 


+ ill  3wl 

3x  + 3y  + 3z  ' U> 


iu  u'  + UCy)  + v'  ^ 

iw  v*  + UCy) 

UCy)  |jp 


1 3p ' 

* p 3x  ’ 

-.1  iEl 

p 3y  * 

1 3p' 

p 3* 


C2) 
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The  boundary  condition  on  walls  Is  that  the  normal  velocities  are  equal 
to  zero,  i.e., 

v'(x*y  ■ +&  , z)  - 0 . (3) 

On  the  airfoil,  the  velocity  normal  to  it  is  prescribed  by  the  given 
movement  of  the  airfoil  and  if  its  amplitude  is  equal  to  W(x,y),  then 

w'(x,y,  z - 0)  - W (x,y)  . (4) 

Upstream  of  the  initial  Mach  cone,  the  flow  is  quiescent.  Because  of  the 
supersonic  nature  of  the  flow,  the  flow  fields  above  and  below  the  airfoil 
is  unrelated  and  hereafter  we  consider  only  the  flow  above  the  airfoil, 
z > 0. 


1 

t 


5 


III.  Solution. 

In  order  to  obtain  a solution  of  (2),  we  first  eliminate  the  x 
dependence  by  taking  the  Fourier  transform  in  that  direction;  next,  the 
application  of  the  Fourier  cosine  transform  In  the  z direction  removes  the 
z dependence  and  the  resultant  equation  involving  only  y is  solved  with 
the  aid  of  the  results  obtained  in  Ref.  2;  the  description  of  these 
successive  steps  will  be  given  below. 

A.  Fourier  transform  in  the  x direction. 

We  take  the  Fourier  transform  of  (2)  defined,  for  example,  by 


<30 


-00 


and  eliminating  the  resulting  u,  v and  w,  one  obtains  the  following 
equations  in  terms  of  p: 


. 3ji  2a  dM(y) 

3z2  oM(y)  - k dy 


§|  + t(cM(y)  - k)2 


a2] 


0 , 


where  k “ 


— and  with  the  boundary  conditions 


(5a) 


- -ip (to  - U(x)a)  W (y)  . 

0 


(5b) 


(5c) 
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B.  Fourier  cosine  transform  in  the  z direction. 


Taking  the  Fourier  cosine  transform  of  p defined  by 


00 

pc  (?>  "Vr  f p 


cos  0 zdz , 


(5)  becomes 


d2Fc(p) 


P)  jy  dF  (P)  _ 0 0 

— - s?frnr  sjf  + «*<>■>  - k>  - <»  + * » Fc  ® 


-ip*  [cM(y>  - • w (y)  > 


with  the  boundary  condition 


dF  (p) 

— § - 0 . 

dy 


y - + £ 


In  order  to  solve  this  Inhomogeneous  equation,  we  first  consider,  according 
to  the  standard  method,  eigenfunctions  of  a corresponding  homogeneous  equation 
given  by 


d2F  2a  dM  dF  . 2 * 2 . „2*,  „ 

7T  - c5(y)“k  57  37  + "“"(r)  ~ k>  - (<*  + e )J  ? - o . 


with  the  boundary  condition 


■ 0 . 


hr  - ± £ 
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To  solve  (7),  one  assumes  that  the  non-uniformity  of. the  flow  inside 
the  duct  is  small , i . e . , 

M(y)  - Mq  + eg(y), 

where  Mq  is  the  averaged  velocity  in  the  duct  defined  by 


(8a) 


Mq  - j H(y)  dy. 


(8b) 


-l 

and  C « 1 . According  to  (8a)  and  (8b),  it  immediately  follows  that 

] 

. g(y)  dy  - 0 . 

-I 


(8c) 


Under  this  assumption,  the  eigenfunctions  of  (7)  are  derived  in  Ref.  2 and 
they  are  given  by 


uQ  (a,  y)  - exp  [-2  ea(oMQ  - k)  J dr}  J g(0  d£  ], 

0 -l 


(9a) 


un  *a’y*  “ 2"  *exp  Zn  + ®Xp  zn* ’ n * lf  2’ 


(9b) 


where  z la  the  complex  conjugate  of  z defined  by 
n “ 


-i  r.r  * lc 


oM„  - h 


g(n)  du 


-i 


4ay 


+ £ tgj— ^ - 20  <*o  - W1 


V 


8 


x [ 


I dn  I g(C)  e"2iY*  (n  “ C)  dC  + — ^ g'(*)  - — g'"<£)  + ] 

J j <2 * *V  (2Yn)5 


-a  -a 


ie 


+ 71  Xn  (y  - A)  ^(aMo  - k)  - a]. 


(9c) 


and  where 


rnr 


Yn  - X * n " !.  2» 


(9d) 


a Y 2 


*n  " AIMqCoMq 


n - k)  - o]  [oMn  - k " 2 (aM0  " k) 1 


x < — i-5-  Ig'(a)  - g'(-a)] — , [g,"(a)  - g”'(-a)]  + — \, 

(2V  «v  ' 


(9e) 


2 2 2 

and  (oMq  - k)  - a - 3 takes  the  following  eigenvalues: 


(c<M0  - k)2  - a2  - g2  - (^)  + 2X  ne  [M0(aM0  - k)  - a]  . 


(9f) 
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(Contrary  to  the  situation  of  Ref.  2,  where  the  eigenvalues  are  assigned  only 
to  a,  here  we  ascribe  it  to  the  constant  consisting  of  combination  of 
a,  0,  Mq  and  k,  which  appears  on  the  left  hand  side.)  Since  (7)  is  a 
Sturm-Llouville  equation,  the  eigenfunctions  uq  are  orthogonal  and  they 
are  Indeed  complete.  Hence  we  expand  ?c(p)  of  the  inhomogeneous  equation 


of  (6a)  and  also  its  right  hand  side  into  series  of  eigenfunctions,  l.e. , 


F (p) 


■ L 


a u (a  , y)  , 
n n 


n * 0 


and 


W (y)  - ^ 


bn  (a  , 6)  ur  (a  , y). 


n * 0 


g(y)  tf  (y)  - ^ cn  (a  , 3)  un  (a  , y). 


n • 0 
£ 


and  where 


W(y)  uQ(a,  S,  y)  dy 


bn  (a  , p)  - 


-£ 


[uQ(a,  3,y)l  dy 


-l 


and 


g(y)  w(y)  un(a,  3,y)  dy 


£n(“’  W 


j [un(a,  6,y)]2  dy 
-£ 


Substituting  these  Into  <6a)  and  making  use  of  the  fact  that  uq 
satisfies  the  homogeneous  equation  (7a) , one  can  determine  an  of  (10a) 
and  F (p)  becomes 


rr 
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■vr 


r r (cMo  - k)  go  . ae£o 
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(OM0  - k)  - a2  - e2  - (Sf>  - 2xne[M0(cM0  - k)  -a] 


+ ae 


z „ , .2  2 02  ,wrN2 

(aM0  - k)  - a - B - (-y) 


un(a  , B). 


(ID 


If  the  velocity  on  the  surface  of  the  airfoil,  (4)  is  such  that 


W(x,y)  - H(x)X(x)  (y  + l) 


21 


(12) 


where  H la  a step  function  and  $ is  the  amplitude  of  oscillation  and  for 
linear  duct  flow  distribution  represented  by 

g(y)  - Ay  , 

then  Xq"  0 and  (lOd)  and  (lOe)  become 


(13) 


n 


Inverting  (11),  we  obtain 


One  can  show  that  if  one  replaces  i- sine  of  (15)  by  cosine,  it  also  satisfies 
(5a  and  b);  and  (5c)  becomes,  in  such  a case 


12 


|z  ■ 0 


Hence  the  solution  Is  a sum  of  these  two  expressions,  i.e. 


p - pa  [(OMq  - k)  Bq  + ae  cQ] 


' (ciMq  - k;  - a 


u0  (a,  y) 


x exp  [lz V(aMn  - k)2  - o2] 


ao 
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'(OM0  - k)2  - a2 


(i.} 


« («.  y) 


x exp  [iz  V(aM0  - k)2  - a2  - (*f)2]  . 


C.  Inversion  of  Fourier  Transform. 


Inverting  (16)  with  respect  to  a and  by  deforming  the  contour 
of  integration,  we  obtain  the  pressure  on  the  airfoil,  z ■ 0,  as 


p'  - ( X (x)  [a^O)  + a2'(0)  + a3  (0)) 


-X'(x)  [a*  (0)  + 2a3*(0)] 


+X"(x)  a3  (0) 


* 1 

‘ j X <*)  ag(x  “ t)  dt  1 , 

0 J 


(17a) 
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a3(x)  - - 5i  0Q(x) 


CO 

Y.  21  v e»w  ■ 


(17e) 


- 1 l 


where 


a 4 4 £ 3 

y - MqJI  + -j-j  At  ek  + -j  A£  , 


v « AX4  ek  MQ  , 


A A 9 

c - I AfeM0  , 


and  where 


-2i  V 

0 (x)  - — p=r  exp  (-ix  — j ) 

n mY2rr  m 


x / cos  h sln0)  exp  (~iX^n  co»0)  VQ(®»  y)  d8  » 


(17f) 


2 1 

and  ^ 1 » 

m in 


vn(8,y)  - un  (a,  y)  , 


which  is  given  in  (9)  and  where 

V 

a ■ X cos  0 + — • 
a > 


t 


If  the  amount  of  shear,  £,  is  equal  to  zero,  one  can  readily 
show  that  the  above  expression  is  reduced  to  the  standard  solution  of 
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two-dimensional  airfoil  oscillating  in  uniform  supersonic  flow. 


16 


1 


TV.  Discussion  of  Results.  * f 

— — 2 I 

Figure  2 and  3 compare  c ■ p'/2  Pu  , computed  from  (17)  and  ; 

P ! 

evaluated  on  the  surface  of  airfoils  placed  in  a linearly  sheared  duct  l 

with  the  one  corresponding  to  uniform  flow.  The  pressure  is  calculated  j 

j 

at  the  top  of  the  duct  and  the  airfoil  is  oscillating  in  torsional  motion  1 

\ 

at  mid- chord  with  unit  amplitude.  Additional  parameters  are  as  follows:  j 

\ 

l j 

— ■ 1.5,  where  JL  is  half-duct  height,  c is  airfoil  chord,  j 


and 

M - 1.3  + 0.0666  (J)  . 

It  clear  that  the  spanwlse  velocity  distribution  corresponding  to  the 
three  dimensional,  radial  velocity  gradient  indeed  induces  significant 
modification  to  the  one  for  radially  uniform  or  two-dimensional  approximation. 


r 


17 


References. 


^Kurosaka,  M.  and  Edelfelt , I.  H. , "Suppressive  Effect  of  Liners  in 
Supersonic  Compressor  Flutter,"  to  be  submitted  to  AIAA.  J. 

^Kurosaka.  M. . "Towards  Simplified  Representation  of  Wave  Propagation 
through  Non-Uniform  Media." 


n |x 


